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ABSTRACT 

Gamma-ray bursts (GRBs) have puzzled astronomers since their accidental discovery in the late sixties. 
The BATSE detector on the COMPTON-GRO satellite has been detecting one burst per day for the last 
six years. Its findings have revolutionized our ideas about the nature of these objects. They have shown 
that GRBs are at cosmological distances. This idea was accepted with difficulties at first. The recent 
discovery of an X-ray afterglow by the Italian/Dutch satellite BeppoSAX has led to a detection of high 
red-shift absorption lines in the optical afterglow of GRB970508 and in several other bursts and to the 
identification of host galaxies to others. This has confirmed the cosmological origin. Cosmological GRBs 
release ~ 10 51 — 10 53 ergs in a few seconds making them the most (electromagnetically) luminous objects 
in the Universe. The simplest, most conventional, and practically inevitable, interpretation of these 
observations is that GRBs result from the conversion of the kinetic energy of ultra-relativistic particles 
or possibly the electromagnetic energy of a Poynting flux to radiation in an optically thin region. This 
generic "fireball" model has also been confirmed by the afterglow observations. The "inner engine" that 
accelerates the relativistic flow is hidden from direct observations. Consequently it is difficult to infer 
its structure directly from current observations. Recent studies show, however, that this "inner engine" 
is responsible for the complicated temporal structure observed in GRBs. This temporal structure and 
energy considerations indicates that the "inner engine" is associated with the formation of a compact 
object - most likely a black hole. 



1. INTRODUCTION 

Gamma-ray bursts (GRBs), short and intense bursts of 
~ lOOkeV-lMeV photons, were discovered accidentally in 
the late sixties by the Vela satellites Q . The mission of 
these satellites was to monitor the "Outer Space Treaty" 
that forbade nuclear explosions in space. A wonderful by- 
product of this effort was the discovery of GRBs. 

The discovery of GRBs was announced in 1973 It 
was confirmed quickly by Russian observations 0] and by 
observations on the IMP-6 satellite Since then, sev- 
eral dedicated satellites have been launched to observe the 
bursts and numerous theories were put forward to explain 
their origin. Claims of observations of cyclotron spectral 
lines and of discovery of optical archival counterparts led in 
the mid eighties to a consensus that GRBs originate from 
Galactic neutron stars. This model was accepted quite 
generally and was even discussed in graduate textbooks 
5, |(| and encyclopedia articles [M |[. 
lie BATSE detector on the COMPTON-GRO 
(Gamma-Ray Observatory) was launched in the spring of 
1991. It has revolutionized GRB observations and con- 
sequently our basic ideas on their nature. BATSE ob- 
servations of the isotropy of GRB directions, combined 
with the deficiency of faint GRBs, ruled out the galac- 
tic disk neutron star model 2 and make a convincing case 
for their extra-galactic origin at cosmological distances . 
This conclusion was recently confirmed by the discovery by 
BeppoSAX [pi of an X-ray transient counterparts to sev- 
eral GRBs. This was followed by a discovery of optical 



|ll , 12l and radio transients [|T3|. Absorption line with a 
rmshrft z = 0.835 were measured in the optical spectrum 
of the counterpart to GRB970508 pi] providing the first 
redshift of the optical transient andthe associated GRB. 
Latter, redshiftcd emission lines from galaxies associated 
with GRB971214 p| (with z = 3.418) and GRB980703 
p| (with z = 0.96H) were discovered. Galaxies has been 
discovered at the positions of other bursts. There is lit- 
tle doubt now that some, and most likely all GRBs are 
cosmological. 

The cosmological origin of GRBs immediately implies 
that GRB sources are much more luminous than previ- 
ously thought. They release ~ 10 51 — 10 53 ergs or more 
in a few seconds, the most (electromagnetically) luminous 
objects in the Universe. This also implies that GRBs are 
rare events. BATSE observes on average one burst per 
day. This corresponds, with the simplest model (assuming 
that the rate of GRBs does not change with cosmologi- 
cal time) to one burst per million years per galaxy. The 
average rate changes, of course, if we allow beaming or a 
cosmic evolution of the rate of GRBs. 

In spite of those discoveries, the origin of GRBs is still 
mysterious. This makes GRBs a unique phenomenon 
in modern astronomy. While pulsars, quasars and X- 
ray sources were all explained within a few years, if not 
months, after their discovery, the origin of GRBs remains 
unknown after more than thirty years. The fact that GRBs 
are a short transient phenomenon which until recently did 
not have any known counterpart, is probably the main rea- 
son for this situation. Our inability to resolve this riddle 
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2 A few GRBs, now called soft gamma repeaters, compose a different phenomenon, are believed to form on galactic neutron stars. 
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also reflects the accidental and unexpected nature of this 
discovery which was not done by an astronomical mission. 
Theoretical astrophysics was not ripe to cope with GRBs 
when they were discovered. 

A generic scheme of a cosmological GRB model has 
emerged in the last few years and most of this review is 
devoted to an exposition of this scheme. The recently 
observed X-ray, optical and radio counterparts were pre- 
dicted by this picture @ £|, EH- EH H • This discovery 
can, to some extent, be considered as a confirmation of 
this model |^2|, [2^, |2^| . According to this scheme the 
observed 7-rays are emitted when an ultra-relativistic en- 
ergy flow is converted to radiation. Possible forms of the 
energy flow are kinetic energy of ultra-relativistic particles 
or electromagnetic Poynting flux. This energy is converted 
to radiation in an optically thin region, as the observed 
bursts are not thermal. It has been suggested that the en- 
ergy conversion occurs either due to the interaction with 
an external medium, like the ISM [ p7| or due to internal 
process, such as internal shocks and collisions within the 
flow ||, H, |§. Recent work @ |lj shows that the ex- 
ternal shock scenario is quite unlikely, unless the energy 
flow is confined to an extremely narrow beam, or else the 
process is highly inefficient. The only alternative is that 
the burst is produced by internal shocks. 

The "inner engine" that produces the relativistic energy 
flow is hidden from direct observations. However, the ob- 
served temporal structure reflects directly this "engine's" 
activity. This model requires a compact internal "engine" 
that produces a wind - a long energy flow (long compared 
to the size of the "engine" itself) - rather than an explosive 
"engine" that produces a fireball whose size is compara- 
ble to the size of the "engine" . Not all the energy of the 
relativistic shell can be converted to radiation (or even 
to thermal energy) by internal shocks |}2], [33|, The 
remaining kinetic energy will most likely dissipate via ex- 
ternal shocks that will produce an "afterglow" in different 
wavelength [^0|. This afterglow was recently discovered, 
confirming the fireball picture. 

At present there is no agreement on the nature of the 
"engine" - even though binary neutron star mergers [p5[ 
are a promising candidate. All that can be said with some 
certainty is that whatever drives a GRB must satisfy the 
following general features: It produces an extremely rela- 
tivistic energy flow containing w 10 51 — 10 52 ergs. The flow 
is highly variable as most bursts have a variable temporal 
structure and it should last for the duration of the burst 
(typically a few dozen seconds). It may continue at a lower 
level on a time scale of a day or so Bq] . Finally, it should 
be a rare event occurring about once per million years in 
a galaxy. The rate is of course higher and the energy is 
lower if there is a significant beaming of the gamma-ray 
emission. In any case the overall GRB emission in 7-rays 
is ~ 10 52 ergs/10 6 years/galaxy. 



We begin (section 
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with a brief review of GRB ob- 
H, |C], H for additional reviews 
45| for a more extensive discussion). We 
then turn to an analysis of the observational constraints. 
We analyze the peak intensity distribution and show how 
the distance to GRBs can be estimated from this data. 
We also discuss the evidence for another cosmological ef- 



fect: time-dilation (section [|). We then turn (section U) 
to discuss the optical depth or the compactness problem. 
We argue that the only way to overcome this problem is 
if the sources are moving at an ultra-relativistic velocity 
towards us. An essential ingredient of this model is the no- 
tion of a fireball - an optically thick relativistic expanding 
electron-positron and photon plasma (for a different model 
see however j46|). We discuss fireball evolution in section 
|6| Kinematic considerations which determine the observed 
time scales from emission emerging from a relativistic flow 
provides important clues on the location of the energy con- 
version process. We discuss these constraints in section 
and the energy conversion stage in sectionH. We review 
the recent theories of afterglow formation ^We examine 
the confrontation of these models with observations and 
we discuss some of the quantitative problems. 

We then turn to the "inner engine" and review the re- 
cent suggestions for cosmological models (section [k]) . As 
this inner engine is hidden from direct observation, it is 
clear that there are only a few direct constraint that can 
be put on it. Among GRB models, binary neutron star 
merger p5j is unique. It is the only model that is based 
on an independently observed phenomenon j48fl, is capable 
of releasing the required amounts of energy~~G9] within a 
very short time scale and takes place at approximately the 
same rate [|o[ |51 , 3 . At present it is not clear if this 
merger can actua ly channel the required energy into a rel- 
ativistic flow or if it could produce the very high energy 
observed in GRB971214. However, in view of the special 
status of this model we discuss its features and the pos- 
sible observational confirmation of this model in section 

GRBs might have important implications to other 
branches of astronomy. Relation of GRBs to other as- 
tronomical phenomena such as UCHERs, neutrinos and 
gravitational radiation are discussed in section The 
universe and our Galaxy are optically thin to low energy 
7-rays. Thus, GRBs constitute a unique cosmological pop- 
ulation that is observed practically uniformly on the sky 
(there are small known biases due to CGRO's observation 
schedule). Most of these objects are located at z « 1 or 
greater. Thus this population is farther than any other 
systematic sample (QSOs are at larger distances but they 
suffer from numerous selection effects and there is no all 
sky QSOs catalog). GRBs are, therefore, an ideal tool to 
explore the Universe. Already in 1986 Paczynski [^3| pro- 
posed that GRBs might be gravitationally lensed. This 
has led to the suggestion to employ the statistics of lensed 
bursts to probe the nature of the lcnsing objects and the 
dark matter of the Universe [gij . The fact that no lensed 
bursts where detected so far is sufficient to rule out a criti- 
cal density of 1O 6 5 M to 10 81 M Q black holes ||. Alter- 
natively we may use the peak-flux distribution to estimate 
cosmological parameters such as Q and A [[36| . The angu- 
lar distribution of GRBs can be used to determine the very 
large scale structure of the Universe |37[ ||^] . The possible 
direct measurements of red-shift to some bursts enhances 
greatly the potential of these attempts. We conclude in 
section |l2] by summarizing these suggestions. 

Over the years several thousand papers concerning 
GRBs have appeared in the literature. With the grow- 
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ing interest in GRBs the number of GRB papers has been 
growing at an accelerated rate recently . It is, of course, 
impossible to summarize or even list all this papers here. 
I refer the interested reader to the complete GRB bibliog- 
raphy that was prepared by K. Hurley |59]. 



2. OBSERVATIONS 

GRBs are short, non-thermal bursts of low energy 7- 
rays. It is quite difficult to summarize their basic features. 
This difficulty stems from the enormous variety displayed 
by the bursts. I will review here some features that I be- 
lieve hold the key to this enigma. I refer the reader to the 
proceedings of the Huntsville GRB meetings Q [O], ff5j 
and to other recent reviews for a more detailed discussion 
fl37|, E)! |3l,|40l|4l. 



2.1. Duration: 

A "typical" GRB (if there is such a thing) lasts about 
lOsec. However, observed durations vary by six orders of 
magnitude, from several milliseconds pOj] to several thou- 
sand seconds plj . About 3% of the bursts are preceded 
by a precursor with a lower peak intensity than the main 
burst [p2[ . Other bursts were followed by low energy X-ray 
tails [ |63| . Several bursts observed by the GINGA detec- 
tor showed significant apparently thermal, X-ray emission 
before and after the main part of the higher energy emis- 
sion m . These are probably pre-discovery detections 
of the X-ray afterglow observed now by BeppoSAX and 
other X-ray detectors. 

The definition of duration is, of course, not unique. 
BATSE's team characterizes it using Tgo (T50) the time 
needed to accumulate from 5% to 95% (from 25% to 75%) 
of the counts in the 50keV - 300keV band. The short- 
est BATSE burst had a duration of 5ms with structure 
on scale of 0.2ms @. The longest so far, GRB940217, 
displayed GeV activity one and a half hours the main 
burst @. The bursts GRB96I027a, GRB961027b, 
GRB96I029a and GRB96I029b occurred from the same 
region in the sky within two days J68[ if this "gang of four" 
is considered as a single very long burst then the longest 
duration so far is two days! These observations may indi- 
cate that some sources display a continued activity (at a 
variable level) over a period of days |6^] . It is also possible 
that the observed afterglow is an indication of a continued 
activity p6| . 

The distribution of burst durations is bimodal. BATSE 
confirmed earlier hints ]70[ | that the burst duration distri- 
bution can be divided into two sub-groups according to 
Tgo: long bursts with Tgo > 2sec and short bursts with 
Tgo < 2sec fn], |7| |7§, |7§ fr| [7§]. The ratio of ob- 
served long bursts to observed short bursts is three to one. 
This does not necessarily mean that there are fewer short 
bursts. BATSE's triggering mechanism makes it less sensi- 
tive to short bursts than to long ones. Consequently short 
bursts are detected to smaller distances |74|, [77[ [78], |79| and 
wc observed a smaller number of short bursts. 



2.2. Temporal Structure and Variability 

The bursts have a complicated and irregular time pro- 
files which vary drastically from one burst to another. Sev- 



eral time profiles, selected from the second BATSE cata- 
log, are shown in Fig. [l| In most bursts, the typical vari- 
ation takes place on a time-scale ST significantly smaller 
than the total duration of the burst, T. In a minority 
of the bursts there is only one peak with no substructure 
and in this case ST ~ T. It turns out that the observed 
variability provides an interesting clue to the nature of 
GRBs. We discuss this in section |?]. We define the ra- 
tio Af = T/ST which is a measure of the variability. Fig. 
[2] depicts the total observed counts (at E > 25keV) from 
GRBI676. The bursts lasted T ~ 100 sec and it had peaks 
of width ST ~ I sec, leading to J\f = 100. 

2.3. Spectrum: 

GRBs are characterized by emission in the few hundred 
keV ranges with a non-thermal spectrum (see Fig. ||) X- 
ray emission is weaker - only a few percent of the energy 
is emitted below lOkeV and prompt emission at lower en- 
ergies has not been observed so far. The current best 
upper limits on such emission are given by LOTIS. For 
GRB970223 LOTIS finds my > II and provides an upper 
limit on the simultaneous optical to gamma-ray fluence 
ratio of < 1.1 x 10~ 4 |8C[ |. Most bursts are accompanied, 
on the other hand, by a high energy tail which contains a 
significant amount of energy - E 2 N(E) is almost a con- 
stant. GRB940217, for example, had a high energy tail up 
to 18 GeVQ. In fact EGRET and COMPTEL (which 
are sensitive to higher energy emission but have a higher 
threshold and a smaller field of view) observations are con- 
sistent with the possibility that all bursts have high energy 

An excellent phenomenological fit for the spectrum was 
introduced by Band et al. p4|: 
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where H = (a — {3)Eq. There is no particular theoret- 
ical model that predicts this spectral shape. Still, this 
function provides an excellent fit to most of the observed 
spectra. It is characterized by two power laws joined 
smoothly at a break energy H . For most observed val- 
ues of a and (5, vF v cx v 2 N(v) peaks at E p — (a + 2)E = 
[{a + 2)/(a - 0)]H. The "typical" energy of the observed 
radiation is E p . That is this is where the source emits the 
bulk of its luminosity. E p defined in this way should not 
be confused with the hardness ratio which is commonly 
used in analyzing BATSE's data, namely the ratio of pho- 
tons observed in channel 3 (100-300keV) to those observed 
in channel 2 (50-100keV). Sometimes we will use a simple 
power law fit to the spectrum: 



N(E)dE oc E~ a dE. 



(2) 



In these cases the power law index will be denoted by a. 
A typical spectra index is a ~ 1.8 — 2 [85[. 

In several cases the spectrum was observed simultane- 
ously by several instruments. Burst 9206022, for example, 
was observed simultaneously by BATSE, COMPTEL and 
Ulysses. The time integrated spectrum on those detectors, 
which ranges from 25keV to lOMeV agrees well with a 
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Band spectrum with: K = 457±30keV, a = -0.86±0.15 
and /3 = -2.5 ± 0.07 (§§. Schaefer ei a/. @ present a 
complete spectrum from 2keV to 500MeV for three bright 
bursts. 

Fig. |i| shows the distribution of observed values of H 
in several samples |84|, |89J . Most of the bursts are 
the range 100 keV < H < 400 keV, with a clear maxi- 
mum in the distribution around H ~ 200keV. There are 
not many soft GRBs - that is, GRBs with peak energy 
in the tens of keV range. This low peak energy cutoff is 
real as soft bursts would have been easily detected by cur- 
rent detectors. However it is not known whether there is a 
real paucity in hard GRBs and there is an upper cutoff to 
the GRB hardness or it just happens that the detection is 
easiest in this (few hundred keV) band. BATSE triggers, 
for example, are based mostly on the count rate between 
50keV and 300keV. BATSE is, therefore, less sensitive to 
harder bursts that emit most of their energy in the MeV 
range. Using BATSE's observation alone one cannot rule 
out the possibility that there is a population of harder 
GRBs that emit equal power in total energy which are not 
observed because of this selection effect [H], [92"[ . 

More generally, a harder burst with the same energy as a 
soft one emits fewer photons. Furthermore, the spectrum 
is generally flat in the high energy range and it decays 
quickly at low energies. Therefore it is intrinsically more 
difficult to detect a harder burst. A study of the SMM 
data |93| suggests that there is a deficiency (by at least a 
factor of 5) of GRBs with hardness above 3MeV, relative 
to GRBs peaking at ~0.5MeV, but this data is consistent 
with a population of hardness that extends up to 2MeV. 

Overall the spectrum is non-thermal. This indicates 
that the source must be optically thin. The spectrum de- 
viates from a black body in both the low and the high 
energy ends: The X-ray paucity constraint rules out opti- 
cally thick models in which the 7-rays could be effectively 
degraded to X-rays |94) . The high energy tails lead to an- 
other strong constraint on physical GRB models. These 
high energy photons escape freely from the source without 
producing electron positron pairs! As we show later, this 
provides the first and most important clue on the nature 
of GRBs. 

The low energy part of the spectrum behaves in many 
cases like a power law: F v oc v a with — | < a < |, []l9], p5f . 
This is consistent with the low energy tail of synchrotron 
emission from relativistic electrons - a distribution of elec- 
trons in which all the population, not just the upper tail, is 
relativistic. This is a direct indication for the existence of 
relativistic shocks in GRBs. More than 90% of the bright 
bursts studied by Schaefer et al. [^7| satisfy this limit. 
However, there may be bursts whose low energy tail is 
steeper |96|. Such a spectrum cannot be produced by a 
simple synchrotron emission model and it is not clear how 
is it produced. 



2.4. Spectral Evolution 

Observations by earlier detectors as well as by BATSE 
have shown that the spectrum varies during the bursts. 
Different trends were found. Golenetskii et al. |nj ex- 
amined two channel data from five bursts observed by the 
KONUS experiment on Venera 13 and 14 and found a 
correlation between the effective temperature and the lu- 



minosity, implying that the spectral hardness is related to 
the luminosity. Similar results were obtained Mitrofanov 
et. al. Ipsj , Norris et al. investigated ten bursts 

seen by instruments on the SMM (Solar Maximum Mis- 
sion) satellite. They found that individual intensity pulses 
evolve from hard-to-soft with the hardness peaking earlier 
than the inte nsity . This was su ppor ted by more recent 
BATSE data po| . Ford et al. |L0| analyzed 37 bright 
BATSE bursts and found that the spectral evolution is a 
mixture of those found by Golenetskii et al. |)7j] and by 
Norris et al. |99| |: The peak energy either rises with or 
slightly proceeds major intensity increases and softens for 
the remainder of the pulse. For bursts with multiple peak 
emission, later spikes tend to be softer than earlier ones. 

A related but not similar trend is shown by the observa- 
tions that the burs ts are narrower at higher energ ies with 
T(v) oc v~ 0A 1 102]. As we show in section |8.7.3| th is be- 
havior is consistent with synchrotron emission | |103| . 

2.5. Spectral Lines 

Both absorption and emission features have been re- 
ported by various experiments prior to BATSE. Absorp- 
tion lines in the 20-40keV range have been observed by sev- 
eral experiments - but never simultaneously. GINGA has 
discovered several cases of lines with harmonic structure 



1 104, 105]. These lines were interpreted as cyclotron lines 
(reflecting a magnetic field of sa 10 12 Gauss) and provid- 
ing one of the strongest arguments in favor of the galactic 
neutron star model. Emission feat ures near 400keV have 
been claimed in other bursts [106 1 . These have been in- 
terpreted as red-shifted 511keV annihilation lines with a 
corresponding red-shift of « 20% due to the gravitational 
field on the surface of the Neutron star. These provided 
additional evidence for the galactic neutron star model. 

So far BATSE has not found any of the spectral features 
(abs o rptio n or emission lines) reported by earlier satellites 
1 107 , 108]. This can be interpreted as a problem with 
previous observations (or with the difficult analysis of the 
observed spectra) or as an unlucky coincidence. Given the 
rate of observed lines in previous experiments it is possible 
(at t he w 5% level) that the two sets of data are consistent 



Recently Meszaros & Rees [11C] suggested that within 
the relativitic fireball model the observed spectral lines 
could have be blue shifted iron X-ray line. 

2.6. Angular Positions 

BATSE is capable of estimating on its own the direc- 
tion to a burst. It is composed of eight detectors that are 
pointed towards different directions in the sky. The rela- 
tive intensity of the counts in the various detectors allows 
us to measure the direction to the burst. The positional er- 
ror of a given burst is the square root of the sum of squares 
of a systematic error and a statistical error. The statistical 
error depends on the strength of the burst. It is as large 
as 20° for a weak burst, and it is negligible for a strong 
one. The estimated systematic error (using a comparison 
of BATSE positions with IPN (Inter Planetary Network) 
localization) is k , 1.6 ° [111]. A different analysis of this 
comparison [112, 1 13 1 suggests that this might be slightly 
higher, around 3°. 
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The location of a burst is determined much better us- 
ing the difference in arrival time of the burst to several 
detectors on different satellites. Detection by two satel- 
lites limits the position to a circle on the sky. Detection 
by three determines the position and detection by four or 
more over-determines it. Even in this case the positional 
error depends on the strength of the bursts. The stronger 
the burst, the easier it is to identify a unique moment of 
time in the incoming signals. Clearly, the accuracy of the 
positional determination is better the longer the distance 
between the satellites. The best positions that have been 
obtained in this way are with the IPN 3 of detectors. For 12 
even ts the positional error boxes are of a few arc-minutes 

EH- 

BcppoSAX Wide Field Camera (WFC) that covers 
about 5% of the sky located a few bursts within 3' (3er). 
BcppoSAX's Narrow Field Instrument (NF1) obtained the 
bursts' positions to within 50". X-ray observations by 
ASCA and ROSAT have yielded error boxes of 30" and 
10" respectively. Optical identification has led, as usual, 
to a localization within 1" . Finally VLBI radio observation 
of GRB970508 has yielded a position within 200/iarcsec. 
The position of at least one burst is well known. 

2.7. Angular Distribution 

One of the most remarkable findings of BATSE was 
the observation that the angular distribution of GRBs' 
positions on the sky is perfectly isotropic. Early stud- 
ies had shown an isotropic GRB distribution [ 1 1 5 1 which 
have ev en l ed to the suggestion that GRBs are cosmo- 
logical [ 1 16 1 . In spite of this it was generally believed, 
prior to the launch of BATSE, that GRBs are associated 
with galactic disk neutron star. It has been expected 
that more sensitive detectors would discover an anisotropic 
distribution that would reflect the planar structure of 
the disk of the galaxy. BATSE's distribution is, within 
the statistical errors, in complete agreement with perfect 
isotropy. For the first 1005 BATSE bursts the observed 
dipole and quadrupole (corrected to BATSE sky expo- 
sure) relative to the galaxy are: (cos0) = 0.017 ± 0.018 
and (sin 2 & - 1/3) = -0.003 ± 0.009. These values are, 
respectively, 0.9a and 0.3cr from complete isotropy |||. 

2.8. Quiescent Counterparts and the historical "No 
Host" Problem 

One of the main obstacles in resolving the GRB mys- 
tery was the lack of identified counterparts in other wave- 
lengths. This has motivated numerous attempts to dis- 



cover GRB counterparts (for a review see [117, 118] ]). This 
is a difficult task - it was not known what to expect and 
where and when to look for it. 

The search for counterparts is traditionally divided to 
efforts to find a flaring (burst), a fading or a quiescent 
counterpart. Fading counterparts - afterglow - have been 
recently discovered by BeppoSAX and as expected this 
discovery has revolutionized GRB studies. This allowed 
also the discovery of host galaxies in seve ral cases, which 
will be discussed in the following section 2.9. Soft X-ray 



a separate component. Flaring has not been discovered in 
other wavelengths yet. Quiescent counter parts were not 
discovered either. 

Most cosmological models suggest that GRBs are in a 
host galaxy. If so then deep searches within the small 
error boxes of some GRBs localized by the IPN system 
should reveal the host galaxy, until the discovery of GRB 
afterglow these searches have yielded only upper limits on 
the magnitudes of possible hosts. This has lead to what 
is called the "No Host" Problem. Schaefer etal. con- 
ducted searches in the near and far infrared lll9| using 



flaring (simultaneous with the GRB) was discovered in sev- 
eral bursts but it is an ambiguous question whether this 
should be considered as a part of the GRB itself or is it 



IRAS, in radio us ing t he VLA |12C] and in archival op- 
tical photographs [121] and have found only upper limits 
and no clear counterpart candidates. Similar results from 
multiple wavelength observations have been ob tain ed by 
Hurley et al. |^]. Vrba, Hartmann & Jennings [122] have 
monitored the error boxes of seven bursts for five year. 
They did not find any unusual objects. As for the "no 
host " problem this authors, as well as L uginbuhl et al. 
|123} | and Larson, McLean & Becklin [ 124 1 concluded, us- 
ing the standard galaxy luminosity function, that there 
are enough dim galaxies in the corresponding GRB error 
boxes which could be the hosts of cosmological burst and 
therefore, there is no "no host" pro blem . 

More recently Larson & McLean [ 125 1 monitored in the 
infrared nine of the smallest error boxes of burst localized 
by the IPN with a typical error boxes of eight arc-min 2 . 
They found in all error boxes at least one bright galaxy 
with K < 15.5. However, the error boxes are too large to 
discern between the host galax y and unrelated background 
galaxies. Schaefer et al. ]126| , searched the error boxes of 
five GRBs using the HST. Four of these are smaller boxes 
with a size of ~ 1 arc-min 2 . They searched but did not find 
any unusual objects with UV excess, variability, parallax 
or proper motion. They have found only faint galaxies. 
For the four small error boxes the luminosity upper lim- 
its of the host galaxy are 10-100 times smaller tha n the 
luminosity of an L„ galaxy. Band & Hartmann |127|1 con- 
cluded that the error boxes of Larsen & McLean | |125| ] are 
too large to discriminate between the presence or the ab- 
sence of host galaxies. However, they fi nd th at the absence 
of host galaxies in the Schaefer et al. [ 126 1 data is signif- 
icant, at the 2 • 10 -6 level. Suggesting that there are no 
bright hosts. 

This situation has drastically changed and the "no 
host" problem has disappeared with afterglow observa- 
tions. These observations have allowed for an accurate 
position determination and to identification of host galax- 
ies for several GRBs. Most of these host galaxies are dim 
with magnitude 24.4 < R < 25.8. This support the con- 
clusions of the earlier studies that GRBs are not associated 
with bright galaxies and definitely not with cores of such 
galaxies (ruling out for example AGN type models). These 
observations are consistent with GRBs rate being either a 
cons tant or being proportional to the star formation rate 
|128} |. According to this analysis it is not surprising that 
most hosts are detected at R ~ 25. However, though these 
two models are consistent with the current data both pre- 
dict the existence of host galaxies brighter than 24 mag, 
which were not observed so far. One could say now that 
the "no host" problem has been replaced by the "no bright 
host" problem. But this may not be a promlem but rather 
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an indication on the nature of the sources. 

The three GRBs with measured cosmological redshifts 
lie in host galaxies with a strong evidence for star forma- 
tion. These galaxies display prominent emission lines from 
line associated with star-formation. In all three cases the 
strength of those lines is hi gh for ga l axies of c omparable 



magnitude and redshift |lq, 129, 130, 131, 128 1. The host 



of GRB980703, for example, show a star forming rate of 
~ l OAffpyr -1 or higher with a lower limit of TM^yx 
1 129 ] . For most GRBs with afterglow the host galaxy was 



detected but no emission or absorption lines were found 
and no redshift was measured. This result is consistent 
with the hypothesis that all GRBs are associated with 
star-forming galaxies. For those hosts that are at red- 
shift 1.3 < z < 2.5 the corresponding emission lines are 
not observed as for this redshift range no s tron g lines are 
found in the optical spectroscopic window [ 1 3 1 1 . 

The simplest conclusion of the above observations is that 
all GRBs are associated with star forming regions. Still 
one has to keep in mind that those GRBs on which this 
conclusion was based had a strong optical afterglow, which 
not all GRBs show. It is possible that the conditions asso- 
ciated with star forming regions (such as high interstellar 
matter density - or the existance of molecular clouds) are 
essential for the appearance of strong optical afterglow and 
not for the appearance of the GRB itself. 

2.9. Afterglow 

GRB observations were revolutionized on Feb ruar y 28, 
1997 by the Italian-Dutch satellite BeppoSAX [ |l32| that 
discovered an X-ray counterpart to GRB970228 0. 
GRB970228 was a double peaked GRB. The first peak 
which lasted ~ 15sec was hard. It was followed, 40 seconds 
later, by a much softer second peak, which lasted some 
~ 40sec. The burst was detected by the GRBM (Gamma- 
Ray Burst Monitor) as well as by the WFC (Wide Field 
Camera). The WFC, which has a 40° x 40° field of 
view detected soft X-rays simultaneously with both peaks. 
Eight hours latter the NFI (Narow Field Instrument) was 
pointed towards the burst's directions and detected a con- 
tinuous X-ray emission. The accurate position determined 
by BeppoSAX enabled the identification of an optical af- 
terglow jll]] - a 20 magnitude poi nt source adjacent to a red 
nebulae. HST observations [133] revealed that the nebula 
adjacent to the source is roughly circular with a diameter 
of 0" .8. The diameter of the nebula is comparable to the 
one of galaxies of similar magnitude found in the Hubble 
Deep Field, especially if one takes into account a possible 
visual extinction in the direction of GRB970228 of at least 



135| 



one magnitude [ 134 1 . 

Follow ing X-ray detect ions by BeppoSAX |lC_ 
ROSAT pi and ASCA pi revealed a decaying X^ray 
flux oc t 1M±0 - 11 (see Fig. 0). The decaying flux can be 
extrapolated as a power law directly to the X-ray flux of 
the second peak (even though this extrapolation requires 
some care in determining when is t = 0) 



The optical emission also depicts a decaying flux [ 138 1 
(see fig. U The source could not be observed from late 
March 97 until early September 1997. When it was ob- 



served again, on Sept. 4th by HST [139, 131 it was found 
that the optical nebulosity does not decay and the point 
source shows no proper motion, refuting earlier sugges- 



tions. The visual magnitude of the nebula on Sept. 4th 
was 25.7 ± 0.25 compared with V = 25.6 ± 0.25 on March 
26th and April 7th. The visual magnitude of the point 
source on Sept. 4th was (V — 28.0 ± 0.25), whi ch i s con- 
sistent with a decay of the flux as t- 1A4±0 - Q5 Jl3^] . In 
spite of extensive efforts no radio emission was detected 
and one can set an uppe r limit of ~ 10/iJy to the radio 
emission at 8.6 Ghz [140 1 . 

GRB970508 was detected by both BATSE in 7 -rays 
[ Till ] an d BeppoSAX in X-rays (lj| on 8 May 1997. The 
7-ray burst lasted for ~ 15sec, with a 7-ray fluence of 



3 x 10 6 ergs/cm -. Var i able emission in X- rays 



tical JTJ, 



143 



144, 145 



146 



op- 



147| and radio (TJ, [l48| fol- 



lowed the 7-rays. The spectrum of the optical transient 
taken by Keck revealed a set of absorption lines associated 
with Fe II and Mg II and O II emission line with a red- 
shift z = 0.835 1 14 1 . A second absorption line system with 
z = 0.767 is also seen. These lines reveal the existence of 



an underlying, dim gala xy h ost. HST images [149, 150 
and Keck observations [ 130 1 show that this host is very 
faint (R = 25.72 ± 0.2 mag), compact (< 1 arcsec) dwarf 
galaxy at z = 0.835 and nearly coincident on the sky with 
the transient. 

The optical light curve peaks at around 2 days af- 
ter the burst. Assuming isotropic emission (and using 
z = 0.835 and H=100km/sec/Mpc) this peak flux cor- 
responds to a luminosity of a few x 10 45 ergs/sec. The flux 
decline shows a con t inuo us power law decay cx t~ 127±om 
[151 



152, 153, 154 



1301 - After about 100 days the light 
curve begun to flatten as th e tra n sient fade d and become 
weaker than the host 1)1 5 3|, 155, 156, 157]. Integration 



of this light curve results in an overall emission of a few 
xlO 50 ergs in the optical band. Radio emission was ob- 
served, first, one week after the burst |lj| (see Fig. ||). 
This emission showed int ensiv e oscillations which were in- 
terpreted as scintillation ]158| . The subsequent disappear- 
ance of these oscillations after about three weeks enabled 
Frail et al. [|l3| to estimate the size of the fireball at this 
stage to be ~ 10 17 cm. This was supported by the indi- 
cation that the radio emission was initially optically thick 
|lj| , which yields a similar estimate to the size 25 1 . 

GRB970828 was a strong GRB that was detected by 
BATSE on August 28, 1997. Shortly afterwards RXTE 
|159| , 160] focused on the approximate BATSE position 
and discovered X-ray emission. This X-ray emission de- 
termined the position of the burst to within an elliptical 
error box with 5' x 2'. However, in spite of enormous ef- 
fort no variable optical counterpart brighter than R=23.8 
that h as ch anged by more than 0.2 magnitude was de- 
tected [161]. There was also no indication of any radio 
emission. Similarly X-ray afterglow was detected from sev- 
eral other GRBs (GRB970615, GRB970402, GRB970815, 
GRB980519) with no optical or radio emission. 

Seventeen GRBs have been detected with arcminute po- 
sitions by July 22, 1998: fourteen by the WFC of Bep- 
poSAX and three by the All-Sky Monitor (ASM) on board 
the Rossi X-ray Timing Explorer (RXTE). Of these sev- 
enteen burst, thirteen were followed up within a day in 
X-rays and all those resulted in good candidates for X- 
ray afterglows. We will not discuss all those here (see 
table 4 for a short summary of some of the properties). 
Worth mentioning are however, GRB971214, GRB980425 
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and GRB980703. 

GRB971214 was a rather strong burst. It was detected 
on December 14.9 UT 1997 [162]. Its optical counterpart 
was observed wit h a magnitude 21.2±0.3 on the I band by 
Halpern et al., [ 163 1 on Dec. 15.47 UT twelve hours after 
the burst. It was observed one day later on Dec. 16.47 
with I magnitude 22.6. Kulkarni et al. Jlj| obtained a 
spectrum of the host galaxy for GRB971214 and found a 
redshift of z=3.418! With a total fluence of 1.09 x 10~ 5 ergs 
cm~ 2 [164 1 this large redshift implies, for isotropic emis- 



sion, Q = 1 and Hq = 65km/sec/Mpc, an energy release of 
~ 10 53 ergs in 7-rays alone 4 . The familiar value of 3 x 10 53 
lla] is obtained for fl = 0.3 and Hq = 0.55km/sec/Mpc. 

GRB980425 was a moderately weak burst with a peak 
flux of 3 ±0.3 x 10 _7 ergs cm -2 sec -1 . It was a single peak 
burst with a rise time of 5 seconds and a decay time of 
about 25 seconds. The burst was detected by BeppoSAX 
(as well as by BATSE) whose WFC obtained a position 
with an error box of 8'. Inspection of an image of this 
error box taken by the New Technology Telescope (NTT) 
revealed a type Ic supernova SN1998bw tha t too k place 
more or less at the same time as the GRB [151|. Since 



the probability for a chance association of the SN and the 
GRB is only 1.1 x 10~ 4 it is likely that this association is 
real. The host galaxy of this supernova (ESO 184-G82) has 
a redshift of z = 0.0085 ± 0.0002 putting it at a distance 
of 38 ± IMpc for H=67km/sec Mpc. The corresponding 
7-ray energy is 5 x 10 47 ergs. With such a low luminosity 
it is inevitable that if the association of this burst with 
the supernova is real it must correspond to a new and rare 
subgroup of GRBs. 

GRB980703 was a very strong burst with an obs erved 
gamma-ray fluence of (4.59±0.42) x 10~ 5 ergs cm -2 [166|. 
Keck observations revealed that the host galaxy has a red- 
shift of z = 0.966. The corresponding energy release (for 
isotropic emiss ion, — 0.2 and Hq = 65km sec /Mpc) 



is - 10 53 ergs g29] 



2.10. Repetition? 

Quashnock and Lamb [167 ] suggested that there is ev- 
idence from the data in the BATSE IB catalog for repe- 
tition of bursts from the same source. If true, this would 
severely constrain most GRB models. In particular, it 
would rule out any other model based on a 'once in life- 
time' catastrophic event. This claim has been refuted by 



several authors 1 16S , 169[ and most notably by the analysis 
of the 2B data |fg] and 3B data Jm|. 

A unique group of four bursts - "the gang of four" - 
emerged from the same position on the sky within two days 
||]. One of those bursts (the third burst GRB961029a) 
was extremely strong, one of the strongest observed by 
BATSE so far. Consequently it was observed by the 
IPN network as well, and its position is known accu- 
rately. The other three (GRB961027a,GRB961027b and 
GRB961029d) were detected only by BATSE. The precise 
position of one burst is within the la circles of the three 
other bursts. However, two of the bursts are almost 3<r 
away from each other. Is this a clear cut case of repetition? 
It is difficult to assign a unique statistical significance to 
this question as the significance depends critically on the 
a priori hypothesis that one tests. Furthermore, the time 

4 This value depends also on the spectral shape of the burst. 



difference between the first and the last bursts is less than 
two days. This is only one order of magnitude longer than 
the longest burst observed beforehand. It might still be 
possible that all those bursts came from the same source 
and that they should be considered as one long burst. 

2.11. Correlations with Abell Clusters, Quasars and 
Supernovae 

Various attempts to search for a correlation between 
GRBs and other astronomical objects led to null result. 
For example, Blumenthal & Hartmann [172] found no 
angular correlation between GRBs and nearby galaxies. 
They concluded that if GRBs arc cosmological then they 
must be located at distances larger than lOOMpc. Other- 
wise, they would have shown a positive correlation with 
the galaxy distribution. 

The only exception is the correlation (at 95% confidence 
leve l) be tween GRBs at the 3B catalog and Abell clusters 



|77|, 172 ] . This correl ation has been recently confirmed by 
Kompaneetz & Stern [ 174 1 . The correlation is strongest for 



a subgroup of strong GRBs whose position is accurately 
known. Comparison of the rich clusters auto-correlation 
with the cross-correlation found suggests that ~ 26 ± 15% 
of the accurate position GRBs sub-sample members are 
located within 600 h _1 Mpc. Recently Schartel et al. [176| 
found that a group of 134 GRBs with position error radius 
smaller than 1.8° are correlated with radio quiet quasars. 
The probability of of such correlation by chance coinci- 
dence is less than 0.3%. 

It should be stressed that this correlation does not im- 
ply that there is a direct association between GRBs and 
Abell clusters, such as would have been if GRBs would 
have emerged from Abell clusters. All that it means is that 
GRBs are distributed in space like the large scale structure 
of the universe. Since Abell clusters are the best tracers 
of this structure they are correlated with GRBs. There- 
fore the lack of excess Abell Clusters in IPN error b oxes 
(which are much smaller than BATSE's error boxes) [|l75| 
does not rule out this correlation. 



2.12. V/V m ax) Count and Peak Flux Distributions 

The limiting fluence observed by BATSE is w 
10 _7 ergs/cm 2 . The actual fluence of the strongest bursts 
is larger by two or three orders of magnitude. A plot of 
the number of bursts vs. the peak flux depicts clearly a 
paucity of weak bursts. This is manifested by the low value 
of (V/V max ), a statistic d esig ned to measure the distribu- 
tion of sources in space [177|. A sample of the first 601 
bursts has (V/V ma x) = -328 ± 0.012, which is 14a away 
from the homogeneous flat space value of 0.5 [178|. Cor- 
respondingly, the peak count distribution is incompatible 
with a homogeneous population of sources in Euclidean 
space. It is compatible, however, with a cosmological dis- 
tribution (see Fig. ITTi). The distribution of short bursts 
has a larger (V/V m a^) and it is com pat ible with a homo- 
geneous Eucleadian distribution JF^ [78|, f79[ . 



3. THE DISTANCE SCALE 
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Table 1: Observational data of several GRBs for which afterglow was detected. The two columns O and R indicate 
whether emission was detected in the optical and radio, respectively. The total energy of the burst is estimated through 
the observed fluence and redshift, assuming spherical emission and a flat Q = 1, A = universe with Hq = 65Km/sec/Mpc. 



3.1. Redshift Measurements. 

The measurements of redshifts of several GRB optical 
counterparts provide the best and the only direct distance 
estimates for GRBs. Unfortunately these measurements 
are available only for a few bursts. 

3.2. The Angular Distribution 

Even before these redshift measurements there was a 
strong evidence that GRBs originate from cosmological 
distances. The observed angular distribution is incompat- 
ible with a galactic disk distribution unless the sources are 
at distances less than lOOpc. However, in this case we 
would expect that (V/V max) = 0.5 corresponding to a ho- 
mogeneous distribution [L77| while the observations yield 
(V/V max ) = 0.33. 

A homogeneous angular distribution could be produced 
if the GRB originate from the distant parts of the galactic 
halo. Since the solar system is located at d = 8.5kpc from 
the galactic center such a population will necessarily have 
a galactic d ipolc of order d/R, where R is a typical distance 
to a GRB f 
constrain t 



179 1 . The lack of an observed dipole strongly 
lis model. Such a distribution of sources is 
incompatible with the distribution of dark matter in the 
halo. The typical distance to the GRBs must of the order 
of lOOkpc to comply with this constraint. For example, if 
one considers an effective distribution that is confined to 
a shell of a fixed radius then such a shell would have to 
be at a distance of lOOkpc i n ord er to be compatible with 
current limits on the dipole [18C ] . 



3.3. Interpretation of the Peak Flux Distribution 

The counts distribution or the Peak flux distribution of 
the bursts observed by BATSE show a paucity of weak 
burst. A homogeneous count distribution, in an Eucle- 
adian space should behave like: N(C) cx C~ 3 / 2 , where 
N(C) is the number of bursts with more than C counts 
(or counts per second). The observed distribution is much 
flatter (see Fig. pi). This fact is reflected by the low 
(V/Vmax) value ofthe BATSE data: there are fewer dis- 
tant sources than expected. 



The observed distribution is compatible with a cosmo- 
logical distribution of sources. A homogeneous cosmolog- 
ical distribution displays the observed trend - a paucity 
of weak bursts relative to the number expected in a Eu- 
cleadian distribution. In a cosmological population four 
factors combine to make distant bursts weaker and by this 
to reduce the rate of weak bursts: (i) K correction - the 
observed photons are red-shifted. As the photon number 
decreases with energy this reduces the count rate of dis- 
tant bursts for a detector at a fixed energy range, (ii) The 
cosmological time dilation causes a decrease (by a factor 
1 + z) in the rate of arrival of photons. For a detector, like 
BATSE, that measures the count rate within a given time 
window this reduces the detectability of distant bursts. 

(iii) The rate of distant bursts also decreases by a factor 
1 + z and there are fewer distant bursts per unit of time 
(even if the rate at the comoving frames does not change). 

(iv) Finally, the distant volume element in a cosmological 
model is different than the corresponding volume element 
in a Eucleadian space. As could be expected, all these 
effects are significant only if the typical red-shift to the 
sources is of order unity or larger. 

The statistics (V/V max ) > is a weighted average of the 
distribution N(> /). Already in 1992 Piran [M com- 
pared the theoretical estimate of this statistics to the ob- 
served one and concluded that the typical redshift of the 
bursts obser ved by BATSE is z max ~ 1. Later Fenimore 
et al. ]200| compared the sensitivity of PVO (that ob- 
serves N(> /) cx /~ 3/2 ) with the sensitivity of BATSE and 
concluded that z max (BATSE) ~ 1 (the maximal z from 
which bursts are detected by BATSE) . This corresponds to 
a peak luminosity of ~ 10 50 ergs/sec. Other calculations 
based on different stat istical methods were per formed by 
Horack & Emslie[186|, Loredo & Wa sserm an [ 181 , |182| , 
Rutledge et al. [ 185|| Cohen fc Piran [183] and Meszaros 
and collaborators^ 187, 188, 189 , |l 90f| a nd others. In par- 
ticular Loredo & Wassermann] 181 , 182] give an extensive 
discussion of the statistical methodology involved. 

Consider a homogeneous cosmological distribution of 
sources with a peak luminosity L, that may vary from 
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one source to another. It should be noted that only the 
luminosity per unit solid angle is accessible by these ar- 
guments. If there is significant beaming, as inferred f25| . 
the distribution of total luminosity may be quite different. 
The sources are emitting bursts with a count spectrum: 
N(v)dv — (L/hv)N(v)dv, where hv, is the average en- 
ergy. The observed peak (energy) flux in a fixed energy 
range, [E m i n , E max \ from a source at a red-shift z is: 



f(L,z) 



(1 + aQ L 
Andf(z) hv 



N[u(l+z)]hu(l+z)hdp (3) 



where di(z) is the luminosity distance 1 184]. 

To estimate the number of bursts with a peak flux larger 
than /, N(> /), we need the luminosity function, ip(L, z): 
the number of bursts per unit proper (comoving) volume 
per unit proper time with a given luminosity at a given 
red-shift. Using this function we can write: 



N(> f ) = 4tt 



■■(f,L) 



df 



dr p {z) 



(1 + z) 3 dz 



dzdL 



(4) 

where the red-shift, z(f,L), is obtained by inverting Eq. 
H and r p (z) is the proper distance to a red-shift z. For a 
given theoretical model and a given luminosity function we 
can calculate the theoretical distribution N(f) and com- 
pare it with the observed one. 

A common simple model assumes that iJj(L, z) = 
4>(L)p(z) - the luminosity does not change with time, but 
the rate of events per unit volume per unit proper time 
may change. In this case we have: 



N(> /) = 4tt / </>(L) 



df 



dr p {z) 



. (5) 

The emitted spectrum, N(v), can be estimated from the 
observed data. The simplest shape is a single power law 
(Eq. |). with a = 1.5 or a — 1.8|85|. More elaborate 
studies have used the Band etal. p4j spectrum or even a 
distribution of such spectra ]185[ | . 

The cosmic evolution function p(z) and the luminosity 
function (j)(L) are unknown. To proceed one has to choose 
a functional shape for these functions and characterize it 
by a few parameters. Then using maximum likelihood, or 
some other technique, estimate these parameters from the 
data. 

A simple characterization of p{z) is: 



p(z) = p {l + z)-f 3 . 



(6) 



Similarly the simplest characterization of the luminosity is 
as standard candles: 



4>{L)=6(L-Lo), 



(7) 



with a single parameter, Lq, or equivalcntly z max , the 
maximal z from which the source is detected (obtained 
by inverting Eq. | for / = f mm and L = L ). 

There are two unknown cosmological parameters: the 
closure parameter, il, and the cosmological constant A. 
With the luminosity function given by Eqs. (5) and (6) we 
have three unknown parameters that determine the bursts' 



distribution: Lq, po, (3. We calculate the likelihood func- 
tion over this five dimensional parameter space and find 
the range of acceptable models (those whose likelihood 
function is not less than 1% of the maximal likelihood). 
We then proceed to perform a KS (Kolmogorov-Smirnov) 
test to check whether the model with the maximal likeli- 
hood is an acceptable fit to the data. 

The likelihood function is practically independent of Q 
in the range: 0.1 < 57 < 1. It is also insensitive to the cos- 
mological constant A (in the range < A < 0.9, in units 
of the critical density) . This simplifies the analysis as we 
are left only with the intrinsic parameters of the bursts' 
luminosity function. 

There is an interplay between evolution (change in the 
bursts' rate) and luminosity. Fig. ^ depicts the likelihood 
function in the (z max , (3) plane for sources with a varying 
intrinsic rate. The banana shaped contour lines show that 
a population whose rate is independent of z ((3 = 0) is 
equivalent to a population with an increasing number of 
bursts with cosmological time (/3 > 0) with a lower Lq 
(lower z max ). This tendency saturates at high intrinsic 
evolution (large /3), for which the limiting z max does not 
go below .5 and at very high Loj for which the limit- 
ing [3 does not decrease below -1.5. This interplay makes 
it difficult to constraint the red shift distribution of GRB 
using the peak flux distribution alone. For completeness 
we quote here "typical" results based on standard cand les, 
no evolution and an Einstein-DeSitter cosmology [183]. 

Recall that (V/V max } of the short bursts distribution is 
rather close to the homogeneous Eucleadian value of 0.5. 
This means that when analyzing the peak flux distribution 
one s hould analyze separately the long and the short bursts 
|183[ ]. For long bursts (bursts with tgo > 2 sec) the likeli- 
hood function peaks at z max — 2.1 (see Fig. Q [ 183|. The 
allowed range at a 1% confidence level is: 1.4 < z max < 3.1 
(ztoi 2 ' = 1-5(^4) f° r a = 2). The maximal red-shift, 

Zmax — 

2.l[+^, corresponds, with an estimated BATSE 



detection efficiency of ~ 0.3, to 2.3 



(-0.7) 



10 6 events per 



galax y per year (for a galaxy density of 10~ 2 h 3 Mpc - ; 
|192[| ). The rate per galaxy is independent of Hq and is 
only weakly dependent on Q. For fi = 1 and A = the 
typical energy of a burst with an observed fluence, F, is 
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(+ii) 



10 



no. 



(_4) ' iu (F/10 7 ergs/cm )ergs. The distance to the 
sources decreases and correspondingly the rate increases 
and the energy decreases if the spectral index is 2 and not 
1.5. These numbers vary slightly if the bursts have a wide 
luminosity function. 

Short bursts are detected only up to a much nearer dis- 
tances: z max (short) = 0.4 +11 , again assuming standard 
candles and no source evolution. There is no significant 
lower limit on z max for short bursts and their distribution 
is compatible with a homogeneous non-cosmological one. 
The estimate of z max (short) corresponds to a comparable 
rate of 6.3(_ 5 6 ) • 10 -6 events per year per galaxy and a 
typical energy of 3^ +39 - ) • 10 49 J 7 '_7 ergs (there are no lower 
limits on the energy or and no upper limit on the rate 
since there is no lower limit on z max (short)) . The fact 
that short bursts are detected only at nearer distances is 
also reflected by the higher (V/V max ) of the population of 
these bursts 1791. 
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Relativ ely w ide luminosity distributions are allowed by 
the data [ 183 1 . For example, the KS test gives a proba- 
bility of 80% for a double peaked luminosity distribution 
with luminosity ratio of 14. These results demonstrate 
that the BATSE data alone allow a variability of one or- 
der of magnitude in the luminosity. 

The above considerations should be modified if the 
rate of GRBs trace the SFR - the star formation rate 
1 193) , [l9j, EH- The SFR has been determined recently 
by two independent studies |l9fj, |l97], |l9§]. The SFR 
peaks at z ~ 1.25. This is a strongly evolving non mono- 
tonic distribution, which is drastically different from the 
power laws considered so far. Sahu et al. [ |194[ find 
that p(z) cx SFR(z) yields N(> /) distribution that is 
compatible with the observed one (for qo = 0.2, Ho = 
50km/sec _1 Mpc _1 ) for a narrow lumi nosit y distribution 
with L 1 = 10 51 ergs/sec. Wijers et al. [195] find that the 
implied peak luminosity is higher L 1 — 8.3 • 10 51 ergs/sec 
and it corresponds to a situation in which the dimmest 
bursts observed by BATSE originate from z ~ 6! 

The direct red-shift measure of GRB970508 (TJ] agrees 
well with e stim a tes made previously using peak-flux count 
statistics ( |200[ |J|, (TJ|]). The red-shift of GRB971214, 
z = 3.418, and of GRB980703, z = 0.966, and the implied 
luminosities disagree with these estimates. A future de- 
tection of additional red-shifts for other bursts will enable 
us to estimate directly the luminosity function of GRBs. 
It will also enable us to determine the evolution of GRBs. 



Krumholz et al. (199| and Hogg & Fruchter |128| find that 
with a wide luminosity function both models of a constant 
GRB rate and a GRB rate following the star formation rate 
are consistent with the peak flux distribution and with the 
observed redshift of the three GRBs. 

3.4. Time Dilation 

Norris et al. J202| . |203| examined 131 long bursts (with 
a duration longer than 1.5s) and found that the dimmest 
bursts are longer by a factor of 2.3 compared to the 
bright ones. With our canonical value of z max = 2.1 the 
bright bursts originate at Zbright ~ 0.2. The corresponding 
expected ratio due to cosmological time dilation, 2.6, is in 
agre ement with this measurement. Fenimore and Bloom 
1 204] find, on the other hand, that when the fact that 
the burst's duration decreases as a function of energy as 
At w E~ a - 5 is included in the analysis, this time dila- 
tion corresponds to z max > 6. This would require a strong 
negative intrinsic evolution: /? ss — 1.5±0.3. Alternatively, 
this might agree with the model in which the GRB rate 
follows the SFR £jj|, 1951 which gives z max ss 6. Cohen & 
Piran [191| suggested a way to perform the time dilation, 
spectral and red shift analysis simultaneously. Unfortu- 
nately current data are insufficient for this purpose. 



4. THE COMPACTNESS PROBLEM AND RELATIVISTIC 
MOTION. 

The key to understanding GRBs lies, I believe, in un- 
derstanding how GRBs bypass the compactness problem. 
This proble m w as realized very early on in one for m by 
Ruderman [205] and in another way by Schmidt [ 206 1 . 
Both used it to argue that GRBs cannot originate from 
cosmological distances. Now, we understand that GRBs 



are cosmological and special relativistic effects enable us 
to overcome this constraint. 

The simplest way to see the compactness problem is to 
estimate the average opacity of the high energy gamma- 
ray to pair production. Consider a typical burst with an 
observed fluence, F. For a source emitting isotropically 
at a distance D this fluence corresponds to a total energy 
release of: 



E = 4ttD 2 F = 10 50 ergs( 
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VSOOOMpc/ VlO~ 7 ergs/cm : 

Cosmological effects change this equality by numerical 
factors of order unity that are not important for our dis- 
cussion. The rapid temporal variability on a time scale 
ST rs 10 msec implies that the sources are compact with 
a size, Ri < cST » 3000 km. The observed spectrum (see 



section contains a large fraction of high energy 7-ray 
photons. These photons (with energy E\) could interact 
with lower energy photons (with energy E2) and produce 
electron-positron pairs via 77 — > e + e~ if \j E\Ei > m e c 2 
(up to an angular factor). Denote by f p the fraction of 
photon pairs that satisfy this condition . Th e average op- 
tical depth for this process is |o|, |o|, |o|: 
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where <tt is the Thompson cross-section. This optical 
depth is very large. Even if there are no pairs to begin with 
they will form rapidly and then these pairs will Compton 
scatter lower energy photons, resulting in a huge optical 
depth for all photons. However, the observed non-thermal 
spectrum indicates with certainty that the sources must 
be optically thin! 

An alternative calculation is to consider the optical 
depth of the highest energy photons (say a GeV photon) 
to pair production with the lower energy photons. The 
observation of GeV photons shows that they are able to 
escape freely. In other words it means t hat this optical 
depth must be much smaller than unity J210| , 211]. This 
consideration leads to a slightly stronger but comparable 
limit on the opacity. 

The compactness problem stems from the assumption 
that the size of the sources emitting the observed radia- 
tion is determined by the observed variability time scale. 
There won't be a problem if the source emitted the en- 
ergy in another form and it was converted to the observed 
gamma-rays at a large distance, Rx, where the system 
is optically thin and r 77 (i?x) < 1. A trivial solution of 
this kind is based on a weakly interacting particle, which 
is converted in flight to electromagnetic radiation. The 
only problem with this solution is that there is no known 
particle that can play this role (see, however |213| ). 

4.1. Relativistic Motion 

Relativistic effects can fool us and, when ignored, lead 
to wrong conclusions. This happened thirty years ago 
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when rapid variability implied "impossible" temperatures 
in extra-galactic radio sources. This puzzle was resolved 
when it was suggested [214, 215] that these objects reveal 
ultra-relativistic expansion. This was confirmed later by 
VLBA measurements of superluminal jets with Lorentz 
factors of order two to ten. This also happened in the 
present case. Consider a source of radiation that is mov- 
ing towards an observer at rest with a rclativistic velocity 
characterized by a Lorentz factor, 7 = 1/ ^/l — v 2 /c 2 ^S> 1. 
Photons with an observed energy hu f, s have been blue 
shifted and their energy at the source was « hv b s /"f- Since 
the energy at the source is lower fewer photons have suffi- 
cient energy to produce pairs. Now the observed fraction 
f p , of photons that could produce pairs is not equal to the 
fraction of photons that could produce pairs at the source. 
The latter is smaller by a factor -f~ 2a (where a is the high 
energy spectral index) than the observed fraction. At the 
same time, relativistic effects allow the radius from which 
the radiation is emitted, R e < j 2 cST to be larger than the 
original estimate, R e < cST, by a factor of j 2 . We have 
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where the relativistic limit on R e was included in the sec- 
ond line. The compactness problem can be resolved if 
the source is moving relativistically towards us with a 
Lorentz factor 7 > io 13 /( 4 + 2q ) « 10 2 . A more detailed 
discussion [210, 211] gives comparable limits on 7. Such 
extreme-relativistic motion is larger than the relativistic 
motion observed in any other celestial source. Extragalac- 
tic super-luminal jets, for example, have Lorentz fact ors of 
~ 10, while the known galactic relativistic jets |216| have 
Lorentz factors of ~ 2 or less. 

The potential of relativistic motion to resolve the com- 
pactness problem was realized in the eigh ties by Goodman 
|217t , Paczyhski @ and Krolik and Pier fl21S| ]. There was, 
however, a difference be twee n the first two approaches and 
the last one. Goodman [217] and Paczyhski [B3[ considered 
relativistic motion in the dynamical context of fireballs, in 
which the relativistic motion is an integr al pa rt of the dy- 
namics of the system. Krolik and Pier [21S] considered, 
on the other hand, a kinematical solution, in which the 
source moves relativistically and this motion is not neces- 
sarily related to the mechanism that produces the burst. 

Is a purely kinematic scenario feasible? In this scenario 
the source moves relativistically as a whole. The radia- 
tion is beamed with an opening angle of 7 -1 . The total 
energy emitted in the source frame is smaller by a factor 
7~ 3 than the isotropic estimate given in Eq. The total 
energy required, however, is at least (Mc 2 + AttFD 2 /j 3 )j, 
where M is the rest mass of the source (the energy would 
be larger by an additional amount E t h"f if an internal 
energy, E t h, remains in the source after the burst has 
been emitted). For most scenarios that one can imag- 
ine Mc 2 7 3> (Att/^ 2 )FD 2 . The kinetic energy is much 
larger than the observed energy of the burst and the pro- 
cess is extremely (energetically) wasteful. Generally, the 



total energy required is so large that the model becomes 
infeasible. 

The kinetic energy could be comparable to the observed 
energy if it also powers the observed burst. This is the 
most energetically-economical situation. It is also the most 
conceptually-economical situation, since in this case the 7- 
ray emission and the relativistic motion of the source are 
related and are not two independent phenomena. This 
will be the case if GRBs result from the slowing down 
of ultra relativistic mat ter. This idea was suggested by 
Meszaros, and Rees |27], 219] in the c onte xt of the slowing 
down of fireball accelerated material [ 22C ] by the ISM and 
by Narayan, et al. p8| and independently by Rees and 
Meszaros |2^] and Paczyhski and Xu |50| in the context 
of self interaction and internal shocks within the fireball. 
It is remarkable that in both cases the introduction of en- 
ergy conversion was motivated by the need to resolve the 
"Baryonic Contamination" problem (which we discuss in 
the next section). If the fireball contains even a small 
amount of baryons all its energy will eventually be con- 
verted to kinetic energy of those baryons. A mechanism 
was needed to recover this energy back to radiation. How- 
ever, it is clear now that the idea is much more general 
and it is an essential part of any GRB model regardless of 
the nature of the relativistic energy flow and of the specific 
_ 2 way it slows down. 

Assuming that GRBs result from the slowing down of a 
relativistic bulk motion of massive particles, the rest mass 
of the ultra-relativistic particles is: 
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where e c is the conversion efficiency and 9 is the opening 
angle of the emitted radiation. We see that the allowed 
mass is very small. Even though a way was found to con- 
vert back the kinetic energy of the baryons to radiation 
(via relativistic shocks) there is still a "baryonic contami- 
nation" problem. Too much baryonic mass will slow down 
the flow and it won't be relativistic. 



4.2. Relativistic Beaming? 

Radiation from relativistically moving matter is beamed 
in the direction of the motion to within an angle 7 -1 . 
In spite of this the radiation produced by relativistically 
moving matter can spread over a much wider angle. This 
depends on the geometry of the emitting region. Let 6 m 
be the angular size of the relativistically moving matter 
that emits the burst. The beaming angle 9 will be 9m if 
9m > 7~ 1 and 7 _1 otherwise. Thus if 9m = 47r - that 
is if the emitting matter has been accelerated spherically 
outwards from a central source (as will be the case if the 
source is a spherical fireball) - the burst will be isotropic 
even though each observer will observe radiation coming 
only from a very small region (see Fig. 0). The radia- 
tion will be beamed into 7 _1 only if the matter has been 
accelerated along a very narrow beam. The opening angle 
can also have any intermediate value if it emerges from 
a beam with an opening angle 9 > 7" 1 , as will be the 
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case if the source is an anisotropic fireball [222 , 223] or an 
electromagnetic accelerator with a modest beam width. 

Beaming requires, of course, an event rate larger by a 
ratio Att/6 2 compared to the observed rate. Observations 
of about one burst per 10 -6 year per galaxy implies one 
event per hundred years per galaxy if sa 7 -1 with 7 given 
by the compactness limit of ~ 100. 

5. AN OVERVIEW OF THE GENERIC MODEL 

It is worthwhile to summarize now the essential features 
of the generic GRB model that arose from the previous 
discussion. Compactness has led us to the requirement of 
relativistic motion, with a Lorentz factor 7 > 100. Ock- 
ham's razor and the desire to limit the total energy have 
lead us to the idea that the observed gamma-rays arise in 
the process of slowing down of a relativistic energy flow, 
at a stage that the motion of the emitting particles is still 
highly relativistic. 

This leads us to the generic picture mentioned earlier 
and to the suggestion that GRBs are composed of a three 
stage phenomenon: (i) a compact inner hidden "engine" 
that produces a relativistic energy flow, (ii) the energy 
transport stage and (iii) the conversion of this energy to 
the observed prompt radiation. One may add a forth stage 
(iv) conversion of the remaining energy to radiation in 
other wavelengths and on a longer time scale - the "after- 
glow" . 

5.1. Models for The Energy Flow 

The simplest mode of relativistic energy flow is in the 
form of kinetic energy of relativistic particles. A vari- 
ant that have been suggested is based on the possibility 
that a fraction of the energy is carried by Poynting flux 
[ 2241 , PM |H^> H> although in all models the power 
must be converted to kinetic energy somewhere. The en- 
ergy flow of ~ 10 50 ergs/sec from a compact object whose 
size is <10 7 cm requires a magnetic field of 10 15 Gauss or 
higher at the source. This large value might be reached 
in stellar collapses of highly magnetized stars or ampli- 
fied from smaller fields magnethohydrodynamically fl69| . 
Overall the different models can be characterized by two 
parameters: the ratio of the kinetic energy flux to the 
Poynting flux and the location of the energy conversion 
stage (ss 10 12 cm for internal conversion or « 10 16 cm for 



External shocks arise from the interaction of the shell 
with external matter. The typical length scale is the Sedov 



length, I = (E/r 



external conversion). This is summarized in Table 5.1. In 
the following section we will focus on the simplest possi- 
bility, that is of a kinetic energy flux. 

5.2. Models for The Energy Conversion 

Within the baryonic model the energy transport is in 
the from of the kinetic energy of a shell of relativistic par- 
ticles with a width A. The kinetic energy is converted to 
"thermal" energy of relativistic particles via shocks. These 
particles then release this energy and produce the observed 
radiation. There are two modes of energy conversion (i) 
External shocks, which are due to interaction with an ex- 
ternal medium like the ISM. (ii) Internal shocks that arise 
due to shocks within the flow when fast moving particles 
catch up with slower ones. Similar division to external and 
internal energy conversion occurs within other models for 
the energy flow. 



c 2)i/3. Therest 

mass energy within 
a sphere of radius I, equals the energy of the shell. Typi- 
cally I ~ 10 18 cm. As we see later (see section |8.7.1 ) rela- 
tivistic external shocks (with a Newtonian reverse shock) 
convert a significant fraction of their kinetic energy at 
i? 7 = Z/7 2 / 3 10 15 — 10 16 cm, where the external mass 
encountered equals 7" 1 of the shell's mass. Relativistic 
shocks (with a relativistic reverse shock) convert their en- 
ergy at i?A = i^A 1 / 4 « 10 16 cm, where the shock crosses 
the shell. 

Internal shocks occur when one shell overtakes another. 
If the initial separation between the shells is 8 and both 
move with a Lorentz factor 7 with a difference of order 
7 these shocks take place at: J7 2 . A typical value is 
10 " 
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5.3. Typical Radii 

In table |5.3| we list the different radii that arise in the 
fireball evolution. 



Figs. OandjT^ (from [228]) depict a numerical solution 
of a fireball from its initial configuration at rest to its final 
Sedov phase. 



6. FIREBALLS 

Before turning to the question of how is the kinetic en- 
ergy of the relativistic flow converted to radiation we ask 
is it possible to produce the needed flows? More specif- 
ically, is it possible to accelerate particles to relativistic 
velocities? It is remarkable that a relativistic particle flow 
is almost the unavoidable outcome of a "fireball" - a large 
concentration of energy (radiation) in a small region of 
space in which there are relatively few baryons. The rela- 
tivistic fireball model was proposed by Goodman [217] and 
by Paczynski |5^]. They have shown that the sudden re- 
lease of a large quantity of gamma ray photons into a com- 
pact region can lead to an opaque photon-lepton "fireball" 
through the production of electron-positron pairs. The 
term "fireball" refers here to an opaque radiation-plasma 
whose initial ener gy i s significantly greater than its rest 
mass. Goodman [217] considered the sudden release of a 
large amount of energy, E, in a small volume, character- 
ized by a radius, Ri. Such a situation could occur in an 
explosion. Paczynski K3| considered a steady radiation 
and electron-positron plasma wind that emerges from a 
compact region of size Ri with an energy, E, released on a 
time scale significantly larger than Ri / c. Such a situation 
could occur if there is a continuous source that operates 
for a while. As it will become clear later both configu- 
rations display, in spite of the seemingly large difference 
betw een t hem, the same qualitative behavior. Both Good- 
man 1 217 1 and Paczynski [[S3| considered a pure radiation 
fireball s in which there are no baryons. Later Shemi & 
Piran [ 220 1 and Paczynski [221] considered the effect of 
a baryonic load. They showed that under most circum- 
stances the ultimate outcome of a fireball with a baryonic 
load will be the transfer of all the energy of the fireball 
to the kinetic energy of the baryons. If the baryonic load 
is sufficiently small the baryons will be accelerated to a 
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Table 2: General Scheme for Energy Transport 
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Table 3: Critical Radii (a) - adiabatic fireball; (b) - radiative fireball 



relativistic velocity with 7 Ri E/M. If it is large the net when T = T v Ri 20keV. The matter opacity, 77,, on the 



result will be a Newtonian flow with v ~ y/2E/M. 

6.1. A simple model 

The evolution of a homogeneous fireball can be under- 
stood by a simple analogy to the Early Universe [220]. 
Consider, first, a pure radiation fireball. If the initial tem- 
perature is high enough pairs will form. Because of the 
opacity due to pairs, the radiation cannot escape. The 
pairs-radiation plasma behaves like a perfect fluid with an 
equation of state p = p/3. The fluid expands under of its 
own pressure. As it expands it cools with T cx (T be- 
ing the local temperature and R the radius). The system 
resembles quite well a part of a Milne Universe in which 
gravity is ignored. As the temperature drops below the 
pair-production threshold the pairs annihilate. When the 
local temperature is around 20keV the number of pairs be- 
comes sufficiently small, the plasma becomes transparent 
and the photons escape freely to infinity. In the meantime 
the fireball was accelerated and it is expanding relativisti- 
cally outwards. Energy conservation (as viewed from the 
observer frame) requires that the Lorentz factor that corre- 
sponds to this outward motion satisfies 7 oc R. The escap- 
ing photons, whose local energy (relative to the fireball's 
rest frame) is ~ 20keV are blue shifted. An observer at 
rest detects them with a temperature of T {> s oc 7T. Since 
T cx R' 1 and 7 cx R we find that the observed temper- 
ature, T bsi approximately equals 7o, the initial tempera- 
ture. The observed spectrum, is however, almost thermal 
1 217] and it is still far from the one observed in GRBs. 

In addition to radiation and e + e~ pairs, astrophysical 
fireballs may also include some baryonic matter which may 
be injected with the original radiation or may be present 
in an atmosphere surrounding the initial explosion. These 
baryons can influence the fireball evolution in two ways. 
The electrons associated with this matter increase the 
opacity, delaying the escape of radiation. Initially, when 
the local temper atur e T is large, the opacity is dominated 



other hand decreases only as R~ 2 , where R is the radius 
of the fireball. If at the point where t v = 1, 77, is still > 1, 
then the final transition to r = 1 is delayed and occurs at 
a cooler temperature. 

More importantly, the baryons are accelerated with the 
rest of the fireball and convert part of the radiation en- 
ergy into bulk kinetic energy. The expanding fireball has 
two basic phases: a radiation dominated phase and a mat- 
ter dominated phase. Initially, during the radiation domi- 
nated phase the fluid accelerates with 7 cx R. The fireball 
is roughly homogeneous in its local rest frame but due to 
the Lorentz contraction its width in the observer frame is 
A w Ri, the initial size of the fireball. A transition to the 
matter dominated phase takes place when the fireball has 
a size 
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by e + e pairs |217|. This opacity, r p , decreases expo- 



nentially with decreasing temperature, and falls to unity 



and the mean Lorentz factor of the fireball is 7 Ri E/Mc 2 . 
We have defined here E52 = £?/10 52 ergs and R47 = 
i?i/10 7 cm. After that, all the energy is in the kinetic en- 
ergy of the matter, and the matter coasts asymptotically 
with a constant Lorentz factor. 

The matter dominated phase is itself further divided 
into two sub-phases. At first, there is a frozen-coasting 
phase in which the fireball expands as a shell of fixed ra- 
dial width in its own local frame, with a width ~ 7 iij ~ 
(E/Mc 2 )Ri. Because of Lorentz contraction the shell ap- 
pears to an observer with a width A » K,. Eventu- 
ally, when the size of the fireball reaches R s = A7 2 Ri 
10 n cm(A/10 7 cm)(7/100) 2 variability in 7 within the fire- 
ball results in a spreading of the fireball which enters the 
coasting-expanding phase. In this final phase, the width 
of the shell grows linearly with the size of the shell, R: 

A^^/T^IOW^)^)" 2 (13) 
for ^>^ = 10 11 cm( I5 |-)( I l 5 ) 2 . 
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The initial energy to mass ratio, 77 = (E/Mc 2 ), deter- 
mines the order of these transitions. There are two critical 



values for 77 [ 220 1 
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and 



3(i T E 1/3 5 1/3 2/3 .--x 

These correspond to four different types of fireballs: 

• (i) A Pure Radiation Fireball (r) pa i r < 77): The 
effect of the baryons is negligible and the evolution is 
of a pure photon-lcpton fireball. When the temper- 
ature reaches T p , the pair opacity t p drops to 1 and 
Tb <C 1. At this point the fireball is radiation domi- 
nated (E > Mc 2 ) and most of the energy escapes as 
radiation. 

• (ii) Electron Dominated Opacity (77b < 77 < 
Vpair)'- In the late stages, the opacity is dominated 
by free electrons associated with the baryons. The 
comoving temperature decreases far below T p before 
t reaches unity. However, the fireball continues to 
be radiation dominated and most of the energy still 
escapes as radiation. 

• (iii) Relativistic Baryonic Fireball (1 < 77 < 77b): 
The fireball becomes matter dominated before it be- 
comes optically thin. Most of the initial energy is 
converted into bulk kinetic energy of the baryons, 
with a final Lorentz factor 7/ ss (E/Mc 2 ). This is 
the most interesting situation for GRBs. 

• (iv) Newtonian Fireball (77 < 1): This is the New- 
tonian regime. The rest energy exceeds the radiation 
energy and the expansion never becomes relativistic. 
This is the situation, for example in supernova explo- 
sions in which the energy is deposited into a massive 
envelope. 



6.2. Extreme- Relativistic Scaling Laws. 

The above summary describes the qualitative features 
of a roughly homogeneous expanding fireball. Surprisingly 
simi lar scaling laws exists also for inhomogeneous fireballs 
| 22S | as well as for relativistic winds ]5^] . Consider a spher- 
ical fireball with an arbitrary radial distribution of radi- 
ation and matter. Under optically thick conditions the 
radiation and the relativistic leptons (with energy density 
e) and the matter (with baryon mass density p) at each 
radius behave like a single fluid, moving with the same 
velocity. The pressure, p, and the energy density, e, are 
related by p — e/3. We can express the relativistic conser- 
vation equations of baryon number, energy and momen- 
tum using characteristic coordinates: r and s = t — r as 
|22fj: 

' \r 2 P u)=°(^-), (16) 



r 2 dr 
r 2 dr ( ' 



ds \7 + u 



i) 



=3/4 



ds \7 + u 



(17) 



J_d_ 

r 2 dr 



d_ 

ds 



de 
ds 



7 + u 
_ de 
dr 



(18) 



where u = u r — \J ^i 2 — 1, and we use units in which c = 1 
and the mass of the particles m = 1. The derivative d/dr 
now refers to constant s, i.e. is calculated along a char- 
acteristic moving outward at the speed of light. After a 
short acceleration phase we expect that the motion of a 
fluid shell will become highly relativistic (7 ^> 1). If we 
restrict our attention to the evolution of the fireball from 
this point on, we may treat 7 as a small parameter and 
set 7 ss u, which is accurate to order 0(j~ 2 ). Then, under 
a wide range of conditions the quantities on the right-hand 
sides of Eqs. |l6fjl8l are significantly smaller than those on 
the left. When we neglect the right hand sides of Eqs. [l6|- 
18 the problem becomes effectively only r dependent. We 
obtain the following conservation laws for each fluid shell: 



r p^f = const. 



r 2 e 3 / 4 7 - 



const., r yp+ —el 7 = const.. 

. ( 19 ) 

A scaling solution that is valid in both the radiation- 
dominated and matter-dominated regimes, as well as in 
the transition zone in between, can be obtained by com- 
bining the conserved quantities in Eq. |l9| appropriately. 
Let to be the time and ro be the radius at which a fluid 
shell in the fireball first becomes ultra-relativistic, with 
7 > few. We label various properties of the shell at this 
time by a subscript 0, e.g. 70, po, and eg. Defining the 
auxiliary quantity D, where 
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we find that 
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These are parametric relations which give r, p, and e of 
each fluid shell at any time in terms of the 7 of the shell 
at that time. The parametric solution ^l] describes both 
the radiation-dominated and matter-dominated phases of 
the fireball within the frozen pulse approximation. That 
is as long as the fireball does not spread due to variation 
in the velocity. 

6.3. The Radiation- Dominated Phase 

Initially the fireball is radiation dominated. For 7 <C 
(eo/po)7o> the first term in Eqj2FJ dominates and we find 
D oc r, 7 oc r, recovering the radiation-dominated scaling: 



7 oc r, 



p oc r 



ear 



(22) 



The scalings of p and e given in Eq. g2| correspond to those 
of a fluid expanding uniformly in the comoving frame. In- 
deed, all three scalings in Eq. |2^ can be derived for a 
homogeneous radiation dominated fireball by noting the 
analogy with an expanding universe. 
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Type 


r/ = hj/Mc 


M 


Pure Radiation 


^Ipair ^ T) 


M < M pair = 10- 12 M Q E^ 2 R]f 2 


Electrons Opacity 


T]b <V < Vpair 


M pair <M < M fc = 2 • 1CT 7 M e Eli 3 R 2 J 3 


Relativistic Baryons 


1 < V < Vb 


M b < M < 5 • 10- a M Q £: 5 2 


Newtonian 


r\ < 1 


5 • l(r 4 M Q £:5 2 < M 



Table 4: Different Fireballs 



Although the fluid is approximately homogeneous in its 
own frame, because of Lorentz contraction it appears as 
a narrow shell in the observer frame, with a radial width 
given by Ar ~ r/7 ~ constant ~ Ri, the initial radius 
of the fireball, or the initial width of the specific shell un- 
der discussion when we consider a continuous wind. We 
can now go back to Eqs. |l6Hl^ and set d/ds ~ j/r. We 
then find that the terms we neglected on the right hand 
sides of these equations are smaller than the terms on the 
left by a factor ~ 1 Ay. Therefore, the conservation laws 
19 and the scalings £2] are valid so long as the radiation- 
dominated fireball expands ultra-relativistically with large 
7. The only possible exception is in the very outermost 
layers of the fireball where the pressure gradient may be 
extremely steep and d/ds may be 3> 7/r. Ignoring this 
minor deviation, we interpret Eq. [l^ and the constancy 
of the radial width Ar in the observer frame to mean that 
the fireball behaves like a pulse of energy with a frozen 
radial profile, accelerating outward at almost the speed of 
light. 

6.4. The Matter-Dominated Phase 

The radiation dominated regime extends out to a radius 
r ~ (eo/po) r o- At larger radii, the first and last terms in 
Eq. EG become comparable and 7 tends to its asymptotic 
value of 7/ = (4eo/3po + l)7o- This is the matter dom- 
inated regime. (The transition occurs when 4e/3 = p, 
which happens when 7 = jj/2.) In this regime, D cx r 2 / 3 , 
leading to the scalings: 



7 



constant, 



p oc r 



,-8/3 



(23) 



The modified scalings of p and e arise because the fireball 
now moves with a constant radial width in the comoving 
frame. (The steeper fall-off of e with r is because of the 
work done by the radiation through tangential expansion.) 
Moreover, since e « /), the radiation has no important 
dynamical effect on the motion and produces no signifi- 
cant radial acceleration. Therefore, 7 remains constant on 
streamlines and the fluid coasts with a constant asymp- 
totic radial velocity. Of course, since each shell moves with 
a velocity that is slightly less than c and that is different 
from one shell to the next, the frozen pulse approximation 
on which Eq. |l^ is based must ultimately break down at 
some large radius. 

6.5. Spreading 

At very late times in the matter-dominated phase the 
frozen pulse approximation begins to break down. In this 
stage the radiation density e is much smaller than the mat- 
ter density p, and the Lorentz factor, 7, tends to a constant 
value 7/ for each shell. We may therefore neglect the term 



-(l/3)(<9e/dr) in Eq. |18| and treat 7 and u in Eqs. |1J-|18| 
as constants. We then find that the flow moves strictly 
along the characteristic, 0ft — r = constant, so that each 
fluid shell coasts at a constant radial speed, /?/ = Uf/jf. 
We label the baryonic shells in the fireball by a Lagrangian 
coordinate i?, moving with a fixed Lorentz factor 7y(_R), 
and let t c and r c represent the time and radius at which 
the coasting phase begins, which corresponds essentially 
to the point at which the fluid makes the transition from 
being radiation dominated to matter dominated. We then 
find 



r(t,R)-r c (R) 



if (R) 



2 7 2 (R) 



-(t-t c (R)) 



[t-t c (r)}. 



(24) 



The separation between two neighboring shells separated 
by a Lagrangian distance AR varies during the coasting 
phase as 



(25) 



~d{dr/dR) 


AR = 


1 c>7/ 


AR 


dt 




7/(-R) 3 OR 





Thus the width of the pulse at time t is Ar(i) 



A lf (t 



tc)/r} 



Ri + (t- t c )/jj, where Ar 



Ri 



is the width of the fireball when it begins coasting, 7/ is 
the mean 7/ in the pulse, and A7/ ~ 7/ is the spread of 
7/ across the pulse. From this result we see that within 
the matter dominated coasting phase there are two sep- 
arate regimes. So long as t — t c < jjRi, we have a 
frozen-coasting phase in which Ar is approximately con- 
stant and the frozen pulse approximation is valid. In this 
regime the scalings in Eq. 2^ are satisfied. However, 
when t — t c > jjRi, the fireball switches to an expanding- 
coasting phase where Ar cx t — t c and the pulse width 
grows linearly with time. In this regime the scaling of p 
reverts to p oc r -3 , and, if the radiation is still coupled to 
the matter, e oc r~ 4 . 

6.6. Optical Depth 

Independently of the above considerations, at some 
point during the expansion, the fireball will become op- 
tically thin. For a pure fireball this happens when the 
local temperature drops to about 20keV at: 



Rpai 



10 10 cm El^R^ 



(26) 



In a matter dominated fireball the optical depth is usu- 
ally determined by the ambient electrons. In this case the 
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fireball becomes optically thin at: 



cttE 



a/2 



iC 2 r] J 



6 • 10 



13 



-(77/100)" 1 . (27) 



From this stage on the radiation and the baryons no 
longer move with the same velocity and the radiation pres- 
sure vanishes, leading to a breakdown of Eqs. |il|]l§| . Any 
remaining radiation will escape freely now. The baryon 
shells will coast with their own individual velocities. If 
the fireball is already in the matter dominated coasting 
phase there will be no change in the propagation of the 
baryons. However, if the fireball is in the radiation dom- 
inated phase when it becomes optically thin, then the 
baryons will switch immediately to a coasting phase. This 
transition radius, R e has another crucial role in the fire- 
ball evolution. It is the minimal radius in which energy 
conversion and generation of the observed GRB can begin. 
Photons produced at R < R e cannot escape. 

6.7. Anisotropic Fireballs 

It is unlikely that a realistic fireball will be spherically 
symmetric. In fact strong deviation from spherical sym- 
metry are expected in the most promising neutron star 
merger model, in which the radiation is expected to emerge 
through funnels along the rotation axis. The initial motion 
of the fireball might be fairly complex but once 7> 1 the 
motion of each fluid element decouples from the motion of 
its neighbors with angular distance larger than 7 . This 
motion can be described by the same asymptotic solution, 
as if it is a part of a spherical shell. We define the angular 
range over which different quantities vary as 0. Addition- 
ally now the motion is not radially outwards and u r 7^ u. 
We define the spread angle <f> as u r = ucos(f>. The spheri- 
cal fireball equations hold locally if: 



and 



7- 



(28) 



7. 



TEMPORAL STRUCTURE AND KINEMATIC 
CONSIDERATIONS 



General kinematic considerations impose constraints on 
the temporal structure produced when the energy of a rela- 
tivistic shell is converted to radiation. The enormous vari- 
ability of the temporal profiles of GRBs from one burst to 
another in contrast to the relatively regular spectral char- 
acteristics, was probably the reason that until recently this 
aspect of GRBs was largely ignored. However, it turns out 
that the observed temporal structure sets a stro ng con- 
straint on the energy conversio n mo dels p0| , 23C]. GRBs 
are highly variable (see section 2.2) and some configura- 
tions cannot produce such temporal profiles. 

7.1. Time-scales 

Special relativistic effects determine the observed dura- 
tion of the burst from a relativistic shell (see Fig. fbi]) . 

• The Radial Time Scale: T ra< n a V- Consider an in- 
finitely thin relativistic shell with a Lorentz factor 
je (the subscript E is for the emitting region). Let 
Re be a typical radius characterizing the emitting 
region (in the observer frame) such that most of the 



emission takes place between Re and 2Re- The ob- 
served duration between the first photo n (e mitted at 
Re) and last one (emitted at 2Re) is ]205|, ul|: 



T 



.dial 



Re/2j 2 e c 



(29) 



• The Angular Time Scale: T ang : Because of rel- 
ativistic beaming an observer sees up to solid angle 
of 7 £ J 1 from the line of sight. Two photons emitted 
at the same time and radius Re, one on the line of 
sight and the other at an angle of 7^ away travel 
different distances to the observer. The diff erenc e 
lead to a delay in the arrival time by [205, |l8|, [230f : 



Tang = Re/2J E C 



(30) 



Clearly this delay is relevant only if the angular 
width of the emitting region, 9 is larger than 7^ . 

In addition there are two other time scales that are de- 
termined by the flow of the relativistic particles. These 
are: 

• Intrinsic Duration: AT: The duration of the flow. 
This is simply the time in which the source that 
produces the relativistic flow is active. AT = A/c, 
where A is the width of the relativistic wind (mea- 
sured in the observer's rest frame) . For an explosive 
source A « R4. However, A could be much larger 
for a wind. The observed duration of the burst must 
be longer or equal to A/c 

• Intrinsic Variability 8T: The time scale on which 
the inner source varies and produces a subsequent 
variability with a length scale <5 = cST in the flow. 
Naturally, ST sets a lower limit to the variability 
time scale observed in any burst. 



Clearly A and 8 must satisfy: 



R, 



< 6 < A 



(31) 



Finally we have to consider the cooling time scale. 



• The Cooling Time Scale: T coo i This is the dif- 
ference in arrival time of photons while the shocked 
material cools measured in the rest frame of an ob- 
server at rest at infinity. It is related to the local 
cooling time, e/P (where e is the internal energy 
density and P is the power radiated per unit vol- 
ume) in the fluid's rest frame by: 



Tcooi = e/PjE 



(32) 



Note that this differs from the usual time dilation 
which gives 7£e/P. 

For synchrotron cooling there is a unique energy 
dependence of th e coo ling time scale on frequency: 
T coo i{v) ex 1/-V2 Jhk| (sec Eq. @). If T cool deter- 
mines the variability we will have ST(v) cx v^ 1 / 2 . 
This is rema rkab ly close to the observed relation: 
ST cx v~ QA 1 102]. Quite generally Tc 00 i i s sh orter 



than the hydrodynamics time scales [103, 232, 
However, during the late stages of an afterglow, T coo i 
becomes the longest time scale in the system. 
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7.2. Angular Spreading and External shocks 

|and|3(i| (using ARe<Re) reveals 



Comparison of Eqs 
that T ang ~ T 
is larger than 7 
arisen due to irregularities in the properties of the shell or 
in the material that it encounters will be spread on a time 
given by T ang . This means that T ang is the minimal time 
scale for the observed temporal variability: ST > T ang . 

Comparison with the intrinsic time scales yields two 
cases: 



radial- As long as the shell's angular width 
1 , any temporal structure that could have 



rn _ J T ang = Re/cj 
-\A/c 



if A <R E /ll (Type-I); 
otherwise (Type-II) . 



(33) 

In Type-I models, the duration of the burst is deter- 
mined by the emission radius and the Lorentz factor. It is 
independent of A. This type of model s inc lude the stan- 
dard "external shock model" ^ EOI HI in which the 
relativistic shell is dec elera ted by the ISM, the relativistic 
magnetic wind model [226] in which a magnetic Poynting 
flux runs into th e ISM, or the scattering of star light by a 
relativistic shell HJ, g3i| . 

In Type-II models, the duration of the burst is deter- 
mined by the thickness of the relativistic shell, A (that 
is by the duration that the source is active and produces 
the relativistic wind). The angular spreading time (which 
depends on the the radius of emission) is shorter and there- 
fore irrelevant. These models include the "internal shock 
model" |2^, E29L ||(| , in which different parts of the shell are 
moving with different Lorentz factor and therefore collide 
with one anot her. A magnetic dominated version is given 
by Thompson |224 |. 

The majority o f G RBs have a complex temporal struc- 
ture (e.g. section |J) with Af = T/5T of order 100. Con- 
sider a Type-I model. Angular spreading means that at 
any given moment the observer sees a whole region of an- 
gular width 7^ . Any variability in the emission due to dif- 
ferent conditions in different radii on a time scale smaller 
than T ang is erased unless the angular size of the emit- 
ting region is smaller than 7^. Thus, such a source can 
produce only a smooth single humped burst with Af = 1 
and no temporal structure on a time-scale ST < T. Put 
in other words a shell, of a Type-I model, and with an 
angular width larger than 7^ cannot produce a variable 
burst with Af 1 . This is the angular spreading problem. 

On the other hand a Type-II model contains a thick shell 
A > Re He anc ^ ^ can produce a variable burst. The vari- 
ability time scale, is again limited ST > T ang but now it 
can be shorter than the total duration T . The duration of 
the burst reflects the time that the "inner engine" oper- 
ates. The variability reflects the radial inhomogeneity of 
the shell which was produced by the source (or the cooling 
time if it is longer than S/c). The observed temporal vari- 
ability provides an upper limit to the scale of the radial 
inhomogeneities in the shell and to the scale in which the 
"inner engine" varies. This is a remarkable conclusion in 
view of the fact that the fireball hides the "inner engine" . 

Can an external shock give rise to a Type-II behavior? 
This would have been possible if we could set the parame- 
ters of the exte rnal s hock model to satisfy Re < 2jecST. 
As discussed in 8.7.1 the deceleration radius for a thin shell 



with an initial Lorentz factor 7 is given by 
Re = h- 2/ \ 



(34) 



and the observed duration is l~/~ 8 / 3 /c. The deceleration is 
gradual and the Lorentz factor of the emitting region je 
is similar to the original Lorentz factor of the shell 7. It 
seems that with an arbitrary large Lorentz factor 7 we can 
get a small enough deceleration radius Re- However , Eq . 
34 is valid only for thin shells satisfying A > l~/~ 8 / 3 ]233[ . 
As 7 increases above a critical value 7 > 7 C = (1/ A) 3 / 8 
the shell can no longer be considered thin. In this situa- 
tion the reverse shock penetrating the shell becomes ultra- 
relativistic and the shocked matter moves with Lorentz 
factor je — 7c < 7 which is independent of the initial 
Lorentz factor of the shell 7. The deceleration radius is 
now given by Re = A 1 / 4 ^ 3 / 4 , and it is independent of the 
initial Lorentz factor of the shell. The behavior of the de- 
celeration radius Re and observed duration as function of 
the shell Lorentz factor 7 is given in Fig. [l5] for a shell 
of thickness A = 3 x 10 12 cm. This emission radius Re is 



always larger than A/7^ 
be of type II. 



thus an external shock cannot 



7.3. Angular Variability and Other Caveats 

In a Type-I model, that is a for a shell satisfying 
A < Re/^e: variability is possible only if the emitting 
regions are significantly narrower than 7^ . The source 
would emit for a total duration T ra( iiai- To estimate the 
allowed opening angle of the emitting region imagine two 
points that emit radiation at the same (observer) time t. 
The difference in the arrival time between two photons 
emitted at (Re,&i) an d {Re, #2) at the same (observer) 
time t is: 



ST 



Re\0 2 



R e 0\ 



0i I R E ose 



(35) 



where ( 
used 61 



2c c c 

is the angle from the line of sight and we have 
h -C 1, 9 = (6»i + 2 )/2 and 89 = \0 2 - 9i\. Since 
an observer sees emitting regions up to an angle 7^ 1 away 
from the line of sight ~ Je ■ The size of the emitting 
region r s — ReS9 is limited by: 

r s < JecST . (36) 

The corresponding angular size is: 

JecST 1 



59 < 



Re 



Af-yi 



(37) 



Note that Fenimore, Madras and Nayakshin [ 230 1 who ex- 
amined this issue, considered only emitting regions that 

are directly on the line of sight with ~ |02 — 0i| and ob- 

1/2 

tained a larger r s which was proportional to i? B . How- 
ever only a small fraction of the emitting regions will be 
exactly on_the line of sight. Most of the emitting regions 
will have ~ 7^, and thus Eq. |3^ yields the relevant 
estimate. 

The above discussion suggests that one can produce 
GRBs with T » T radial « Re/cj% and ST = T/Af 
if the emitting regions have angular size smaller than 
1/AfjE ~ 10~ 4 . That is, one needs an extremely nar- 
row jet. Relativistic jets are observed in AGNs and even 
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in some galactic objects, however, their opening angles 
are of order of a few degrees almost two orders of magni- 
tude larger. A narrow jet with such a small opening angle 
would be able to produce the observed variability. Such a 
jet must be extremely cold (in its local rest frame); other- 
wise its internal pressure will cause it to spread. It is not 
clear what could produce such a jet. Additionally, for the 
temporal variability to be produced, either a rapid modu- 
lation of the jet or inhomogeneities in the ISM are needed. 
These two options are presented in Fig. [ll| 

A second possibility is that the shell is relatively "wide" 
(wider than 7 £ J 1 ) but the emitting regions are small. An 
example of this situation is schematically described in Fig. 
|l7j This may occur if either the ISM or the shell itself 
are very irregular. This situation is, however, extremely 
inefficient. The area of the observed part of the shell is 
nRg/jg. The emitting regions are much smaller and to 
comply with the temporal constraint their area is irr 2 . For 
high efficiency all the area of the shell must eventually ra- 
diate. The number of emitting regions needed to cover 
the shell is at least (Re/jet s ) 2 ■ In Type-I models, the 
relation Re = 2cjgT holds, and the number of emitting 
region required is 4AT 2 . But a sum of 4J\f 2 peaks each of 
width 1 /Af of the total duration does not produce a com- 
plex time structure. Instead it produces a smooth time 
profile with small variations, of order \j\JlM <C 1, in the 
amplitude. 

In a highly variable burst there cannot be more than 
M sub-bursts of duration ST = T/Af. The corresponding 
area covering factor (the fraction of radiating area of the 
shell) and the corresponding efficiency is less than X/4AT. 
This result is independent of the nature of the emitting re- 
gions: ISM clouds, star light or fragments of the shell. This 
is the case, for example, in the Shaviv & Dar model [235 
where a relativistic iron shell interacts with the starlight 
of a stellar cluster (a spherical shell interacting with an 
external fragmented medium). This low efficiency poses a 
series energy crisis for most (if not all) cosmologic al m odels 
of this kind. In a recent paper Fenimore et al. [ 231 1 con- 
sider other ways, which are based on low surface covering 
factor, to resolve the angular spreading problems. None 
seems very promising. 

7.4. Temporal Structure in Internal shocks. 

Type-II behavior arises naturally in the internal shock 
model. In this model different shells have different Lorcntz 
factors. These shells collide with each other and con- 
vert some of their kinetic energy into internal energy 
and radiate (Fig. |l8|). If the emission radius is suffi- 
ciently small angular spreading will not erase the tem- 
poral variability. This requires Re < 2j^,cST. This con- 
dition is always satisfied as internal shocks take place at 
Re = R s ~ <^7 2 with je — 7- Since S < A we have 
T = A/c > T ang = R E /2j%c = 5/c = ST. Clearly multi- 
ple shells are needed to account for the observed temporal 
structure. Each shell produces an observed peak of du- 
ration ST and the whole complex of shells (whose width 
is A) produces a burst that lasts T = A/c. The angular 
spreading time is comparable to the temporal variability 
produced by the "inner engine" . They determine the ob- 
served temporal structure provided that they are longer 
than the cooling time. 



Before concluding that internal shocks can actually pro- 
duce GRBs we must address two issues. First, can internal 
shocks produce the observed variable structure? Second, 
can it be done efficiently? We will addr ess the first ques- 
tion here and the second one in section B.6 where we dis- 



cuss energy conversion in internal shocks. 

Mochkovitch, Maitia and Marques (3^] and Kobayashi, 
Piran & Sari |33| used a simple model to calculate the tem- 
poral profiles produce by an internal shock. According to 
this model the relativistic flow that produces the shocks is 
characterized by multiple shells, each of width I and with 
a separation L. The shells are assigned random Lorentz 
factors (varying in the range [7mm,7mra]) and random 
density or energy. One can follow the motion of the shell 
and calculate the time of the binary (two-shell) collisions 
that take place, until all the shells that could collide have 
collided and the remaining flow has a monotonically de- 
creasing velocity. The energy generated and the emitted 
radiation in each binary encounter are then combined to 
a synthetic sample of a temporal profile. 

The emitted radiation from each binary collision will be 
observed as a single pulse characterized by an amplitude 
and a width. The amplitude depends on the amou nt o f 
energy converted to radiation in a given collision (see B.6). 
The time scale depends on the cooling time, the hydrody- 
namic time, and the angular spreading time. The internal 
energy is radiated via synchrotron emission and inverse 
Compton scattering. In most of cases, the electrons' cool- 
ing t i me i s much shorter than the hydrodynamic time scale 



|10||232 



It 



so we consider only the latter two. 
The hydrodynamic time scale is determined by the time 
that the shock crosses the shell, whose width is d. In fact 
there are two shocks: a forward shock and a reverse shock. 
Detailed calculations |53| reveals that this time scale (in 
the observer's rest frame) is of order of the light crossing 
time of the shell, that is: d/c. 

Let the distance between the shells be S. A collision 
takes place at 5/j 2 and the angular spreading yields an 
observed pulse whose width is ~ 5/c. If S > d the overall 
pulse width 5T is determined by angular spreading. The 
shape of the pulse become asymmetric with a fast rise and 
a slower decline (Fig. 1^) which GRBs typically show (see 
section 2.2). The amplitude of an individual pulse depends 
on the energy pr odu ced in the collision, which we calculate 
latter in section 8.6. 

Typical synthetic temporal profile are shown in Fig. 
Clearly, internal shocks can produce the observed highly 
variable temporal profiles, provided that the source of the 
relativistic particles, the "inner engine" , produces a suffi- 
ciently variable wind. Somewhat surprisingly the result- 
ing temporal profile follows, to a large extent, the shape 
of the pulse emitted by the source. Long bursts require 
long relativistic winds that last hundred seconds, with a 
rapid variability on a time scale o f a fr action of a second 
Thus, unlike previous worries J7S" 
is some direct information on the 



236] we find that there 
; 'inncr engine" of GRBs: 
It produces the observed complicated temporal structure. 
This severely constrains numerous models. 



8. ENERGY CONVERSION 
8.1. Slowing Down of Relativistic Particles 
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The cross section for a direct electromagnetic or nu- 
clear interaction between the relativistic particles and the 
ISM particles is too small to convert efficiently the kinetic 
energy to radiation. The fireball particles can be slowed 
down only via some collective interaction such as a colli- 
sionless shock. Supernova remnants (SNRs) in which the 
supernova ejecta is slowed down by the ISM show that col- 
lisionless shocks form in somewhat similar circumstances. 
One can expect that collisionless shocks will form here as 
well §7[ @. 

GRBs are the relativistic analogues of SNRs. In both 
cases the phenomenon results from the conversion of the 
kinetic energy of the ejected matter to radiation. Even 
the total energy involved is comparable. One crucial dif- 
ference is the amount of ejected mass. SNRs involve a 
solar mass or more. The corresponding velocity is several 
thousands kilometers per second, much less than the speed 
of light. In GRBs, the masses are smaller by several orders 
of magnitude and with the same energy the matter attains 
ultra-relativistic velocities. A second crucial difference is 
that while SNRs result from the interaction of the ejecta 
with the ISM, GRBs result from internal collisions. The 
interaction of the ejecta in GRBs with the ISM produces 
the "afterglow" that follows some GRBs. 

The interaction between the SNR ejecta and the ISM 
takes place on scales of several pc and it is observed for 
thousands of years. The internal interaction of the rela- 
tivistic matter in GRBs takes place on a scale of several 
hundred astronomical units and special relativistic effects 
reduce the observed time scale to a fraction of a second. 
The interaction with the ISM that leads to the "afterglow" 
takes place on a scale of a tenth of a pc. Once more special 
relativistic effects reduce the observed time scale to several 
days. 

In the following sections I discuss the slowing down of 
matter due to relativistic shocks. The discussion begins 
with a general review of relativistic inelastic collisions and 
continues with the relativistic shock conditions. Then I 
turn to the radiation processes: synchrotron emission and 
Inverse Compton. After the general discussion I apply the 
general results to internal shocks, to external shocks and 
to the afterglow. 



8.1.1. Relativistic Inelastic Collisions 

Consider a mass (denoted by the subscript r) that 
catches up a slower one (denoted s). The two masses col- 
lide and merge to form a single mass (denoted m) . Energy 
and momentum conservation yield: 



is needed to yield 7 m ~ 7r /2 and £ ~ m r /2. The external 
mass needed to convert half of the kinetic energy is smaller 



r 7 r + m s 7 s = {m r + m s + £/c 2 )y m (38) 



•\/lr - 1 + m s Vis - 1 = ( m r + m s +£/c 2 )V% 



1. 



where £ is the internal energy generated in the collision 
(in the rest frame of the merged mass). 

There are two interesting limits. First let m s be at rest: 
7 S = 1. This is the case in external shocks, or in a shock 
between relativistic ejecta and a non-relativistic material 
that was ejected from the source before it exploded. Eqs. 
39 reveal that a mass: 




than the original mass by a factor of 7 r 

The second case corresponds to an internal collision be- 
tween shells that are moving at different relativistic veloc- 
ities: 7r^7s > 1. Eqs. p9| yield: 



7n 



(40) 



The internal energy (in the frame of an external observer) 
of the merged shell, Ei nt = r y m £, is the difference of the 
kinetic energies before and after the collision: 



E int = m r c 2 ("f r - 7m) + m s c 2 (y s - 7 m ). 



(41) 



The conversion efficiency of kinetic energy into internal 
energy is [p|, @: 



1 



(m r +m s )7 m 
(m r ^ r + m s j s ) 



(42) 



As can be expected a conversion of a significant fraction of 
the initial kinetic energy to internal energy requires that 
the difference in velocities between the shells will be signif- 
icant: 7 r 3> 7 S and that the two masses will be comparable 



8.1.2. Shock Conditions 

Quite often the previous estimates based on the approx- 
imation that the whole shell behaves as a single object are 
good enough. However, the time scale between the inter- 
action of different parts of the shell with the ISM may be 
relatively long (compared to the time scale to collect an 
external mass M/7) and in this case one has to turn to the 
hydrodynamics of the interaction. This calculation takes 
into account the shocks that form during the collision. 

Consider the situation a cold shell (whose internal en- 
ergy is negligible compared to the rest mass) that overtakes 
another cold shell or moves into the cold ISM. Generally, 
Two shocks form: an outgoing shock that propagates into 
the ISM or into the external shell, and a reverse shock 
that propagates into the inner shell, with a contact dis- 
continuity between the shocked material (see Fig. ^l|) . 
Two quantities determine the shocks' structure: 7, the 
Lorentz factor of the motion of the inner shell (denoted 
4) relative to the outer one - or the ISM (denoted 1) and 
the ratio between the particle number densities in these 
regions, / = n^jn\. 

There are three interesting cases: (i) Ultra-relativistic 
shock (7^1) with / > ■y 2 . This happens during the early 
phase of an external shock or during the very late external 
shock evolution when there is only a single shock. We call 
this configuration "Newtonian" because the reverse shock 
is non-relativistic (or mildly relativistic). In this case the 
energy conversion takes place in the forward shock: Let 72 
be the Lorentz factor of the motion of the shocked fluid 
relative to the rest frame of the fluid at 1 (an external ob- 
server for interaction with the ISM and the outer shell in 
case of internal collision) and let 73 be the Lorentz factor 
of the motion of this fluid relative to the rest frame of the 
relativistic shell (4). 



(39) 



72 ~ 7 



73 



1. 



(43) 
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The particle and energy densities (n, e) in the shocked re- 
gions satisfy: 

U2 « 47m, ; e = e 2 = 47 2 riim p c 2 ; n 3 = 7n 4 , ; e 3 = e. 

(44) 

(ii) Later during the propagation of the shell the den- 
sity ratio decreases and / < j 2 . Both the forward and 
the reverse shocks are relativis tic. The shoc k equations 
between regions 1 and 2 yield: [§3§ |39[ |||: 



72 = / 1/4 7 1/2 /V2 ; n 2 = A-y 2 ni ; e = e 2 = A^nim p c 2 , 
Similar relations hold for the reverse shock: 



73 = r 1/4 7 1/2 /V2 ; 7,3 = 473^4. 



(46) 



In addition we have e 3 = e and 73 = (7/72 + 72/t)/2 
which follow from the equality of pressures and velocity 
on the contact discontinuity. Comparable amounts of en- 
ergy are converted to thermal energy in both shocks when 
both shocks are relativistic. But only a negligible amount 
of energy is converted t o th ermal energy in the reverse 
shock if it is Newtonian p33[ . 

(iii) Internal shocks are characterized by / w 1 - both 
shells have similar densities, and by a Lorentz factor of 
order of a few (2 < 7 < 10) describing the relative motion 
of the shells. In this case, for an adiabatic index (4/3) we 
have: 



72 = V(7 2 + l)/2 « max[l, ; 
n- 2 = (47 2 + 3)ni « 47 2 ni ; e 2 = 7 2 n 2 m p c 2 



(47) 



Both shocks are mildly relativistic and their strength de- 
pends on the relative Lorentz factors of the two shells. 



8.1.3. Lorentz Factors in Different Emitting Regions 

Before concluding this section and turning to the radi- 
ation mechanisms we summarize briefly the different rel- 
ativistic motions encountered when considering different 
emitting regions. The relativistic shocks are characterized 
by Jsh that describes the shock's velocity as well as the 
"thermal" motion of the shocked particles. It is measured 
relative to a rest frame in which the unshocked material is 
at rest. The Lorentz factor of the forward shock is usually 
different from the Lorentz factor of the reverse shock. The 
emitting region - the shocked material - moves relativis- 
tically relative to an observer at rest at infini ty. T his is 
characterized by a Lorentz factor, je- Table 3.1.3 sum- 



marizes the different values of "fsh and je for external and 
internal shocks and for the afterglow. 



8.2. Synchrotron Emission from Relativistic Shocks 
8.2.1. General Considerations 



The most likely rad iation 
chrotron emission [24C, O , 



pr oces s 



222. L03 



in GRBs is syn- 
Thc observed low 



energy spect ra p rovide an indication that this is indeed 
the case (951 p41 |. 

The parameters that determine synchrotron emission 
are the magnetic field strength, B, and the electrons' en- 
ergy distribution (characterized by the minimal Lorentz 



factor j e ,rnin and the index of the expected power-law elec- 
tron energy distribution p). These parameters should be 
determined from the microscopic physical processes that 
take place in the shocks. However, it is difficult to es- 
timate them from first principles. Instead we define two 
dimcnsionlcss parameters, 6 r and e fi , that incorporate our 



103 | 



ignorance and uncertainties (222 , 

The dimensionless parameter eb measures the ratio of 
the magnetic field energy density to the total thermal en- 
ergy e: 

UB Dl (48) 



87re' 



so that, after substituting the shock conditions we have: 



B 



1/2 1/2 



(49) 



The r e ha ve been different attempts to estimate es 



240 , 226]. We keep it as a free parameter. Additionally 



we assume that the magnetic field is randomly oriented in 
space. 

The second parameter, e e , measures the fraction of the 
total thermal energy e which goes into random motions of 
the electrons: 

e e = ^ . (50) 



8.2.2. The Electron Distribution 

We call consider a "typical" electron as one that has the 
average 7 e of the electrons distribution: 



<7e 



HI,, 



(51) 



Collisionless acceleration of electrons can be efficient if 
they are tightly coupled to the proto ns and the magnetic 
field by mean of plasma waves [243]. Since the electrons 
receive their random motions through shock-heating, we 
assume (following the treatment of non-relativistic shocks) 
that they develop a power law distribution of Lorentz fac- 
tors: 



N (le) ~ 7e P for 7e > le 



(52) 



We require p > 2 so that the energy does not diverge at 
large j e . Since the shocks are relativistic we assume that 
all the electrons participate in the power-law, not just a 
small fraction in the tail of the distribution as in the New- 
tonian case. An indication that this assumption is correct 
is given by the lower energy spectrum observed in some 
GRBs f95|, |4l). The minimum Lorentz factor, r fe,mim 

of 

Dution is related to e e and to the total energy 



^5 
the distri 



y s hnm p c" 



7e 



TO p p - 



m e p - 



1 



Celsh 



p-1 



<7e 



(53) 



where 7^ is the relative Lorentz factor across the corre- 
sponding shock. 

The energy index p can be fixed by requiring that the 
model should be able to explain the high energy spectra of 
GRBs. If we assume that most of the radiation observed in 
the soft gamma-rays is due to synchrotron cooling, then it 
is straightforward to relate p to the power-law index of the 
observed spectra of GRBs, (3. The mean spectral index of 
GRBs at high photon energies f3 w —2.25, [pij corresponds 
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Shock type 


IE 




External 


Newtonian 


Forward 


7 


7 


Reverse 


7 


1 


Relativistic 


Forward 


7C 3/4 


7£ 3/4 


Reverse 


7C 3/4 




Internal 


7 


\/7mt/2 ~ a few 


Afterglow 




7(0 



Table 5: Lorentz Factors 



to p k 2.5. This agrees, as we see later (9.3.2) with the 
value inferred from afterglow observations (p ~ 2.25). We 
assume this value of p in what follows. The corresponding 
ratio that appears in Eq. 53 (p — 2)/(p — 1) equals 1/3 and 
we have 7 e , mm = 6107^. 

The shock acceleration mechanisms cannot accelerate 
the electrons to arbitrary high energy. For the maximal 
electron's energy, with a corresponding je.max, the accel- 
eration time equals to the cooling time. The accelera- 
tion time is determ ined by the Larmor radius Rl and the 
Alfven velocity va |244| : 



cRr 



(54) 



This time scale should be compared with the synchrotron 
cooling time 7 e m e c 2 / 'P syn (in the local frame). Using 



4, Eq. 



to estimate B and Eq. 511 below to 



estimate P syn one finds: 



247repe „ „ „ s < 
B - 3.7 x 10 8 



3/2 
B 



Bo~t1. 



sh 



1/2 

7sfc% 



(55) 



cooling This value is quite large and generally it does not 
effect the observed spectrum in the soft gamma ray range. 



8.2.3. Synchrotron Frequency and Synchrotron Power 

The typical energy of synchrotron photons as well as the 
synchrotron cooling time depend on the Lorentz factor j e 
of the relativistic electron under consider atio n and on the 
strength of the magnetic field (see e.g. [245]). Since the 



emitting material moves with a Lorentz factor je the pho- 
tons are blue shifted. The characteristic photon energy in 
the observer frame is given by: 



r , s hqeB 2 

{hv S yn)obs = —lelE- 



m e c 



(56) 



The power emitted by a single electron due to syn- 
chrotron radiation in the local frame is: 



P,, 



-OT&JBle , 



(57) 



where or is the Thomson cross section. The cooling time 
of the electron in the fluid frame is then 7 e m e c 2 / P . The 
observed cooling time t syn is shorter by a factor of 75: 



tsyn (7e) 



3m e c 



4cr T £/B7e7E 



(58) 



Substituting the value of 7, 
cooling rate Eq, 



from equation pq into the 



58 we obtain the cooling time scale as a 



function of the observed photon energy: 



3 2Trchm e q e 

tsyn\V) — 4 / — — ■= 



(59) 



Since 7 e does not appear explicitly in this equation t syn 
at a given observed frequency is independent of the elec- 
trons' energy distribution within the shock. This is pro- 
vided, of course, that there are electrons with the required 
7 e so that there will be emission in the frequency consid- 
ered. As long as there is such an electron the cooling time 
is "universal". This equation shows a characteristic scal- 
ing of t syn (y) cx v^ 1 ! 2 . This i s no t very different from the 
observed relation ST cx v~ 0A |l02}| . However, it is not clear 
if the cooling time and not another time scale determined 
the temporal profile. 

The cooling time calculated above sets a lower limit to 
the variability time scale of a GRB since the burst cannot 
possibly contain spikes that are shorter than its cooling 
time. Observations of GRBs typically show asymmetric 
spikes in the intensity variation, where a peak generally 
has a fast rise and a slower exponential decline (FRED). A 
plausible explanation of this observation is that the shock 
heating of the electrons happens rapidly (though episod- 
ically), and that the rise time of a spike is related to the 
heating time. The decay time is then set by the cooling, 
so that the widths of spikes directly measure the cooling 
time. 



8.2.4. The Integrated Synchrotron Spectrum 

The instantaneous synchrotron spectrum of a single elec- 
tron with an initial energy 7 e m e c 2 is a power law with 
F v cx v 1 !^ up to Vsynirie) and an exponential decay above 
it. If the electron is energetic it will cool rapidly until it 
will reach j e , c . This is the Lorentz factor of an electron 
that cools on a hydrodynamic time scale. For a rapidly 
cooling electron we have to consider the time integrated 
spectrum above h syn {^f e _ c ): F v cx u^ 1 / 2 from v S yn{le,c) up 

tO V syn {"ie)- 

To calculate the overall spectrum due to all the electrons 
we n eed to integrate over j e . Our discussion here follows 
|246( j . We consider a power-law electron distribution with 
a power index p and a minimal Lorentz factor 7e,mm 

(see 

Eq. |52| ). Overall we expect a broken power law spectrum 
with a break frequency around the synchrotron frequency 
of the lowest energy electrons y S yn{le,min)- These power 
law indices depend on the cooling rate. The most energetic 
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electrons will always be cooling rapidly (independently of 
the behavior of the "typical electron" ) . Thus the highest 
spectrum is always satisfy: 



F v = N[("f(v)]m e c 2 (~f(v)d-f/dv oc iy- p/2 . 



(60) 



similarly the low energy electrons will always be slow cool- 
ing and thus the lowest part of the spectrum will behave 
like F v oc v 1 ^. 

For slow cooling we have the instantaneous spectrum: 
F y oc v 1 ^ for the lower part of the spectrum. For the 
upper part we have 

F» = N[(j(u)]P[(j(u)]d'y/dp oc jy-^ 1 ^ 2 , (61) 

where "f(v) is the Lorentz factor for which the synchrotron 
frequency equals v. The most energetic electrons will cool 
rapidly even when the overall system is in slow cooling. 
These electrons emit practically all their energy m e c 2 7, at 
their synchrotron frequency. Thus the uppermost part of 
the spectrum will satisfy: 

For fast cooling we have F v oc v^ 1 / 2 for the lower part 
and F v oc v~ p / 2 for the upper part. Here at the lower end 
the least energetic electrons will be cooling slowly even 
when the typical electron is cooling rapidly. Thus we will 
have /„ oc v 1 ^ 3 in the lowest part of the spectrum. 

The critical parameter that determines if the electrons 
are cooling fast or slow is 7e iC , the Lorentz factor of an elec- 
tron that cools on a hydrodynamic time scale. To estimate 
7e, c we compare t syn (Eq. |58|) with thydi the hydrodynamic 
time scale (in the observer's rest frame): 



3m P c 



^(JTUBlEthyd 



(62) 



Fast cooling occurs if 7 e , c < 7e,rmn- All the electrons cool 
rapidly and the electrons' distribution effectively extends 
down to 7e, c . If 7e, c > 7e,mm only the high energy tail of 
the distribution (those electrons above 7 e . c ) cool and the 
system is in the slow cooling regime. 

For the GRB itself we must impose the condition of 
fast cooling: the relativistic shocks must emit their energy 
effectively - otherwise there will be a serious inefficiency 
problem. Additionally we won't be able to explain the 
variability if the cooling time is too long. The electrons 
must cool rap idly and release all their energy. In this case 
7e,min > 7e,c ]l03| and all the electrons cool down roughly 
to 7e, c - The observed flux, F v , is given by: 



F u oc { {v/v c )- x ' 2 F v 

.max-, ^rn ^> V ^> Vq. 

{VrnlVc)~ 1/2 (v/v m )~ P/2 F v ^nax-, V > V m , 



where V m = V S yn(le,min) 



y syn 



(7e, c ) and F v 



(63) 
is the 



observed peak flux. 

It is most likely that during the latter stages of an ex- 
ternal shock (that is within the afterglow phase - pro- 
vided that it arises due to external shocks) ther e wi ll be 
a transition from fast to slow cooling pi], |23|, ^7], |247| , p5[ . 
When 7 C > 7 e ,mm, only those electrons with j e > -f c can 
cool. We call this slow cooling, because the electrons with 
7e ~ le.min, which form the bulk of the population, do 



not cool. Integration over the electron distribution gives 
in this case: 

(^/t'm)~ (P ~ 1)/2 f 1 /,ma a; , V c > V > V m , 

(v c /v m )~ (p ~ 1)/2 (v/v c )~ P/2 F^ maXl V > V c . 

(64) 

For fast cooling v c < v m . We find that the peak flux is 
at v c while the peak energy emitted (which corresponds to 
the peak of vF v ) is at v m . For slow cooling the situation 
reverses v m < v c . The peak flux is at v m while the peak 
energy emitted is at v c . 

Typical spectra corresponding to fast and slow cooling 
are shown in Fig. ^2[ The light curve depends on the hy- 
drodynamic evolution, which in turn determines the time 
dependence of v mi v c and F UtTnax . 

For fast cooling the power emitted is simply the power 
given to the electrons, that is e e times the power generated 
by the shock: 

dE , . 

Pfast = e e-jj-- ( 65 ) 

For slow cooling the emitted power is determined by the 
ability of the electrons to radiate their energy: 



Pslow — ^eFsyn\fc,min} 



(66) 



where, N e is the number of electrons in the emitting region 
and P S yn(le,min), the synchrotron power of an electron 
with 7e iTO i„, is given by Eq. |7|. 

8.3. Synchrotron Self Absorption 

An important effect that we have ignored so far is the 
possibility of self absorption. This is irrelevant during the 
GRB itself. One of the essential features of the GRB spec- 
trum is that it is produced in an optically thin region. 
However, self absorption may app ear at la t e tim e and typ- 



ically in radio emission 1l8l, 247 



248 



249 1 . When it 



appears it will cause a steep cutoff of the low energy spec- 
trum, either as the commonly known v°l 2 or as v 2 . 

To estimate the self absorption frequency we need the 
optical depth along the line of sight. A simple approxi- 
mation is: a' v ,Ri/ji where a' v , is the absorption coefficient 
fl2lF 



(P + 2) 
%ixm e v 12 



(67) 



The self absorption frequency v a satisfies: a' v ,R/i = 1. It 
can be estimates only once we have a model for t he hydro - 
dynamics and how do R and 7 change with time [ [248| , |249| . 

The spectrum below the the self- absorption frequency 
depend s on the electron distribution. One obtaine the well 
known [ 245 |, v 5 ^ 2 slop when the synchrotron frequency of 
the electron emitting the self absorbed radiation is inside 
the self absorption range. One obtains a slop of v 2 2 if 
there is self absorption, but the radiation in that range is 
due to the low energy tail of electrons radiating effectively 
at higher energies. For this latter case, which is more ap- 
propriate for GRB afterglow we find that H, SHI : 



F v oc v 2 [k B T e /{ imp c 2 )]R 2 



(68) 
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where R is the radius of the radiating shell and the factor 
k B T e /(jm p c 2 ) describes the degree of electron equiparti- 
tion in the plasma shock-heated to an internal energy per 
particle 'ym p c 2 and moving with Lorentz factor 7. 

8.4. Inverse Compton Emission 

Inverse Compton (IC) scattering may modify our anal- 
ysis in several ways. IC can influence the spectrum even 
if the system is opti cally thin (as it must be) to Compton 
scattering (see e.g. |245fl ). In view of the high energies in- 
volved we assume that only one IC scattering takes place. 
After this scattering the photon's energy is so high that 
in the electron's rest frame it is above the Klein-Nishina 
energy and in this case the decrease in the Compton cross 
section makes this scattering unlikely. The effect of IC de- 
pends on the Comptonization parameter Y — 7 2 r e . For 
fast cooling one can show [ 103 1 that Y satisfies: 



Y 
Y 



if e e < e B 
if e e > e B . 



(69) 



IC is unimportant if Y < 1 and in this case it can be 
ignored. 

If Y > 1, which corresponds to e e > e B and to Y = 
\J e e /e B then a large fraction of the low energy synchrotron 
radiation will be up scattered by IC and a large fraction of 
the energy will be emitted via the IC processes. If those IC 
up scattered photons will be in the observed energy band 
then the observed radiation will be IC and not synchrotron 
photons. Those IC photons might be too energetic, that 
is their energy may be beyond the observed energy range. 
In this case IC will not influence the observed spectra di- 
rectly. However, as IC will take a significant fraction of the 
energy of the cooling electrons it will influence the obser- 
vations in two ways: it will shorten the cooling time (the 
emitting electrons will be cooled by both synchrotron and 
IC process). Second, assuming that the observed 7-ray 
photons results from synchrotron emission, IC will influ- 
ence the overall energy budget and reduce the efficiency of 
the production of the observed radiation. We turn now to 
each of this cases. 

Consider, first, the situation in which Y > 1 and the IC 
photons are in the observed range so that some of the ob- 
served radiation may be due to IC rather than synchrotron 
emission. This is an interesting possibility since one might 
expect that the IC process will ease the requirement of 
rather large magnetic fields that is imposed by the syn- 
chrotron process. We show here that, somewhat surpris- 
ingly, this cannot be the case. 

An IC scattering boosts the energy of the photon by 
a factor 7*. Typical IC photons will be observed at the 
energy: 



(hvic)obs 



hq e B 4 

"7e7£ 



12MevBi G (7B100) 



7e 



m e c ■" ■ ' |_ (jn p /m e ) 

(70) 

where Big = B/lGauss and 7^100 = 7b/ioo- The Lorentz 
factor of electrons radiating synchrotron photons which 
are IC scattered on electrons with the same Lorentz factor 
and have energy hv in the observed range is the square 
root of the 7 e required to produce synchrotron radiation 
in the same frequency. The required value for 7 e is rather 



low relative to what one may expect in an external shock 
(in which ^ e ,ext ~ £e(jnp/me)lsh)- In internal shocks we 
expect lower values (j e ,mt ~ £e(^p/w e )) but in this case 
the cquipartition magnetic: field is much stronger (of the 
order of few thousand Gauss, or higher). Thus IC might 
produce the observed photons in internal shocks if e B is 
rather small (of order 10 -5 ). 

These electrons are cooled both by synchrotron and by 
IC. The latter is more efficient and the cooling is enhanced 
by the Compton parameter Y . The cooling time scale is: 



tic 



0~T 



reioo 



)~ 3/4 (/iz//100keVr 1/4 



(71) 

As we see in the following discussion for external shocks, 
tic(W0keV), the IC cooling time if the IC radiation is in 
the observed range (soft gamma-rays) is too long, while for 
internal shocks tjc{W0keV) is marginal. However, even if 
IC does not produce the observed 7-ray photons it still in- 
fluences the process if Y > 1. It will speed up the cooling 
of the emitting regions and shorten the cooling time, t syn 
estimated earlier (Eq. [59]) by a factor of Y. Additionally 
IC also reduces the efficiency by the same factor, and the 
efficiency becomes extremely low as described below. 



8.5. Radiative Efficiency 

The efficiency of a burst depends on three factors: First 
only the electrons' energy is available. This yields a fac- 
tor e e . Second, if e B < e e there is an additional factor of 
min[l, \/e B / e e ] if the IC radiation is not observed. Third, 
there is a specific Lorentz factor, of an electron which 
emits synchrotron (or IC) radiation in the lOOkeV energy 
band. Therefore, only the energy radiated by electrons 
with 7 e > % is observed as soft 7-rays. Assuming a powe r 
law electron distribution with an index p = 2.5 (see [ 103 1 ) 
this gives a factor of ('y e ,min/%) 1 ' 2 (which is valid of course 
provided that ^e.min < 7e)- The total efficiency is the mul- 
tiplication of those three factors and it is given by: 



V /e B/e e ](7e,rmn/7e) 1/2 



(72) 



The efficiency depends first of all on the electrons' en- 
ergy density and to a lesser extend on the magnetic energy 
density. Both should be close to equipartition in order that 
the efficiency will be large. Additionally, in order that 
there will be photons in the lOOkeV range 7e jrra in should 
be smaller than 7. However efficient production of soft 
7-rays requires that "f e ,min will not be too small compared 
with 7. This estimate is, of course, different if the observed 
7-rays are produced by inverse Compton scattering. 

8.6. Internal Shocks 

Internal shocks are the leading mechanism for energy 
conversion and production of the observed gamma-ray ra- 
diation. We discuss, in this section, the energy conversion 
process, the typical radiation frequency and its efficiency. 
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8.6.1. Parameters for Internal Shocks 

Internal shocks take place when an inner shell overtakes 
a slower outer shell. Consider a fast inner shell with a 
Lorentz factor j r that collides with a slower shell whose 
Lorentz factor is j s . If 7 r ^7s ~ 7 then the inner shell will 
overtake the outer one at: 



Rs « 7 2 (5 « 10 14 cm <5 10 7 2 00 



(73) 



where S is the initial separation between the shells in the 
observer's rest frame and Siq = S/10 10 cm and 7100 = 
7/IOO. Clearly internal shocks are relevant only if they 
appear before the external shock that is pr oduce d as the 
shell sweeps up the ISM. We show in section 8.7.1 that the 
necessary condition for internal shocks to occur before the 
external shock is: 

C 3/2 > C ■ (74) 

where £ and C are two dimensionless parameters. The pa- 
rameter, £, characterizes the interaction of the flow with 
the extern al medium and it is defined in Eq. ^2| (see sec- 
The second parameter, £, characterizes the 



tion 



.7.1) 



variability of the flow: 



< 1. 



(75) 



We have seen in section 7.4 that for internal shocks the 
duration of the burst T « A/c and the duration of indi- 
vidual spik es ST w S/c. The observed ratio J\f defined in 
section 2J2 must equal 1/C and this sets £ » 0.01. 



The overall duration of a burst produced by internal 
shocks equals A/c. Thus, whereas external shocks require 
an extremely large value of 7 to produce a very short burst, 
internal shocks can produce a short burst with a modest 
value of the Lorentz factor 7. This eases somewhat the 
baryon purity constraints on source models. The condi- 
tion [74| can be turned into a condition that 7 is sufficiently 
small: 

-i/ 2 T,„- 3 / a /" (76) 



7 < 2800 Co.or 1/z T lc T 3/8 4 /8 , 



where we have used T = A/c and we have defined T\q = 
T/lOs and £0.01 = C/0-01. It follows that internal shocks 
take place in relatively "low" 7 regime. Fig. ^ depicts the 
regimes in the physical parameter space (A, 7) in which 
various shocks are possible. It also depicts an example of 
a T = A/c = Is line. 

Too low a value of the Lorentz factor leads to a large 
optical depth in the internal shocks region. Using Eq. |27| 
for R e , at which the optical depth for Compton scatteri ng 
of the photons on the shell's electrons equals one, Eq. [7^ 
for Rg and the condition R e < Rs we find: 
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(77) 



In addition, the radius of emission should be large 
enough so that the optical depth for 77 — > e + e~ will be less 
than unity (r 77 < 1). There are several ways to consider 
this constraint. The strongest constraint is obtained if one 
demands that the optical depth of an observ ed h i gh e nergy, 
e.g. lOOMeV photon will be less than unifrHpl0|, Fol- 
lowing these calculations and using Eq. |73jto express Rs 
we find: 

7 > 570(Co.oi7i r 1/4 • (78) 



This constraint, which is due to the 77 interaction, is gen- 
erally more important than the constraint due to Compton 
scatte ring : that is r 77 > r e . 

Eq. [76L and the more restrictive Eq. 
a relatively narrow range: 



constrains 7 to 



570 CoToi^io" 174 < IE < 2800 (79) 



This can be translated to a rather narrow range of emission 
radii: 



10 15 cm Clm T w 2 <Rs< 2-5 ■ 10 16 cm T^ 4 l% 4 . (80) 

In Fig. |24|, we plot the allowed regions in the 7 and S 
parameter space. Using the less restrictive r e limit 77 we 
find: 
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< 1E < 2800 Co.oi 72 Til*. 



Three main conclusions emerge from the discussion so 
far. First, if the spectrum of the observed photons extends 
beyond lOOMeV (as was the case in the bursts detected by 
EGRET [Q) and if those high energy photons are emit- 
ted in the same region as the low energy ones then the 
condition on the pair production, r 77 , Eq. [78| is stronger 
than the condition on Compton scattering Eq. ^l|. This 
increases the required Lorentz factors. Second, the Comp- 
ton scattering limit (which is independent of the observed 
high energy tail of the spectrum) poses also a lower limit 
on 7. However, this is usually less restrictive then the 
r 77 limit. Finally, one sees in Fig. that optically thin 
internal shocks are produced only in a narrow region in 
the (S, 7) plane. The region is quite small if the stronger 
pair production limit holds. In this case there is no single 
value of 7 that can produce peaks over the whole range of 
observed durations. The allowed region is larger if we use 
the weaker limits on the opacity. But even with this limit 
there is no single value of 7 that produces peaks with all 
durations. The IS scenario suggests that bursts with very 
narrow peaks should not have very high energy tails and 
that very short bursts may have a softer spectrum. 



8.6.2. Physical Conditions and Emission from Internal 

Shocks 

Provided that the different parts of the shell have com- 
parable Lorentz factors differing by factor of ~ 2, the in- 
ternal shocks are mildly relativistic. The protons' thermal 
Lorentz factor will be of order of unity, and the shocked 
regions will still move highly relativistically towards the 
observer with approximately the initial Lorentz factor 7. 
In front of the shocks the particle density of the shell is 
given by the total number of baryons E/jrripC 2 divided by 
the co- moving volume of the shell at the radius Rs which is 
47ri?^A7. The particle density behind the shock is higher 
by a factor of 7 which is the limiting compression possible 
by Newtonian shocks (assuming an adiabatic index of rela- 
tivistic gas, i.e., 4/3). We estimate the pre-shock density of 
the particles in the shells as: [i?/(7m p c 2 )]/(47r(57 2 ) 2 7A). 
We introduce 7^ as the Lorentz factor of the internal 
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shock. As this shock is relativistic (but not extremely rel- 
ativistic) 7 int is of order of a few. Using Eq. [l8] for the 
particle density n and the thermal energy density e behind 
the shocks we hnd: 
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eint = (Jint I TUntrripC . 
We have defined here A 12 = A/10 12 cm. Using Eq. 
find: 

l/2„l/2 



(81) 



(82) 
ilwe 



Bi n t = 6-10 5 Gauss ej E< 



52 7ioo- 3 A 12 - 1 / 2 ^ 10 - 1 ( 7mt /2) 1 / 2 . 

(83) 

Using Eqs. ||, || || and fell we can estimate the typi- 
cal synchrotron frequency from an internal shock. This is 
the synchrotron frequency of an electron with a "typical" 
Lorentz factor: 

(u \ i hqeB 2 

(nv S yn)obs\( 7c ) = 7e7 
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(84) 
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The corresponding observed synchrotron cooling time is: 



tsyn\(~f e ) = 1-3 • 10" 



sec e^S 2 ^ 



Ai2ll m E b2 -\ lmt /2)-\ 
(85) 

Using Eq. |33j we can express 7 e ,mm in terms of 7;„t to 
estimate the minimal synchrotron frequency: 
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The energy emitted by a "typical electron" is around 
220keV. The energy emitted by a "minimal energy" elec- 
tron is about one order of magnitude lower than the typical 
observed energy of ~ lOOkeV. This should correspond to 
the break energy of the spectrum. This result seems in a 
good agreement with the observations. But this estimate 
might be misleading as both e# and e e might be signif- 
icantly lower than unity. Still these values of (hv syn ) t, s 
are remarkably close to the observations. One might hope 
that this might explain the observed lower cutoff of the 
GRB spectrum. Note that a lower value of es or e e might 
be compensated by a higher value of jmt- This is advan- 
tageous asshocks with higher 7;„ t are more efficient (see 
section 8.6.4). 

The synchrotron cooling time at a given frequency (in 
the observer's frame) is given by: 



tsyn(hv) = 2 X 10" 



sec e 



-3/4 
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(87) 
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We recover the general trend t syn cx {hv)~ x l 2 of syn- 
chrotron emission. However if (as we expect quite gen- 
erally) this cooling time is much shorter than T ang it does 
not determine the width of the observed peaks. It will 
correspond to the observed time scales if, for example, cb 
is small. But then the "typical" photon energy will be far 
below the observed range. Therefore, it is not clear this 
relation can explain the observed dependence of the width 
of the bursts on the observed energy. 



.6.3. Inverse Compton in Internal Shocks 



The calculations of section 3.4 suggest that the typ- 
ical Inverse Compton (IC) (actually synchrotron - self 
Compton) radiation from internal shocks will be at en- 
ergy higher by a factor 7 2 then the typical synchrotron 
frequency. Since synchrotron emission is in the keV range 
and 7e,mm ~ m p /m e , the expected IC emission should 
be in the GeV or even TeV range. This radiation might 
contribute to the prompt very high energy emission that 
accompanies some of the GRBs [B3| . 

However, if the magnetic field is extremely low: eb ~ 
10~ 12 then we would expect the IC photons to be in the 
observed ~ lOOkeV region: 



hvw-int = 800 keV {e B /l^ 12 ) 1/2 
x a- 1 A- 1/2 7io 2 [ le /(mp/me)} 4 EU 2 [ 7int /2] ^ . 



Using Eqs. |7lj and |S3| we find that the cooling time for 
synchrotron-self Compton in this case is: 

T ICi nt = 1 sec e^/^es/lO- 12 )- 1 / 2 (89) 
x5 10 2 A 12 7 1 5 00 [7 e /(™p/me)]- 1 i; 52 1 (7 m t/2)- 1 . 

This is marginal. It is too large for some bursts and possi- 
bly adequate for others. It could possibly be adjusted by 
a proper choice of the parameters. It is more likely that 
if Inverse Compton is important then it contributes to the 
very high (GeV or even TeV) signal that accompanies the 
lower energy GRB (see also [ [250f ). 



8.6.4. Efficiency in Internal Shocks 

The elem entary unit in the internal shock model (see 
section 



7.4) is a a binary (two shells) encounter between 



a rapid shell (denoted by the subscript r) that catches up 
a slower one (denoted s). The two shells merge to form 
a single shell (denoted m). The system behaves like an 
inelastic collision between two masses m r and m s . 

The efficiency of a single c ollisio n between two shells was 



calculated earlier in section 8.1.1. For multiple collisions 



the efficiency depends on the nature of the random distri- 
bution. It is highest if the energy is distributed equally 
among the different shells. This can be explained analyti- 
cally. Consider a situation in which the mass of the shell, 
rrii is correlated with the (random) Lorentz factor, 7$ as 
m,i oc 7?. Let all the shells collide and merge and only 
then emit the thermal energy as radiation. Using con- 
servation of energy and momentum we can calculate the 
overall efficiency: 



e = 1 - £7? /j^r 1 ^ 1 



(90) 



Averaging over the random variables ji, and assuming a 
large number of shells N — > 00 we obtain: 
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This formula explains qualitatively the numerical results: 
the efficiency is maximal when the energy is distributed 
equally among the different shells (which corresponds to 
V = -l). 

In a realistic situation we expect that the internal en- 
ergy will be emitted after each collision, and not after all 
the shells have merged. In this case there is no simple 
analytical formula. However, numerical calculations show 
that the efficiency of this process is low (less than 2%) if 
the initial spread in 7 is only a factor of two |Q . However 
the efficiency could be much higher p3| . The most effi- 
cient case is when the shells have a comparable energy but 
very different Lorentz factors. In this case (77 = —I, and 
spread of Lorentz factor Jmax/lmin > 10 3 ) the efficiency 
is as high as 40%. For a moderate spread of Lorentz factor 
Imax/lmin = 10, with j] = -1, the efficiency is 20%. 

The efficiency discussed so far is the efficiency of con- 
version of kinetic energy to internal energy. One sho uld 
multiply this by the radiative efficiency, discussed in 
(Eq. [72|) to obtain the overall efficiency of the process. 
The resulting values may be rather small and this indi- 
cates that some sort of beaming may be required in most 
GRB models in order not to come up with an unreasonable 
energy requirement. 



8.6.5. Summary - Internal Shocks 

Internal shocks provide the best way to explain the ob- 
served temporal structure in GRBs. These shocks, that 
take place at distances of ~ 10 15 cm from the center, con- 
vert two to twenty percent of the kinetic energy of the flow 
to thermal energy. Under reasonable conditions the typ- 
ical synchrotron frequency of the relativistic electrons in 
the internal shocks is around lOOkeV, more or less in the 
observed region. 

Internal shocks require a variable flow. The situation in 
whi ch an inner shell is faster than an outer shell is unsta- 
ble ]25l| . The instability develops before the shocks form 
and it may affect the energy conversion process. The full 
implications of this instability are not understood yet. 

Internal shocks can extract at most half of the shell's 
energy j3^, [33[ . Highly relativistic flow with a kinetic 
energy and a Lorentz factor comparable to the original 
one remains after the internal shocks. Sari & Piran | p0| 
pointed out that if the shell is surrounded by ISM and col- 
lisionless shock occurs the relativistic shell will dissipate 
by "external shocks" as well. This predicts an additional 
smooth burst, with a comparable or possibly greater en- 
ergy. This is most probably the source of the observed 
"afterglow" seen in some counterparts to GRBs which we 
discuss later. This leads to the Internal-External scenario 



1 252, B(J, Eq] in which the GRB itself is produced by an 
Internal shock, while the "afterglow" that was observed to 
follows some GRBs is produced by an external shock. 

The main concern with the internal shock model is its 
low efficiency of conversion of kinetic energy to 7-rays. 
This could be of order twenty percent under favorable 
conditions and significantly lower otherwise. If we assume 
that the "inner engine" is powered by a gra vitat ional bind- 
ing energy of a compact object (see section 10. 1) a low effi- 
ciency may require beaming to overcome an overall energy 
crisis. 



8.7. Shocks with the ISM - External shocks 

We turn now to the interaction of a relativistic shell with 
the ISM. We have seen in section 7.4 that external shocks 



cannot produce bursts with a complicated temporal struc- 
ture. Still it is worthwhile to explore this situation. First, 
there are some smooth bursts that might be produced in 
this way. Second, one needs to understand the evolution 
of external shocks in order to see why they cannot satisfy 
the condition Re/ J 2 < A. Third, it is possible that in 
some bursts emi ssion is observed from both internal and 
external shocks |253| ] . Finally, as we see in the following 
section ^| the observed afterglow is most likely produced 
by external shocks. 



8.7.1. Newtonian vs. Relativistic Reverse Shocks 

The interaction between a relativistic flow and an exter- 
nal medium depends, like in SNRs, on the Sedov length, 
I = (E /riismmpC 2 ) 1 / 3 . The ISM rest mass energy within a 
volume I 3 equals the energy of the GRB: E. For a canon- 
ical cosmological burst with E » 10 52 ergs and a typical 



ISM density ni 



lparticle/cnr we have / ss 10 cm. A 



second length scale that appears in the problem is A, the 
width of the relativistic shell in the observer's rest frame. 



There are two possible types of external shocks [ 233 1 . 
They are characterized according to the nature of the re- 
verse shock: Newtonian Reverse Shock (NRS) vs. Rela- 
tivistic Reverse Shock (RRS). If the reverse shock is rela- 
tivistic (RRS) then it reduces significantly the kinetic en- 
ergy of each layer that it crosses. Each layer within the 
shell loses its energy independently from the rest of the 
shock. The energy conversion process is over once the re- 
verse shock crosses the shell (see Fig. |l3|). A Newtonian 
or even mildly relativistic reverse shock (NRS) is com- 
paratively weak. Such a shock reduces the energy of the 
layer that it crosses by a relatively small amount. Signifi- 
cant energy conversion takes place only after the shock has 
crossed the shell several time after it has been reflected as 
a rarefraction wave from the inner edge (see Fig. |l^) . The 
shell behaves practically like a single object and it loses its 
energy only by the time that it accumulates an external 
mass equal to M/7. 

The question which scenario is taking place depends on 
the parameters of the shell relative to the parameters of 
the ISM. As we see shortly it depends on a single dimc n- 
sionless parameter £ constructed from I, A and 7: |23 



(92) 



As the shell propagates outwards it is initially very 
dense and the density ratio between the shell and the ISM, 
/ = 714/711, is extremely large (more specifically / > 7 2 ). 
The reverse shock is initially Newtonian (see Eq. B3|) . Such 
a shock converts only a small fraction of the kinetic energy 
to thermal energy. As the shell propagates the density 
ratio, /, decreases (like R~ 2 if the width of the shell is 
constant and like R~ 3 if the shell is spreading). Eventu- 
ally the reverse shock becomes relativistic at Rn where 
/ = 7 2 . The question where is the kinetic energy con- 
verted depends on whether the reverse shock reaches the 
inner edge of the shell before or after it becomes relativis- 
tic. 
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There are four different radii that should be consid- 
ered. The following estimates assume a spherically sym- 
metric shell, or that E and M are energy and rest mass 
divided by the fraction of a sphere into which they are 
launched. The reverse shock becomes relativistic at Rn, 
where / = 714/ni = 1: 



R 



N 



z 3/2 /A l/2 7 2 



(93) 



Using the expression for the velocity of the reverse shock 
into the shell (Eq. [l6]) we find that the reverse shock 
reaches the inner edge of the shell at Ra ]233| : 



R A = I 3 /* A 1 / 4 



(94) 



A third radius is i? 7 , where the shell collects an ISM mass 
of M/7 [^7], ^8|. For NRS this is where an effective energy 
release occurs: 
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where we defined ri\ = n,;. sr „/l particle/cm . Finally we 
have Rs = S^f 2 , (see Eq. |7^). The different radii are 
related by the dimensionless parameter £, and this deter- 
mines the character of the shock: 

Rs/c = R A /e /2 = R~,e = R N /e (96) 

If £ > 1 then: 

Rs < Ra < R-y < Rn- (97) 

The reverse shock remains Newtonian or at best mildly rel- 
ativistic during the whole energy extraction process. The 
reverse shock reaches the inner edge of the shock at Ra 
while it is still Newtonian. At this stage a reflected rar- 
efraction wave begins to move forwards. This wave is, in 
turn, reflected from the contact discontinuity, between the 
shell's material and the ISM material, and another reverse 
shock begins. The overall outcome of these waves is that 
in this case the shell acts as a single fluid element of mass 
M « E/jc 2 that is interacting collectively with the ISM. 
It follows from Eq. |3^ that an external mass m = M/"f 
is required to reduce 7 to 7/2 and to convert half of the 
kinetic energy to thermal energy. Energy conversion takes 
place at i? 7 Comparison of R 1 with R e (equation ^) shows 
that the optical depth is much smaller than unity. 

If the shell propagates with a constant width then 
Rn/€ = R 1 = V£Ra (see Fig. ||) and for £ > 1 the 
reverse shock remains Newtonian during the energy ex- 
traction period. If there are significant variations in the 
particles velocity within the shell it will spread during the 
expansion. If the typical variation in 7 is of the same or- 
der as 7 then the shell width increases like -R/7 2 . Thus A 
changes with time in such a manner that at each moment 
the current width, A(i), satisfies A(t) ~ max[A(0), fi/7 2 ]. 
This delays the time that the reverse shock reaches the 
inner edge of the shell and increases Ra- It also reduces 
the shell's density which, in turn, reduces / and leads to a 
decrease in Rn. The overall result is a triple coincidence 



and a significant energy conversion in the reverse shock 
as well. This means that due to spreading a shell which 
begins with a value of £ > 1 adjusts itself so as to satisfy 
£ = !• 

For £ > 1 we find that T radia i ~ T ang ~ R y /j 2 > A. 
Therefore, NRS can produce only smooth bursts. The 
bursts' duration is determined by the slowing down time 
of the shell. In section [7] we have shown that only one time 
scale is possible in this case. Given the typical radius of 
energy conversion, i? 7 this time scale is: 

ST « T obs « i? 7 /( 7 Ic) « i? 7 /(7 2 c) « 170sec J B 5 1 2 /3 nr 1/ Voo /3 ! 

(98) 

If 7 or A are larger then £ < 1 . In this case the order is 
reversed: 

R N < R-y < Ra- (99) 

The reverse shock becomes relativistic very early (see Fig. 
25). Since j s h = 72 *C 7 the relativistic reverse shock 
converts very efficiently the kinetic energy of the shell to 
thermal energy. Each layer of the shell that is shocked 
loses effectively all its kinetic energy at once and the time 
scale of converting the shell's kinetic energy to thermal 
energy is the shell crossing time. The kinetic energy is 
consumed at Ra, where the reverse shock reaches the in- 
ner edge of the shell. Using Eq. |94| for Ra and Eq. |^] we 
find that at Ra 
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Note that 7^ is independent of 7. The observed radial or 
angular time scales are: 



T, 



radial 



Tang « #a/tI « A/c = 30 sec A 12 . (101) 



Thus even for RRS we find that 6T ~ T and there is only 
one time scale. This time scale depends only on A and it is 
independent of 7! Spreading does not affect this estimate 
since for £ < 1 spreading does not occur before the energy 
extraction. 

In the following discussions we focus on the RRS case 
and we express all results in terms of the parameter £. 
By setting £ < 1 in the expressions we obtain results cor- 
responding to RRS, and by choosing £ = 1 in the same 
expressions we obtain the spreading NRS limit. We shall 
not discuss the case of non-spreading NRS (£ 3> 1), since 
spreading will always bring these shells to the mildly rel- 
ativistic limit (£ ~ 1). Therefore, in this way, the same 
formulae are valid for both the RRS and NRS limits. 

if £ 
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that internal shocks 
cs. If £ < 1 then the 
1 
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Rj pa Ra with a mildly relativistic reverse shock 



> 1 it follows from Eq 
will take place before external shoe 
condition for internal shocks Rs < Ra becom es Eq . 
£3/2 ^ £ As we have seen earlier (see section 8.6.1) this 
sets an upper limit on 7 for internal shocks. 



8.7.2. Physical Conditions in External Shocks 

The interaction between the outward moving shell and 
the ISM takes place in the form of two shocks: a forward 
shock that propagates into the ISM and a reverse shock 
that propagates into the relativistic shell. This results in 
four distinct regions: the ISM at rest (denoted by the sub- 
script 1 when we consider properties in this region), the 
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shocked ISM material which has passed through the for- 
ward shock (subscript 2 or f), the shocked shell material 
which has passed through the reverse shock (3 or r), and 
the unshocked material in the shell (4). See Fig. gl]. The 
nature of the emitted radiation and the efficiency of the 
cooling processes depend on the conditions in the shocked 
regions 2 and 3. Both regions have the same energy density 
e. The particle densities 712 and 713 are, however, differ- 
ent and hence the effective "temperatures," i.e. the mean 
Lorentz factors of the random motions of the shocked pro- 
tons and electrons, are different. 

The bulk of the kinetic energy of the shell is converted 
to thermal energy via the two shocks at around the time 
the shell has expanded to the radius Ra- At this radius, 
the conditions at the forward shock are as follows, 

72 = 7£ 3/4 > n 2 = 47 2 ni, e 2 = A-^lmrripC 2 , (102) 
while at the reverse shock we have 



and 
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Substitution of 7^ = 72 = 7^ 3 ^ 4 in Eq. yields: 
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If the magnetic field in region 2 behind the forward shock 
is obtained purely by shock compression of the ISM field, 
the field would be very weak, with cb <C 1. Such low 
fields are incompatible with observations of GRBs. We 
therefore consider the possibility that there may be some 
kind of a turbulent instability which may bring the mag- 
netic field to approximate equipartition. In the case of 
the reverse shock, magnetic fields of considerable strength 
might be present in the pre-shock shell material if the orig- 
inal exploding fireball was magnetic. The exact nature of 
magnetic field evolution during fireba ll exp ansion depends 
on several assumptions. Thompson 224 found that the 
magnetic field will remain in equipartition if it starte d off 
originally in equipartition. Meszaros, Laguna & Rees [ 240 
on the other hand estimated that if the magnetic field was 
initially in equipartition then it would be below equipar- 
tition by a factor of 10 -5 by the time the shell expands 
to i?A- It is uncertain which, if either, is right. As in the 
forward shock, an instability could boost the field back to 
equipartition. Thus, while both shocks may have < 1 
with pure flux freezing, both could achieve eb — > 1 in the 
presence of instabilities. In principle, eb could be different 
for the two shocks, but we limit ourselves to the same cb 
in both shocks. 

In both regions 2 and 3 the electrons have a power law 
distribution with a minimal Lorentz factor J em in given 
by Eq. [33] with the corresponding Lorentz factors for the 
forward and the reverse shock. 

8.7.3. Synchrotron Cooling in External Shocks 

The typical energy of synchrotron photons as well as the 
synchrotron cooling time depend on the Lorentz factor j e 
of the relativistic electrons under consideration and on the 
strength of the magnetic field. Using Eq. for 7 e ,mm we 
find the characteristic synchrotron energy for the forward 
shock: 



[JlV 'syn) obs |- 



160 keV es /2 ee(72/100) 4 nj /2 , (105) 



= 0.085 sec es 1 e f T 1 (72/100)- 4 nJ- 



(106) 



The characteristic frequency and the corresponding cool- 
ing time for the "typical" electron are larger by a factor of 
[(p-2)/(p- 1)] 2 and shorter by a factor of [(p-2)/(p- 1)] 2 , 
correspondingly. 

These photons seems to be right in the observed soft 
gamma-ray range. Howev er, o ne should recall that the fre- 



quency calculated in Eq. 105 depends on the forth power 
of 72. An increase of the canonical 72 by a factor of 3 (that 
is 72 = 300 instead of 72 = 100) will yield a "typical" syn- 
chrotron emission at the 16MeV instead of 160keV. The 
Lorentz factor of a "typical electron" in the reverse shock 
is lower by a factor £ 3 / 2 . Therefore the observed energy is 
lower by a factor £ 3 while the cooling time scale is longer 
by a factor £~ 3 / 4 . 

Alternatively we can check the conditions in order that 
there are electrons with a Lorentz factor j e that be emit- 
ting soft gamma-rays with energies ~ lOOkeV. Using Eq. 
56 we calculate 7 e : 
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Electrons with 7 e = % are available in the shocked mate- 
rial if 7e,min < 7e- This corresponds to the condition 
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in the reverse shock, and the condition 
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e e 1 < 0.8 e"; 
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(109) 



in the forward shock. Since by definition e e < 1, we see 
that the reverse shock always has electrons with the right 
Lorentz factors to produce soft gamma-ray synchrotron 
photons. However, the situation is marginal in the case of 
the forward shock. If 7 > 100 and if the heating of the 
electrons is efficient, i.e. if e e \ s ~ 1, then most of the elec- 
trons may be too energetic. Of course, as an electron cools, 
it radiates at progressively softer energies. Therefore, even 
if 7min is initially too large for the synchrotron radiation 
to be in soft gamma-rays, the same electrons would at a 
later time have 7 e ~ "f e and become visible. However, the 
energy remaining in the electrons at the later time will also 
be lower (by a factor 7/7min), which means that the burst 
will be inefficient. For simplicity, we ignore t his ra diation. 

Substituting the value of "f e from equation 107 into the 
cooling rate Eq. 58 we obtain the cooling time scale as a 
function of the observed photon energy to be 
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Eq. 110 is valid for both the forward and reverse shock, 



and is moreover independent of whether the reverse shock 
is relativistic or Newtonian. 

The cooling time calculated above sets a lower limit to 
the variability time scale of a GRB since the burst cannot 
possibly contain spikes that are shorter than its cooling 
time. However, it is unlikely that this cooling time actu- 
ally determines the observed time scales. 



8.8. The Internal - External Scenario 

Internal shocks can convert only a fraction of the total 
energy to radiation (3|, |69(| . After the flow has pro- 
duced a GRB via internal shocks it will interact via an 
external shock with the surrounding medium pb]| . This 
will produce the afterglow - a signal that will follow the 
GRB. The idea of an afterglow in other wavelengths was 
suggested earlier jlT], [l8|, ^lj but it was suggested as a fol- 
low up of the, then standard, external shock scenario. In 
this case the afterglow would have been a direct continu- 
ation of the GRB activity and its properties would have 
scaled directly to the properties of the GRB. 

According to internal-external models (internal shocks 
for the GRB and external shocks for the afterglow) dif- 
ferent mechanisms produce the GRB and the afterglow. 
Therefore the afterglow should not be scaled directly to 
the properties of the GRB. This was in fact seen in the 
recent afterglow observations |2^, M. In all models of 
external shocks the observed time satisfy t cx R/j% and 
the typical frequency satisfies v cx 7* . Since most of the 
emission takes place at practically the same radius and all 
that we see is the variation of the Lorentz factor we ex- 
pect quite generally p5[ : v cx t 2±l . The small parameter 
l reflects the variation of the radius and it depends on the 
specific assumptions made in the model. We would expect 



that trc/t-f 



5 and topt/t^ 



300. 



The observations of 
GRB970508 show that (topt/t^observed ~ 10 4 . This is in 
a clear disagreement with the single external shock model 
for both the GRB and the afterglow. 

Under quite general conditions the initial typical syn- 
chrotron energy for either the forward or the reverse ex- 
ternal shock may fall in the soft GRB band. In this case 
the initial stage of the afterglow migh t overlap the 7-ray 
emission from the internal shock |£53[. The result will be 
superposition of a rapidly varying signal on top of a long 
smooth and softening pulse. This possibility should be 
explored in greater detail. 



9. AFTERGLOW 

It is generally believed that the observed afterglow re- 
sults from slowing down of a relativistic shell on the ex- 
ternal ISM. The afterglow is produced, in this case, by 
an external shock. A second alternative is of "continu- 
ous emission". The "inner engine" that powers the GRB 
continues to emit energy for much longer duration with 
a lower amplitude [ p6[ and may produce the earlier part 
(first day or two in GRB970228 and GRB970508) of the 
afterglow. It is most likely that both processes take place 
to some extent p6|. We discuss in this section theoretical 
models for the production of the afterglow focusing on the 
external shock model. 
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9.1. Hydrodynamics of a Slowing Down Relativistic Shell 

Within the external shock model there are several pos- 
sible physical assumptions that one can make. The "stan- 
dard" model assumes adiabatic hydrodynamics (energy 
losses are negligible and do not influence the hydrodynam- 
ics), slow cooling (the electrons radiate a small fraction of 
the energy that is generated by th e shock) and synchrotron 

. However there are 



emission [l|, |2_ 
other possibilities. First, the electrons' energy might be 
radiated rapidly. In this case the radiation process is fast 
and the observed flux is determined by the rate of energy 
generation by the shock. If the electrons carry a signifi- 
cant fraction of the total internal energy fast cooling will 
influence the hydrodynamics which will not be adiabatic 
any more. In this case we have a radiative solution |2^, ^4) 
which differs in its basic scaling laws from the adiabatic 
one . The different possibilities are summarized in Table 
9~l| 



9.1.1. A Simple Collisional Model 

We consider first a simple model for the slowing down of 
the shell. In this model the slowing down is described by a 
series of infinitesimal inelastic collisions between the shell 
and infinitesimal external masses. We assume a homoge- 
neous shell described by its rest frame energy M (rest mass 
and thermal energy) and its Lorentz factor 7. Initially, 
Eq = Moc 2 7o- The shell collides with the surrounding 
matter. We denote the mass of the ISM that has already 
collided with the shell by m(R). As the shell propagates 
it sweeps up more ISM mass. Additional ISM mass ele- 
ments, dm, which are at rest collides inelastically with the 
shell. 

Energy and momentum conservation yield: 



d;~) 



and 



dm 



= (7 - l)dm, 



(111) 



(112) 



where dE is the thermal energy produced in this collision. 
We define e as the fraction of the shock generated thermal 
energy (relative to the observer frame) that is radiated. 
The incremental total mass satisfies: 

dM = (1 - e)dE + dm = [(1 - e) 7 + e]dm. (113) 

These equations yields analytic relations between the 
Lorentz factor and the total mass of the shell: 



(7-IX7+I) 



l-2e 



( 7 0-l)(70 + I) 1 " 26 

and between m{R) (and therefore R) and 7 
m(R) 



(M/M )- 2 , (114) 



M 

1 



= -(7o-l) 1/2 (7o + l) 1/2 



(115) 



x / ( 7 '-i)- 3 / 2 ( 7 ' + i)-3/2+^y. 

-'70 

These relations completely describe the hydrodynamical 
evolution of the shell. 
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Adiabatic Hydrodynamics 


Radiative Hydrodynamics 


Slow Cooling 


Arbitrary e e 


Impossible 


Fast Cooling 


e e < 1 


e e « 1 



Table 6: Afterglow Models 



Two basic features can be seen directly from Eq. |116| . 
First, we can estimate the ISM mass m that shoul d be 
swept to get significant deceleration. Solving Eq. |116] 
with an upper limit 7o/2 and using 7q > 1 we obtain the 
well known result: a mass m = Mo/ (2jo) is required to 
reach 7 = 70/2. Apparently this result is independent of 
the cooling parameter e. 

A second simple result can be obtained in the limit that 
70 >7> 1: 



m(R) = 



(2 - 5)70 \7o 



7 



-2+e 



(116) 



so that 7 oc R 3 ^ 2 e \ For e = this yields the well known 
adiabatic result: 



47T q 22 

—R nrripC 7 = E , 



(117) 



and 7 oc R~ 3/2 |||, |§ || || |247). For e = 1 this yields 
the completely radiative result: 



47T 



-R 3 nm p c 2 jj = E , 



(118) 



and 7 oc R- 3 g3§, ||, ||. 

For comparison with observations we have to calculate 
the observed time that corresponds to different radii and 
Lorentz factors. The well known formula 



obs 



R 
2^0 



(119) 



is valid only for emission along the line of sight from 
a sh ell that propagates with a constant velocity. Sari 
253] pointed out that as the shell decelerates this formula 
should be used only in a differential sense: 



dt b s — 



dR 

27V 



(120) 



Eq. |120] should be combined with the relation 117 or hi 



and integrated to get the actual relation between observed 
time and emission radius. For an adiab atic expansi on, f or 

[p53|. Eq. is 



example, this yields: 



R/16Y 



valid only along the line of sight. The situation is compli- 
cated further if we recall that the emission reaches the ob- 
serve from an angle of order 7 -1 around the line of sight. 
Averagin g on all angles yields another numerical factor 
1 254, 255, 25E] and altogether we get 



t. 



R 



oh* 



C 7 72C 



(121) 



where the value of the numerical factor, c ga , depends on 
the details of the solution and it varies between ~ 3 and 



~ 7. Using Eqs. 121 and 117 or 118 we obtain the follow- 
ing relations between R, 7 and t: 



pu\ ~ / {^Et I -Km p nc) 
(> -\ {Act/ L) 1 / 1 L, 



1/4 



7(*) = 



nnm v c 



f (3.B/256 
\ (Act/L)- 3 /\ 



,5 f 3\l/8 



ad, 
rad, 

ad, 
rad, 



(122) 



(123) 



where L = {3E / AimnipC 2 ^) 1 / 3 is the radius where the ex- 
ternal mass equals the mass of the shell. 

One can procee d an d us e th e relation between R and 
7 and t bs (Eqs. |l22| and 123) to estimate the physical 
conditions at the shocked material using Eqs. |4J. Then 
one can estimate the emitted radiation from this shock us- 
ing Eqs. |56| and |57j. However, before doin g so we explore 
the Blandford-McKee self similar solution |238| , which de- 
scribes more precisely the adiabatic expansion. This solu- 
tion is inhomogeneous with a well determined radial pro- 
file. The matter at the front of the shell moves faster than 
the average speed. This influences the estimates of the 
radiation emitted from the shell. 



9.1.2. The Blandford-McKee Self-Similar Solution 



Blandford & McKee [238] discovered a self-similar solu- 
tion that describes the adiabatic slowing down of an ex- 
tremely relativistic shell propagating into the ISM. Using 
several simplifications and some algebraic m anip ulations 
we rewrite the Blandford-McKee solution as [253]: 



n(r,t) = 4n 7 (t) [l + 16 7 (i) 2 (l - r/R)] 5/4 , 

7 (r,t) = 7 (*)[l + 167(*) 2 (l-»-/i2)]~ 1/a , 

e(r, t) = 4nm p c 2 ~/(t) 2 [l + 16 7 (<) 2 (1 - r/R)] ~ 17 (l24) 

where n(r, t), e(r,t) and j(r, t) are, respectively, the den- 
sity, energy density and Lorentz factor of the material be- 
hind the shock (not to be confused with the ISM density 
n) and j(t) = j(R(t)) is the Lorentz factor of material just 
behind the shock. n{r,t) and e(r, t) are measured in the 
fluid's rest frame while 7(7*, t) is relative to an observer 
at rest. The total energy in this adiabatic flow equals 
E = Eq, the initial energy. The scaling laws of R(t) and 
7(i) that follow from these profiles and from the condition 
that the total energy in the flow equals E is: 



R(t) 
lit) 



17 Et 



1/4 



= 3.2 x 10 16 cm EU 4 n~ 1/4 tl/\ 



17 E 



4 V Tmm p c 5 t 3 



1/8 



260£ 5 1 2 / V 1/8 t; 3/8 .(125) 



122 



The scalings 125 are consistent with the scalings 
123 which were derived using conservation of energy 
momentum. They provide the exact numerical factor that 



and 
and 
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cannot be calculated by the simple analysis of section |9.1.l| . 
These equations can serve as a starting point for a detailed 
radiation emission calculation and a comparison with ob- 
servations. 

The Blandford-McKee solution is adiabatic and as such 
it does not allow for any energy losses. With some sim- 
plifying assumptions it is possible to derive a self-similar 
radiative solution in which an arbitrary fractio n of the en- 
ergy generated by the shock is radiated away 257]. 



9.2. Phases in a Relativistic Decelerating Shell 

There are several phases in the deceleration of a rela- 
tivistic shell: fast cooling (with either radiative or adi- 
abatic hydrodynamics) is followed by slow cooling (with 
adiabatic hydrodynamics). Then if the shell is non spheri- 
cal its evolution changes and a phase of sideways expansion 
and much faste r slo w down begins when the Lorentz fac- 
tor reaches [258]. Finally the shell becomes Newtonian 
when enough mass is collected and 7 ~ 1 . In the following 
we estimate the time scale for the different transitions. We 



define j e 



c^e e (m p /m e )j and t obs = (1 + z)R/Ac t cj 2 



such that the factors c 7 and Ct reflect some of the uncer- 
tainties in the model. The canonical values of these factors 
are: c 7 « 0.5 and Ct ~ 1. 

The deceleration begins in a fast cooling phase. If e e is 
close to unity than this cooling phase will also be radiative. 
The first transition is from fast to slow cooling. There are 
several different ways to estimate this transition. One can 
compare the cooling time scale to the hydrodynamic time 
scale; alternatively one can calculate the fast cooling rate 
(given by the rate of energy generation by the shell) and 
compare it to the slow cooling rate (given by the emissiv- 
ity of the relativistic electrons). We have chosen here to 
calculate this time as the time when the "typical electron" 
cools - that is when v r — v m : 



tf a 



210 days e%eiE 52 ni, 

4.6 days e 7 J 5 el /5 E^ 5 j~ 4 /5 nl /5 



ad, 
rad. 



(126) 



All methods of estimating tt s give the same dependence 
on the parameters. However, the numerical factor is quite 
sensitive to the definition of this transition. 

If the solution is initially radiative the transition from 
fast to slow cooling and from a radiative hydrodynamics 
to adiabatic hydrodynamics takes place at: 

t rad ^ ad = 1.3 day Sj E 5 4 2 /5 n? /5 4/ 5 4 /5 ((l + z)/2) 12 / 5 (127) 

(c 7 /0.5) 14 / 5 cr 12/5 ( 7 o/100)- 4 / 5 . 

During a radiative evolution the energy in the sho ck d e- 
creases with time. The energy that appears in Eqs. 122 in 
the radiative scalings is the initial energy. When a radia- 
tive shock switches to adiabatic evolution, it is necessary 
to use the reduced energy to calculate the subsequent adi- 
abatic evolution. The energy £7,52 which one should use 
in the adiabatic regime is related to the initial -£^52 of the 
fireball by 



E f t52 = 0.022e 



-3/5 



B 



-3/5 £4/5 



,527l00 



4/5 -2/5 



(128) 



If the shell is not spherical and it has an ope ning angle: 
8, then the evolution will change when 7 ~ B~ x [258|. Ear- 
lier on the jet expands too rapidly to expand sideways and 



it evolves as if it is a part of a spherical shell. After this 
stage the jet expands sideways and it accumulates much 
more mass and slows down much faster. This transition 
will take place, quite generally, during the adiabatic phase 
at: 

t g sa 0.5daysS 5 1 2 /3 nr 1/3 ((l + z)/2)(6»/0.1) 8/3 cr 1 (129) 

The shell eventually becomes non relativistic. This hap- 
pens at: R » I = {AEq/ AirnimpC 2 ) 1 ^ for an adiabatic 
solution. This corresponds to a transition at: 



tNR.ad »J/c» 300 days E^ n x 



1/3 -1/3 



(130) 



A radiative shell loses energy faster and it becomes non 

relativistic at R = L = l/^ 3 = (4i? /47rn ism m p c 2 7o) 1/ ' 3 . 
This will take place at: 

t N R,rad~M days £; 5 1 2 /3 nr 1/3 ( 7 o/100)- 1 / 3 . (131) 

However, the earlier estimate of the transition from fast 
to slow cooling suggests that the shell cannot remain ra- 
diative for such a long time. 



9.3. Synchrotron Emission from a Relativistic 
Decelerating Shell 

We proceed now to estimate the expected instantaneous 
spectrum and light curve from a relativistic decelerating 
shell. The task is fairly simple at this stage as all the 
ground rules have been set in the previous sections. We 
limit the discussion here to a spherical shock propagat- 
ing into a homogeneous external matter. We consider two 
extreme limits for the hydrodynamic evolution: fully ra- 
diative and fully adiabatic. If e e is somewhat less than 
unity during the fast cooling phase (t < tf s ) then only a 
fraction of the shock energy is lost to radiation. The scal- 
ings will be intermediate between the two limits of fully 
radiative and fully adiabatic discussed here. 

For simplicity we assume that all the observed radiation 
reaches the observer from the front of the shell and along 
the line of sight. Actually to obtain the observed spectrum 
we should integrate over the shell's profile and over differ- 
ent angles relative to the line of sight. A full calculation 
|259j ] of the integrated spectrum over a Blandford-McKee 
profile shows that that different radial points from which 
the radiation reaches the observers simultaneously con- 
spire to have practically the same synchrotron frequency 
and therefore they emit the same spectrum. Hence the ra- 
dial integration over a Blandford-McKee profile does not 
change the observed spectrum (note that this re sult differs 
from the calculation for a homogeneous shell ]255[ ). On 
the other hand the contribution from angles away from 
the line of sight is important and it shapes the observed 
spectrum, the light curve and the shape of the Afterglow 
(see Fig|§). 

The instantaneous synchrotron spe ctra f rom a relativis- 
tic shock were described in section 3.2.4 . They do not 
depend on the hydrodynamic evolution but rather on the 
instantaneous conditions at the shock front, which deter- 
mines the break energies v c and v m . The only assumption 
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made is that the shock properties are fairly constant over 
a time scale comparable to the observation time t. 

Using the adiabatic shell conditions (Eqs. 12^- 123) , Eqs. 
44 for the shock conditions, Eq. |56| for the synchrotron en- 
ergy and Eq. p3 for the "cooling energy" we find: 



I'm 



= 2.7 x 10 12 Hz e 



3/2 £T 1/2 " 



B 



"1/2 



52 



= 5.7x10 



14 



tt l/2 2pl/2,- 

Hz e B e e E 52 t d 



3/2 



F v 



= 1.1 x 10° /xJ € 1 J 2 E 52 n l ''D 



1/2; 



(132) 



where £<j is the time in days, D28 = -D/10 28 cm and we 
have ignored cosmological redshift effects. Fig. ^2] depicts 
the instantaneous spectrum in this case. 
For a fully radiative evolution we find: 



F„ 



1 1 v in 13 tlr, c -3/2 p -4/7 4/7 -13/14 -2/7 

1.3 x 10 Hz e B ii 52 7 100 n 1 t d , 
1.2 x 10 14 

l3. 



Hz e 1/2 e 2 £; 47 V 47 V l7l V 1277 
nz e B e e £/ 52 7 100 tj 



4.5 x 10^JeV 2 J E; 5 8 / 7 7 ro 8 o /7 n5 714 J D 28 2 ^ 3/ (l33) 



where we have scaled the initial Lorentz factor of the ejecta 
by a factor of 100: 7100 = 70/IOO. These instantaneous 
spectra are also shown in Fig. |2^. 

9.3.1. Light Curves 

The light curves at a given frequency depend on the tem- 
poral evolution of the break frequencies v m and v c and the 
peak power N e P syn ("f e ^ min ) (see Eq. |6|). These depend, 
in turn, on how 7 and N e scale as a function of t. 

The spectra presented in Fig. |2^ show the positions of 
v c and v m for typical parameters. In both the adiabatic 
and radiative cases v c decreases more slowly with time 
than v m . Therefore, at sufficiently early times we have 
v c < v m , i.e. fast cooling. At late times we have v c > v m , 
i.e., slow cooling. The transit ion b etween the two occurs 
when v c = v m at tf s (see Eq. |126| ). At t = i/ s , the spec- 
trum changes from fast cooling (Fig. |22]a) to slow cooling 
(Fig. p2|b). In addition, if e e w 1, the hydrodynamical 
evolution changes from radiative to adiabatic. However, if 
e e <C 1, the evolution remains adiabatic throughout. 

Once we know how the break frequencies, v c , v m , and 
the peak flux F v _ max vary with time, we can calculate the 
light curve. Consider a fixed frequenc y (e. g. v = 10 1 5 ^i5 



Hz). From the first two equations in (132) and (132) we 
see that there are two critical times, t c and t m , when the 
break frequencies, v c and v m , cross the observed frequency 
v: 



tc = 



7.3 x 10- 6 days e^E^ufv^ , 



ad, 



2.7 x 10- 7 days e B 21/i E^ loa n; 13/4 ^ /2 , rad, 

(134) 

ad, 



_ / 0.69 days e'fef 3 2$ V' 3 



[ 0.29 days eJ^el^E^ 3 ^ 3 ^ 12 ^ 1/24 , ra. 

(135) 

There are only two possible orderings of the three critical 
times, t c , t m , tf s , namely t fs > t m > t c and tj s < t m < t c . 
We define the critical frequency, vq = v c (tf s ) — v m {tf s ): 



VQ = 



f 1.8 x lQ ll e B 5/2 e- l E^n^ /2 Hz, 



8.5 x 10^e 



12 -19/10 -2/5 p-4/5 4/5 -11/10 



-^52 7l00 n l 



ad, 

Hz, rad. 
(136) 



When v > vq, we have tt s > t m > t c and we refer to the 
corresponding light curve as the high frequency light curve. 
Similarly, when v < i/q, we have tf s < t m < t c , and we 
obtain the low frequency light curve. 

Fig. |27| a depicts a typical high frequency light curve. 
At early times the electrons cool fast and v < v m and 
v < v c . Ignoring self absorption, the situation corre- 
sponds to segment B in Fig. |2^, and the flux varies 
as F v ~ F v , max (y /^ c ) 173 '■ If the evolution is adiabatic, 
F Vi max is constant, and F v ~ t 1 / 6 . In the radiative case, 
F-u.max ~ t -377 and F v ~ t -1 / 3 . The scalings in the other 
segments, which correspond to C, D, H in Fig. |22|, can be 
derived in a similar fashion and are shown in Fig. |27]a. 

Fig. p7j b shows the low frequency light curve, corre- 
sponding to v < vq. In this case, there are four phases in 
the light curve, corresponding to segments B, F, G and H. 
The time dependences of the flux are indicated on the plot 
for both the adiabatic and the radiative cases. 

For a relativistic electron distribution with a power dis- 
tribution 7~ p the uppermost spectral part behaves like 
j,-p/2 The corresponding temporal index (for adiabatic 
hydrodynamics) is — 3p/4. In terms of the spectral index 
a, this yields the relation F v cx £( 1-3Q! )/2. Alternatively for 
slow cooling there is also another frequency range (between 
v m and v c ) for which the spectrum is given by ^~(p _1 )/2 
and the temporal decay is — 3(p — l)/4. Now we have 
F v oc £~ 3q / 2 , Note that in both cases there is a specific re- 
lation between the spectral index and the temporal index 
which could be tested by observations. 



9.3.2. Parameter Fitting for GRBs from Afterglow 
Observation and GRB970508 



Shortly after the observation of GRB970228 Meszaros et 
al. |22| showed that the decline in the intensity in X-ray 
and several visual bands (from B to K) fit the afterglow 
model well (see fig. ^8|) . The previous discussion indicates 
that this agreement shows that the high energy tail (or 
late time behavior) is produced by a synchrotron emission 
from a power law distribution. 

There are much more data on the afterglow of 
GRB970508. The light curves in the different optical 
bands generally peak around two days. There is a rather 
steep rise before the peak which is followed by a long power 
law decay (see figs. In the optical band the ob- 

served power law decay for GRB970508 is -1.141 ± 0.014 
1 151]. This implies for an adiabatic slow cooling model 
(which Meszaros et al. Q use) a spectral index a — 
—0.761 ± 0.009. H owev er, the observed spectral index is 
a = -1.12 ± 0.04 p60| . A fast cooling model,for which 
the spectral index is p/2 and the temporal behavior is 
3p/4— 1/2, fits the data better as the temporal power law 
implies that p — 2.188 ± 0.019 while the spec tral index 
implies, consistently, p = 2.24 ± 0.08 (246| ||| (see Fig. 
29 . Unfortunately, at present this fit does not tell us much 
about the nature of the hydrodynamical processes and the 
slowing down. 

Using both the optical and the radio data on can try to 
fit the whole spectrum and to obtain the unknown parame- 
ters that determine the fireball evolution 248 , 249 1. Wijers 
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and Galama p4q| have attempted to do so using the spec- 
trum of GRB970508. They have obtained a reason able 
set of parameters. However, more detailed analysis [249 
reveals that the solution is very sensitive to assumptions 
made on how to fit the observational data to the theoret- 
ical curve. Moreover, the initial phase of the light curve 
of GRB970508 does not fit any of the theoretical curves. 
This suggests that at least initially an additional process 
might be taking place. Because of the inability to obtain 
a good fit for this initial phase there is a large uncertainty 
in the parameters obtained in this way. 

No deviation in the observed decaying light curve from 
a single power law was observed for GRB970508, until it 
faded below the level of the surrounding nebula. This sug- 
gests that there was no significant beaming in this case. If 
the outflow is in the form of a jet the temporal behavior 
will change dras tically when the opening angle of the jet 
equals 1/7 J258[ . 

9.4. New Puzzles from Afterglow observations 

Afterglow observations fit well the fireball picture that 
was developed for explaining the GRB phenomena. The 
available data is not good enough to distinguish between 
different specific models. But in the future we expect to be 
able to distinguish between those models and even to be 
able to determine the parameters of the burst E and 70 (if 
the data is taken early enough), the surrounding ISM den- 
sity and the intrinsic parameters of the relativistic shock 
e e , £b and p. Still the current data is sufficient to raise 
new puzzles and present us with new questions. 

• Why afterglow accompany some GRBs and 
not others? 

X-ray, Optical and radio afterglows have been ob- 
served in some bursts but not in others. According 
to the current model afterglow is produces when the 
ejecta that produced the GRB is shocked by the sur- 
rounding matter. Possible explanations to this puz- 
zle invoke environmental effects. A detectable after- 
glow might be generated efficiently in some range of 
ISM densities and inefficiently in another. High ISM 
densities would slow down of the ejecta more rapidly. 
This could make some afterglows detectable and oth- 
ers undetectable. ISM absorption is another alterna- 
tive. While most interstellar environments are opti- 
cally thin to gamma-rays high density ISM regions 
can absorb and attenuate efficiently x-rays and op- 
tical radiation. 

• Jets and the Energy of GRB971214 

How can we explain the 10 53 ergs required for 
isotropic emission in GRB971214? As we discuss 
in the next section this amount is marginal for most 
models that are based on the formation of a com- 
pact source. This problem can be resolved if we in- 
voke beaming, with 9 ~ 0.1. However, such beaming 
would results in a break at the light curve when the 
local Lorentz factor would reach a value of 1/9. Such 
a break was not seen in other afterglows for which 
there are good data. Note that recently Perna & 
Loeb 126311 inferred from the lack of radio transients 



that GRB beamns cannot be very narrow. If typi- 
cal GRBs are beamed, the beam width 9 should be 
larger than 6°. 

• GRB980425 and SN1998bw 

SN1998bw (and the associated GRB980425) is a fac- 
tor of a hundred nearer than a typical GRB (which 
are expected to be at z ~ 1). The corresponding 
(isotropic) gamma-ray energy, ~ 5x 10 47 ergs, is four 
order of magnitude lower than a regular burst. This 
can be in agreement with the peak flux distribution 
only if the bursts with such a low luminosity compose 
a very small fraction of GRBs. This leads naturally 
to the question is there an observational coincidence 
between GRBs and SNs? To which there are con- 



flicting answers [264, 265, |266| 



10. MODELS OF THE INNER ENGINE 

We turn now to the most difficult part: the nature of the 
beast that produces the GRB-modeling of the Inner En- 
gine. We examine a few general considerations in section 



10.1 and then we turn to the Binary Neutron star merger 



model in 10.2 



10.1. The "Inner Engine" 

The fireball model is based on an "inner engine" that 
supplies the energy and accelerates the baryons. This "en- 
gine" is well hidden from direct observations and it is im- 
possible to determine what is it from current observations. 
Unfortunately, the discovery of afterglow does not shed ad- 
ditional direct light on this issue. However it adds some 
indirect evidence from the association of the location of 
the bursts in star forming regions. 

Once the cosmological origin of GRBs was established 
we had two direct clues on the nature of the "inner en- 
gine": the rate and the energy output. GRBs occur at a 
rate of about one per 10 6 years per galaxy ]5(| and the to- 
tal energy is ~ 10 52 ergs. These estimates assume isotropic 
emission. Beaming with an angle 9 changes these estimates 
by a factor Air/9 2 in the rate and 9 2 /4ir in the total energy 
involved. These estimates are also based on the assump- 
tion that the burst rate does not vary with cosmic time. 
The ob servations that GR B hosts are star forming galax- 
ies pi, p|, pM p|, |28| indicates that the rate of GRBs 



may follow the star formation rate |193|, 194, 195]. In this 
case the bursts are further and they take place at a lower 
rate and have significantly higher energy output. 

The fireball model poses an additional constraint: the 
inner engine should be capable of accelerating ~ 10 _5 M Q 
to relativistic energies. One can imagine various scenar- 
ios in which 10 52 ergs are generated within a short time. 
The requirement that this energy should be converted to 
a relativistic flow is much more difficult as it requires a 
"clean" system with a very low but non zero baryonic load. 
This requirement suggests a preference for models based 
on electromagnetic energy transfer or electromagnetic en- 
ergy generatio n as t hese cou ld m ore naturally satisfy this 
condition (see |267[ p24, 



^ |27]]). Paczyhski M] has re- 
cently suggested a unique nydrodynamical modelin which 
10 54 ergs are dumped into an atmosphere with a decreasing 
dens ity profile. This is a cosmological variant of Colgate's 
|262[| galactic model. This would lead to an acceleration 
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of fewer and fewer baryons and eventually to a relativistic 
velocities. Overall one could say that the "baryonic load" 
problem is presently the most bothersome open question 
in the "fireball model" . 

The recent realization that energy conversion is most 
likely via internal shocks rather than via external shocks 
provides additional information about the inner engine: 
The relativistic flow must be irregular (to produce the in- 
ternal shocks), it must be variable on a short time scale 
(as this time scale is seen in the variability of the bursts), 
and it must be active for up to a few hundred seconds and 
possibly much longer |}6| - as this determines the observed 
duration of the burst. These requirements rule out all ex- 
plosive models. The engine must be compact (~ 10 7 cm) 
to produce the observed variability and it must operate for 
a million light crossing times to produce a few hundred- 
second signals. 

There are more than a hundred GRB models p68|. At 
a certain stage, before BATSE, there were probably more 
models than observed bursts. Most of these models are, 
however, galactic and those have been ruled out if we ac- 
cept the cosmological origin of GRBs. This leaves a rather 
modest list of viable GRB mod els: binar y n eutron st ar 
mergers - NS 2 Ms - M (see also p6f, |70[ $Ug% W% ), 
failed supernova |273| ] , white dwarf collapse | p67| and hy- 
pernova pq] . All these are based on the formation of a 
compact object of one type or another and the release of 
its binding energy. With a binding energy of ~ 5 x 10 53 
ergs or higher, all these models have, in principle, enough 
energy to power a GRB. However they face similar diffi- 
culties in channeling enough energy to a relativistic flow. 
This would be particularly difficult if indeed 10 53 ergs are 
needed, as some recent burst have indicated. Paczyhski's 
hypernova is an exception as in this model all the energy is 
channeled initially to a non-relativistic flow and only later 
a small fraction of it is converted to relativistic baryons. 
All these models are consistent with the possibility that 
GRBs are associated with star forming regions as the life 
time of massive stars is quite short and even the typical 
life time of a neutron star binary (~ 10 8 yr) is sufficiently 
short to allow for this coincidence. - 

Other models are based on an association of GRBs with 
massive black holes as socia ted with Quasars or AGNs in 
galactic centers (e.g. [274]). These are ruled out as all 
GRBs with optical afterglow are not associated with such 
objects. Furthermore, such objects do not appear in othe r 
small GRB error boxes searched by Schaefer et al. [ 126 1 . 
From a theoretical point of view it is difficult to explain 
the observed energy and time scales with such objects. 

10.2. NS 2 Ms: Binary Neutron Star Mergers 



Binary neutron star mergers (NS 2 Ms) |35| or , wi th a 
small variant: neutron star-black hole mergers ]271| are 
probably the best candidates for GRB sources. These 
mergers take place because of the decay of the binary 
orbits due to gravitational radiation emission. A NS 2 M 
results, most likely, in a rotating bla ck hole [276|. The 
process releases « 5 x 10 53 ergs [277]. Most of this en- 



The discovery of othe r binary pulsars, and in particular 
of PSR 1534+12 (27|, has shown that PSR 1913+16 is 
not unique and that such systems are common. These 
observations suggest that NS 2 Ms take place at a rate of 
10~ 6 events per year per galaxy EG, 



52 1 . This rate 



is comparable to the simple estimate of the GRB event 
rate (assu min g no beaming and no cosmic evolution of the 



rate) M |272|, |183| 



It has been suggested [279, 280 1 that most neutron star 
binaries are born with very close orbits an d hence with 
very short lifetimes (see however, ]281| , |282fl ). If this idea 
is correct, then the merger rate will be much higher. This 
will destroy, of course, the nice agreement between the 
rates of GRBs and NS 2 Ms. Consistency can be restored 
if we invoke beaming, which might even be advantageous 
as far as the energy budget is concerned. Unfortunately, 
the short lifetime of those systems, which is the essence of 
this idea means that at any given moment of time there 
are only about a hundred such systems in the Galaxy 
(compared to about 10 5 wider neutron star binaries). This 
makes it very hard to confirm or rule out this speculation. 
We should be extremely lucky to detect such a system. 

It is not clear yet how NS 2 Ms form. The question is 
how does the system survive the second supernova event? 
The binary system will be disrupted it this explosion ejects 
more than half of its total mass. There are two competing 
scenarios for the formation of NS 2 Ms. In one scenario the 
first neutron star that forms sinks into the envelope of its 
giant companion and its motion within this envelope lead 
to a strong wind that carries away most of the secondary's 
mass. When the secondary reaches core collapse it has 
only a small envelope and the total mass ejected is rather 
small. In a second scenario the second supernova explosion 
is asymmetric. The asymmetric explosion gives a velocity 
of a few hundred km/sec to the newborn neutron star. In a 
fraction of the cases this velocity is in the right direction to 
keep the binary together. Such a binary s ystem wil l hav e 
a comparable center of mass velocity [EM 283, 284, |285|. 



ergy escapes as neutrinos and gravitational radiation, but 
a small fraction of this energy suffices to power a GRB. 
The discovery of the famous binary pulsar PSR 191 3+1 6 
p8| demonstrated that this decay is taking place [ 275 1 . 



This second scenario has several advantages. First it ex- 
plains both the existence of binar y ne u tron stars and the 
existence of high velocity pulsars |283| , 286 1. Second, and 
more relevant to GRBs, with these kick velocities these bi- 
naries could escape from their parent galaxy, provided that 
this galaxy is small enough. Such escaping binaries will 
travel a distance of ~ 200 kpc (v/200 km/sec) (T/10 9 yr) 
before they merge. The GRB will occur when the system 
is at a distance of the order of hundred kpc from the par- 
ent galaxy. Clearly there is no "no host" problem in this 
case p8[ ]. 

While a NS 2 M has enough energy available to power 
a GRB it is not clear how the GRB is produced. A 
central question is, of course, how does a NS 2 M gener- 
ate the relativistic wind required to power a GRB. Most 
of the binding ener gy (w hich is around 5 x 10 53 ergs es- 
capes as neutrinos [277]. Eichler et al. |35| suggested 
that about one thousandth of these neutrinos annihilate 
and produce pairs that in turn produce gamma-rays via 
vv — > e + e~ — > 77. This idea was criticiz ed o n several 
grounds by different authors. Jaroszynksi ]288[ | pointed 
out that a large fraction of the neutrinos wil l be swallowed 
by the black hole that forms. Davies et al. ]276| ] and Ruf- 
fert & Janka [289, 290, 291] who simulated neutron star 
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mergers suggested that the central object won't be warm 
enough to produce a significant neutrino flux because the 
merger is nearly adiabatic |69| . The neutrinos are also 
emitted over a diffusion time of several seconds, too long 
to explain the rapid variations observed in GRB Jr9[ , but 
to short to e xpla i n th e observed GRB durations. Wilson 
& Mathews [292, 292] included approximate general rela- 
tivistic effects in a numerical simulation of a neutron star 
merger. They found that the neutron stars collapse to a 
single black hole before they collide with each other. This 
again will suppress the neutrino emission from the merger. 
However, the approximation that they have used has been 
criticized by various authors and it is not clear yet that 
the results are valid. Others suggested that the neutrino 
wind will carry too many baryons. However, it seems that 
the most severe problem with this model stems from the 
fact that the prompt neutrino burst could produce only a 
single smooth pulse. This explosive burst is incompatible 
with the internal shocks scenario. 

An alternative source of energy within the NS 2 M is the 
accretio n pow er of a disk that forms around the black hole 
p§| , |S9[ 227] . Various numeri cal simulations of neutron 
star mergers |276j [289] |290[ p9l| find that a ~ O.LM© forms 
around the central black hole. Accretion of this disk on the 
central black hole may take a few dozen seconds |3!|. It 
may produce the wind needed to produce internal shocks 
that could produce, in turn a GRB. 

How can one prove or disprove this, or any other, GRB 
model? Theoretical studies concerning specific details of 
the model can, of course, make it more or less appeal- 
ing. But in view of the fact that the observed radiation 
emerges from a distant region which is very far from the 
inner "engine" I doubt if this will ever be sufficient. It 
seems that the only way to confirm any GRB model will 
be via detecting in time-coincidence another astronomi- 
cal phenomenon, whose source could be identified with 
certainty. Unfortunately while the recent afterglow obser- 
vations take us closer to this target they do not tell us 
what are the sources of GRBs. We still have to search for 
additional signals. 

NS 2 Ms have two accompanying signals, a neutrino sig- 
nal and a gravitational radiation signal. Both signals are 
extremely difficult to detect. The neutrino signal could be 
emitted by some of the other sources that are based on a 
core collapse. Furthermore, with present technology de- 
tection of neutrino signals from a cosmological distance is 
impossible. On the other hand the gravitational radiation 
signal has a unique characteristic form. This provides a 
clear prediction of coincidence that could be proved or fal- 
sified sometime in the not too distant future when suitable 
gravitational radiation detectors will become operational. 

10.3. Binary Neutron Stars vs. Black Hole - Neutron 
star Mergers 

We have grouped together binary neutron star mergers 
and black- hole neutron star mergers. At present several 
neutron star binaries are known while no black hole - neu- 
tron star binary was found. Still on theoretic al g rounds 
one should expect a similar rate for both events |5C| ] . Some 
even suggest that there are more black hole - neutron 



star binaries than neutron star- neutron star binaries [282|. 
There is a lot of similarity between the two processes which 



are both driven by gravitational radiation emission and 
both result in a single black hole. First, unless the mass of 
the black hole is of rather small the neutron star will not 
be tidally disrupted before it is captured by the black hole. 
Even if such a tidal disruption will take place then while 
in the binary neutron star merger we expect a collision, 
in the black hole - neutron star merger we expect at most 
a tidal disruption followed by infall of the debris on the 
black hole. 

This could lead to a situation in which one of the 
two events will produce a GRB and the other will not. 
Presently it is too difficult to speculate which of the two 
is the right one. One should recall, however, that there is 
a marked difference between the gravitational signature of 
those events and thus hopefully when we discover a coinci- 
dence between a GRB and a gravitational radiation signal 
we would also be able to find which of the two mechanism 
is the right one. 



11. OTHER RELATED PHENOMENA 

It is quite likely that other particles (in addition to 7- 
rays) are emitted in these events. Let f XJ be the ratio of 
energy emitted in other particles relative to 7-rays. These 
particles will appear as a burst accompanying the GRB. 
The total fluence of a "typical" GRB observed by BATSE, 
F 7 is 10~ 7 ergs/cm 2 , and the fluence of a "strong" burst is 
about hundred times larger. Therefore we should expect 
accompanying bursts with typical fluences of: 
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(137) 

where E x is the energy of these particles. This burst will 
be spread in time and delayed relative to the GRB if the 
particles do not move at the speed of light. Relativistic 
time delay will be significant (larger than 10 seconds) if 
the particles are not massless and their Lorentz factor is 
smaller than 10 8 ! similarly a deflection angle of 10 -8 will 
cause a significant time delay. 

In addition to the prompt burst we should expect a con- 
tinuous background of these particles. With one 10 51 ergs 
GRB per 10 6 years per galaxy we expect ~ 10 4 events 
per galaxy in a Hubble time (provided of course that the 
event rate is constant in time). This corresponds to a 
background flux of 
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(138) 



For any specific particle that could be produced one 
should calculate the ratio j xl and then compare the ex- 
pected fluxes with fluxes from other sources and with the 
capabilities of current detectors. One should distinguish 
between two types of predictions: (i) Predictions of the 
gene ric fireball m odel which inc lude low energy cosmic rays 
122C|1, UCHERs [pi E>95], E>96|] and high energy neutrinos 



297 ] and (ii) Predictions of specific models and in partic- 
ular the NS 2 M model. These incl ude lo w energy neutrinos 
[277] and gravitational waves pj], 301 1 . 
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11.1. Cosmic Rays 



Already in 1990, Shemi & Piran |22C|] pointed out 
that fireball model is closely related to Cosmic Rays. 
A "standard" fireball model involved the acceleration of 
~ 10~ 7 Af Q of baryons to a typical energy of lOOGeV per 
baryon. Protons that leak out of the fireball will become 
low energy cosmic rays. However, a comparison of the 
GRB rate (one per 10 6 years per galaxy) with the ob- 
served low energy cosmic rays flux, suggests that even if 
fcR--y ~ 1 this will amount only to 1% to 10% of the 
observed cosmic ray flux at these energies. Cosmic rays 
are believed to be produced by SNRs. Since supernovae 
are ten thousand times more frequent than GRBs, unless 
GRBs are much more efficient in producing Cosmic Rays 
in some specific energy range their contribution will be 
swamped by the SNR contribution. 



5 — lOMeV) neutrinos |p77| . The total energy is quite 
large ~ a few xl0 53 ergs, leading to /i ow energy v--y ~ 10. 
However, this neutrino signal will be quite similar to a su- 
pernova neutrino signal, and at present only galactic SN 
neutrinos can be detected. Supernovae are ten thousand 
times more frequent then GRBs and therefore low energy 
neutrinos associated with GRB constitute an insignificant 
contribution to the background at this energy range. 



11.6. Black Holes 

An NS 2 M results, inevitably, in a black hole [^76| . Thus 
a direct implication of the NS 2 M model is that GRBs sig- 
nal to us (indirectly) that a black hole has just (with the 
appropriate time of flight in mind) formed. 



11.2. UCHERs - Ultra High Energy Cosmic Rays 



Waxman [ 294 1 and Vietri [ 295 1 have shown that the ob- 
served flux of UCHERs (above 10 19 eV) is consistent with 
the idea that these are produced by the fireball sh ocks 
provided that fucHERs-^ ~ 1- Milgrom & Usov [296 



pointed out that the error boxes of the two highest en- 
ergy UCHERs contain strong GRBs - suggesting an as- 
sociation between the two phenomena. The relativistic 
fireball shocks that appear in GRBs are among the few 
astronomical objects that satisfy the con ditions for shock 
acceleration of UCHERs. Waxman [29? ] has shown that 
the spectrum of UCHERs is consistent with the one ex- 
pected from Fermi acceleration within those shocks. 

11.3. High Energy Neutrinos 



Waxman & Bahcall |297| suggested that collisions be- 
tween protons and photons within the relativistic fireball 
shocks produce pions. These pions produce high energy 

neutrinos with E v ~ 10 14 eV and /high energy v f > 0.1. 

The flux of these neutrinos is comparable to the flux of at- 
mospheric neutrinos but those will be correlated with the 
position of strong GRBs. This signal might be detected in 
future km 2 size neutrino detectors. 

11.4. Gravitational Waves 

If GRBs are associated with NS 2 Ms then they will be 
associated with gravitational waves and low energy neu- 
trinos. The spiraling in phase of a NS 2 M produces a clean 
chirping gravitational rad iation signal. This si gnal is the 
prime target of LIGO |§99| and VIRGO, the two Q large 
interferometers that are build now in the USA and in Eu- 
rope. The observational scheme of these detectors is heav- 
ily dependent on digging deeply into the noise. Kocha- 
neck & Piran [ 301 1 suggested that a coincidence between a 
chirping gravitational radiation signal from a neutron star 
merger and a GRB could enhance greatly the statistical 
significance of the detection of the gravitational radiation 
signal. At the same time this will also verify the NS 2 M 
GRB model. 



11.5. Low Energy Neutrinos 

Most of the energy generated in any core collpase event 
and in particular in NS 2 M is released as low energy (~ 



12. COSMOLOGICAL IMPLICATIONS 

Cosmological GRBs seem to be a relatively homoge- 
neous population of sources with a narrow luminosity func- 
tion (the pe ak lumino sity of GRBs varies by less than a 
factor of 10 |183|, |186|) that is lo cated at relatively high 
redshifts T 
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The universe a nd ou r 
Galaxy are transparent to MeV 7-rays (see e.g. |305[ ) . 
Hence GRBs constitute a unique homogeneous population 
of sources which does not suffer from any angular distor- 
tion due to absorption by the Galaxy or by any other ob- 
ject. Could GRBs be the holy grail of Cosmology and 
provide us with the standard candles needed to determine 
the cosmological parameters Hq, fi, and A? Lacking any 
spectral feature, there is no indication of the redshift of in- 
dividual bursts. The available number vs. peak luminosity 
distribution is is not suitable to distinguish between differ- 
ent cosmological models even when the sour ces are perfect 
standard candles with no source evolution [183]. 

The situation might be different if optical afterglow ob- 
servations would yield an independent redshift measure- 
ment of a large number of bursts. If the GRB luminosity 
function is narrow enough this might allow us, in the fu- 
ture, to determine the cosmological closure parameter SI 
using a peak-flux vs. red-shift diagram (or the equiva- 
lent more common magnitude - red-shift diagram). For 
example a hundred bursts with a measured z are needed 
to es timate $1 with an accuracy of ctq = 0.2, if o~l/L = 1 
|20|. 

Currently, the rate of detection of bursts with counter- 
parts is a few per year and of those detected until now 
only two have a measured red-shift. This rate is far too 
low for any cosmological measurement. However, there is 
an enormous potential for improvements. For example, 
systematic measurements of the red-shift of all bursts ob- 
served by BATSE (w 300 per year) would yield an indepen- 
dent estimate of f2, with ctji = 0.1, even if the luminosity 
function is wide, (er^/L = 0.9), within one year. 

Direct redshift measurements would also enable us to 
dete rmine the cosmological evolution of the rate of GRBs 
201]. Most current cosmological GRB models suggest that 
the GRB rate follows ( with a rather short time delay) the 
rate of star formation [306]. Consequently measurements 
of the rate of GRBs as a function of the red-shift will 
provide an independent tool to study star formation and 
galactic evolution. 
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It is also expected that the bursts' sources follow the 
matter distribution. Then GRBs can map the large scale 
structure of the Universe on scales that c ann ot be spanned 
directly otherwise. Lamb & Quashnock |307| have pointed 
out that a population of several thousand cosmological 
bursts should show angular deviations from isotropy on a 
scale of a few degrees. This would immediately lead to new 
interesting cosmological limits. So far there is no detected 
anisotropy in the 1112 bursts of the BATSE 3B catalog 
1 308]. But the potential of this population is clear and 
quite promising. A more ambitious project would be to 
measure the multipole moments of the GRB distribution 
and from this to estimate cosmological parameters p8| . 
However, it seems that too many bursts are required to 
overcome the signal to noise ratio in such measurements. 

GRBs can also serve to explore cosmology as a back- 
ground pop ulation which could be lensed by foreground 
objects |5J|. While standard gravitational lensed object 
appears as several images of the same objects, the low an- 
gular resolution of GRB detectors is insufficient to distin- 
guish between the positions of different images of a lensed 
GRB. However, the time delay along the different lines 
of sight of a gravitationally lensed burst will cause such 
a burst to appear as repeated bursts with the same time 
profile but different intens ities from practically the same 
position on the sky. Mao [309] estimated that the proba- 
bility for lensing of a GRB by a regular foreground galaxy 
is 0.04 %-0 .4%. Hence the lack of a confirmed lensed event 
so far ]310| is not problematic yet. In the future, the statis- 
tics of lensed bursts could probe the nature of the lensing 
objects and the dark matter of the Universe |i3]. The fact 
that no lensed bursts have been detected so far is suffi- 
cient to rule out a critical density (tt = 1) of 10 6 5 M Q to 
10 81 M Q black holes [ p5| . Truly, this was not the leading 
candidate for cosmological dark matter. Still this result is 
a demonstration of the power of this technique and the po- 
tential of GRB lensing. The statistics of lensing depends 
on the distance to the lensed objects which is quite un- 
certain at present. The detection of a significant number 
of counterparts whose red-shift could be measured would 
improve significantly this technique as well. 

13. SUMMARY AND CONCLUSIONS 

Some thirty years after the discovery of GRBs a generic 
GRB model is beginning to emerge. The observations of 
isotropy, peak flux distribution and time dilation indicated 
that GRBs are cosmological. The measurement of a red- 
shift provided a final confirmation for this idea. All cos- 
mological models are based on the fireball mode. The 
discovery of the afterglow confirmed this general scheme. 

The Fireball-Internal-External Shocks model seems to 
have the necessary ingredients to explain the observations. 
Relativistic motion, which is the key component of all Fire- 
ball Model provided the solution for the compactness prob- 
lem. The existence of such motion was confirmed by the 
radio afterglow observations in GRB970508. Energy con- 
version via internal shocks can produce the observed highly 
variable light curves while the external shock model agrees, 
at least qualitatively, with afterglow observations. There 
are some indications that the two kind of shocks might 
combine and operate within the GRB itself producing two 
different components of the signal. This fireball model has 



some fascinating immediate implications on accompany- 
ing UCHER and high energy neutrino signals. An obser- 
vations of these phenomena in coincidence with a GRB 
could provide a final confirmation of this model. 

In spite of this progress we are still far from a complete 
solution. There are many open questions that has to be 
resolved. Within the internal-external shocks model there 
is a nagging efficiency problem in conversion of the initial 
kinetic energy to the observed radiation. If the overall effi- 
ciency is too low the initial energy required might be larger 
than 10 53 ergs and it is difficult to imagine a source that 
could provide so much energy. Beaming might provide a 
solution to this energy crisis. However, so far there is no in- 
dication for the corresponding break in the afterglow light 
curve, which is essential in any relativistic beaming model 
when 7 ~ These last two facts might be consistent 

if GRB970508 took place within a very low density ISM 
- an issue that should be explored further. Afterglow ob- 
servations agree qualitatively, but not quantitatively with 
the model. Better observations and more detailed theo- 
retical modeling are needed. Another nagging open ques- 
tion is what determines the appearance of afterglow. Why 
there was no X-ray afterglow in the very strong 970111? 
Why was optical afterglow observed in GRB970228 and in 
GRB970508 but it was not seen in others (in particular in 
GRB970828 @)? Finally we turn to the GRB itself and 
wonder why is the observed radiation always in the soft 7- 
ray band? Is there an observational bias? Are there other 
bursts that are not observed by current instruments? If 
there are none and we do observer all or most of the bursts 
why is the emitted radiation always in the soft gamma-ray 
range? Why it is insensitive to a likely variability in the 
Lorentz factors of the relativistic flow and to variability of 
other parameters in the model. 

While there are many open questions concerning the 
fireball and the radiation emitting regions the first and 
foremost open question concerning GRBs is what are the 
inner engines that power GRBs? In spite of all the recent 
progress we still don't know what produces GRBs. My 
personal impression is that binary neutron mergers are 
the best candidates. But other models that are based on 
the formation of a compact object and release a significant 
amount of its binding energy on a short time scale are also 
viable. A nagging question in all these models is what pro- 
duces the the "observed" ultra-relativistic flow? How are 
~ 10 _5 M Q of baryons accelerated to an ultra-relativistic 
velocity with 7 ~ 100 or larger? Why is the baryonic load 
so low? Why isn't it lower? There is no simple model for 
that. An ingenious theoretical idea is clearly needed here. 

However, I believe that theoretical reasoning won't be 
enough and only observations can provide a final resolu- 
tion of the questions what is are the sources of GRBs? 
The binary neutron star merger model has one specific ob- 
servational prediction: A coincidence between a (near by 
and therefore strong) GRB and a characteristic gravita- 
tional radiation signal. Luckily these events have a unique 
gravitational radiation signature. The detection of these 
gravitational radiation events is the prime target of three 
gravitational radiation detector that are being built now. 
Hopefully they will become operational within the next 
decade and their observations might confirm or rule out 
this model. Such predictions, of an independently ob- 
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served phenomena are clearly needed for all other com- 
peting models. 

GRBs seem to be the most relativistic phenomenon dis- 
covered so far. They involve a macroscopic relativistic mo- 
tion not found elsewhere before. As cosmological objects 
they display numerous relativistic cosmological phenom- 
ena. According the the NS 2 M model they are associated 
with the best sources for gravitational radiation emission 



and more than that they signal, though in directly, the 
formation of a new black hole. 
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Fig. 1. — Total number of counts vs. time for several bursts from the BATSE Catalogue. Note the large diversity of 
temporal structure observed. 



GRB1676 




time since triger (s) 

Fig. 2. — Counts vs. time for BATSE burst 1676. The bursts lasted T ~ 60 sec and it had peaks of width 5T ~ 1 sec, 
leading to J\f w 60. 
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Fig. 3. — Observed spectrum of BATSE' burst 228. 
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Fig. 4. — N(H) - the number of bursts with hardness , H , in the Band et al. faj/ sample (dashed- dotted line) and in the 
Cohen et al. sample (solid line) togather with a theoretical fit of a distribution above H = 120KeV with 7 ~ —0.5 (a 
slowly decreasing numbers of GRBs per decade of hardness ) 
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Fig. 5. — Decay of the X-ray afterglow from GRB970228, from Jlq/. Shown is source flux at the 2-10KeV range. The 
data is fitted with a power law t~ 132 . 
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Fig. 6.— Decay of the optical afterglow in GRB070228, GRB9700508 and GRB971214. A clear power law decay can be 
seen in all cases. 
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Fig. 9. — The likelihood function (levels 33%, 10%, 3.3% 1% etc..) in the (/3,L) plane for standard candles a — 1.5, 
= 1, and evolution given by p(z) = + . Superimposed on this map is the luminosity of GRB970508, (solid curve), 
GRB971214 (dashed curve) and GRB980703 (dahsed- dotted curv e). We have used /175 = 1. 

From pOll . 
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Fig. 10. — The observed long burst peak flux distribution and three theoretical cosmological distributions with Q = 1, A = 0, 
a = —1.5, standard candles and no source evolution: L ~ 3.4 • 10 50 ergs/sec (solid line: best fit), L — 7.2 • 10 50 ergs/sec 
(dashed line: lower 1% bound), L — 1.4- 10 50 ergs/sec (dashed- dotted line: upper 1% bound) 
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Fig. 11. — Radiation from a relativistic beam with a width 9. Each observer will detect radiation only from a very narrow 



beam with a width T 



The overall angular size of the observed phenomenon can vary, however, with T 



< 



< 4ir. 




Fig. 12. — Fireball evolution from its initial formation at rest to the final Newtonian Sedov solution. The energy extraction 
is due to the interaction with the ISM via a relativistic forward shock and a Newtonian reverse shock. We have used for 
this calculations £ = 43, Eq — 10 52 [erg], 70 = 50 Rq = 3 x 10 10 [cm]. Shown are the average value of the Lor entz factor 
(thick solid line), the value at th e forward shock (thin solid line), the maximal value (dotted line) and an analytic estimate 
(dashed dotted line). From (22i]. 
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Fig. 13. — Fireball evolufrom its initial formation at rest to the final Newtonian Sedov solution. The energy extraction is 
due to the interaction with the ISM via relativistic forward and reverse shocks. The parameters for this computation are: 
£ = 0.1, Eq = 10 52 [erg], 70 = 10 4 , i?o = 4.3 x 10 9 [cm]. Shown are the average value of the Lorentz factor (thick solid 
line), the value at t he f orward shock (thin solid line), the maximal value (dotted line) and an analytic estimate (dashed 
dotted line). From 





Fig. 14. — Different time scales in terms of the arrival time of four photons: tA,ts, tc> ond tjj. T rac ii a i — tc ~ ^a; 



T, 



angular 



tD - tA, A/c = t B ~t A . 
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Fig. 15. — The deceleration radius R e and the Lorentz factor of the shocked shell 7 e as functions of the initial Lorentz 
factor 7, for a shell of fixed width A = 3 x 10 12 cm. For low values ofj, the shocked material moves with Lorentz factor 
7e ~ 7. However as 7 increases the reverse shock becomes relativistic reducing significantly the Lorentz factor 7 e < 7. 
This phenomena prevents the "external shock model" from being Type-II. 




Fig. 16. — A very narrow jet of angular size considerably smaller than 7" 1 , for which the angular spreading problem does 
not exist. The duration of the burst is determined by the deceleration distance AR e , while the angular time is assumed 
small. The variability could now be explained by either variability in the source which leads to a pulsed jet (a) or by a 
uniform jet interacting with an irregular LSM (b ) 
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Fig. 17. — A shell with angular size 7" 1 (the angular size is highly exaggerated). The spherical symmetry is broken by 
the presence of bubbles in the ISM. The relative angular size of the shell and the bubbles is drawn to scale assuming that 
a burst with N = 15 is to be produced. Consequently N — 15 bubbles are drawn (more bubbles will add up to a smooth 
profile). The fraction of the shell that will impact these bubbles is small leading to high inefficiency. As N increases the 
efficiency problem becomes more severe ~ N^ 1 , from $0 
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Fig. 18. — The internal shock scenario. The source produces multiple shells as shown in this figures. The shells will have 
different Lorentz factors. Faster shells will catch up with slower ones and will collide, converting some of their kinetic 
energy to internal energy. This model is Type-II and naturally produces variable bursts, from Mw 
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Fig. 19. — A peak produced by a collision between two shells. The luminosity plotted versus the arrival time. The solid 
line corresponds to R = cSt e and the dotted line corresponds to R = 0, from 
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Fig. 20. — Luminosity vs. observer time, for different synthetic models a: jmin = 100, jmax = 1000, N = 100, r\ = —1 
and L/l = 5 b: lmm = 100, lmax = 1000, N = 100, rj = 1 and L/l = 5 c: 7mi „ = 100, lmax = 1000, N = 20 , V = -1 
and L/l — 5 d: -f, n i n — 100, 7 max = 1000, N = 20 , ry = — 1 and L/l = 1 e: -f rn i n = 100, j ma x — 1000, N — 100, random 
energy with E m n X = 1000 and L/l = 5 f: 7 m i„ = 100, 7 maa: = 1000, = 100, random density with p max = 1000 and 
L/l = 5. From 0- 
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Fig. 21. — The Lorentz factor j, the density p and the pressure p in the shocks. There are four regions: the ISM (region 
1), the shocked ISM (region 2), the shocked shell (region 3) and the unshocked shell (region 4), which are separated by the 
forward sh ock (FS), the contact discontinuity (CD) and the reverse shock (RS). The initial parameters are the same as in 
From /H|/. 
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fast cooling 




v(Hz) 

Fig. 22. — Synchrotron spectrum of a relativistic shock with a power-law electron distribution, (a) Fast cooling, which 
is expected at early times (t < to). The spectrum consists of four segments, identified as A, B, C, D. Self-absorption is 
important below v a . The frequencies, v m , v c , v a , decrease with time as indicated; the scalings above the arrows correspond 
to an adiabatic evolution, and the scalings below, in square brackets, to a fully radiative evolution, (b ) Slow cooling, which 
is ex pected at late times (t > to). The evolution is always adiabatic. The four segments are identified as E, F, G, H. From 
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Fig. 23. — Different scenarios in the A (in cm) - 7 plane for £ = 5/ A = 0.01. Relativistic ES occur for large A and large 
7 - upper right -above the £ = 1 /me (dark gray and light gray regions). Newtonian ES occur below £ = 1 - lower left - 
white region. IS occur, if there are sufficient variation in 7 below the £ = £ 2 / 3 line (light gray and white regions). The 
equal duration T = Is curve is shown for Newtonian ES (solid line) a relativistic ES (dotted line) and IS (dashed line). 
Note that a relativistic ES and an internal sh ock w ith the same parameters have the same overall duration T but different 
temporal substructure depending on S. From p23Bl. 



MO 




10 



10 



10 



10 



10 



10 



12 



Fig. 24. — Allowed regions for internal shocks in the 5 (in cm), 7 plane. Note that the horizontal 6 axis also corresponds 
to the typical peak duration, 5t multiplied by c. Internal shocks are impossible in the upper right (light gray) region. The 
lower boundary of this region depends on £ = 6/ A and are marked by two solid curves, the lower one for £ = 1 and the 
upper one for £ = 0.01. Also shown are r 77 = 1 for an observed spectrum with no upper bound (dotted line), r 77 = 1 
for an observed spectrum with an upper bound of lOOMeV (dashed line) and r e = 1 ( dashed- dotted) . The optically thin 
internal shock region is above the r = 1 curves and below the £ = £ 2 / 3 (solid) lines. From fftSU^ . 
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Fig. 25. — (a) Schematic description of the different radii for the case £ > 1. The different distances are marked on 
a logarithmic scale. Beginning from the inside we have AR, the initial size of the shell, R v , the radius in which a 
fireball becomes matter- dominated (see the following discussion), R c , the radius where inner shells overtake each other 
and collide, Ra, where the reverse shock reaches the inner boundary of the shell, and Rp, where the kinetic energy of the 
shell is converted into thermal energy, (b) Same as (a) for £ < 1. i?r does not appear here since it is not relevant. Rn 
marks the place where the reverse shock becomes relativistic. From 
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Fig. 26. — Calculated spectra or light curve from a Blandford-McKee solution. F u (t) is plotted as a function of <j) = 
Const, x vi?!" 1 . Thus if we consider a constant time this figure yields the spectrum while if we consider a fixed frequency 
it yields the light curve. The solid curve depictes emi ssion from the full fireball. The dashed line depictes the spectrum 
resultsing from emission along the line of sight. From [ 25i J. 
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Fig. 27. — Light curve due to synchrotron radiation from a spherical relativistic shock, ignoring the effect of self- absorption, 
(a) The highfrequency case (v > Vq). The light curve has four segments, separated by the critical times, t c , t m , i - The 
labels, B, C, D, H, indicate the correspondence with spectral segments in Fig. [13. The observed flux varies with time 
as indicated; the scalings within square brackets are for radiative evolution (wh ich i s restricted to t < to) and the other 
scalings are for adiabatic evolution, (b) The low frequency case (v < isq). From 124.6] - 
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Fig. 28. — Decline in the afterglow of GRB970228 in different wavelength and theoretical model, from 
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Fig. 29.— The X-ray to radio spectrum of GRB 970508 on May 21.0 UT (12.1 days after the event). The fit to the 
low-frequency part, <X4.86-86GHz = 0.44 i 0.07, is shown as well as the extrapolation from X-ray to optical (solid lines). 
The local optical slope (2.1-5.0 days after the event) is indicated by the thick solid line. Also indicated is the extrapol ation 
F v oc v~ o e (lines). Indicated are the rough estimates of the break frequencies v a , v m and v c for May 21.0 UT from \26(\f 



